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INTRODUCTION 

This  report  describes  the  results  of  a  six-month  analytic 
effort,  and  some  bench  tests  and  experiments,  conducted  to  deter¬ 
mine  the  areas  of  applicability,  and  the  limitations  of,  nonintrusive 
(noninterfering)  ultrasonic  techniques  to  measure  gas  mass  flow  rates. 

In  this  introductory  section  we  will  first  state  the  problem  in 
terms  of  AEDC's  requirements.  Then  we  will  identify  a  simple  duct 
as  a  starting  point  for  analysis  and  experiments.  We  will  then  indicate 
one  possible  extension  of  measurements  on  a  single  duct,  to  measure¬ 
ments  on  a  number  of  adjacent  ducts,  analogous  to  a  venturi  plate, 
but  with  minimal  pressure  drop.  Lastly,  we  will  mention  some  other 
applications  of  interest.  The  introductory  section  will  conclude  with 
a  table  that  tells  where  in  the  report  the  analyses  of  the  eleven 
contractually- specified  items  are  to  be  found  (see  page  10). 

The  rest  of  the  report  itself  is  organized  into  principal  sections 
(see  also,  Table  of  Contents)  covering  these  topics: 

Ultrasonic  Flowmeter  Measurement  Approaches 

Signal,  Noise,  Response  Time  and  Accuracy 

Transducers 

Other  Analytic  Considerations 
Experimental  Results 

Ultrasonic  Mass  Flowmeter  Design  for  Duct  or  Pipe, 

~250-  to  500-mm  Diameter,  Mach  No.  <  0.  5 

Conclus  ions 

Recommendations 

Acknowledgments 

References 

It  will  be  shown  that,  as  a  result  of  the  above  analysis  and  ex¬ 
periments,  we  are  able  to  adapt  one  of  our  ultrasonic  liquid  flowmetering 
designs  to  gas  flow,  suitable  for  tests  in  the  existing  Propulsion  Wind 
Tunnel  (PWT)  facilities  at  AEDC. 
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Table  1.  Location  in  report  of  analytic  considerations. 


Item  Page 

1.  Precision  and  bias  (uncertainty  <  0.  5%).  54,  55,  60 

2.  Response <  1  sec.  29,43,60,44 

3.  Stability:  system  must  continuously  conform 
to  a  particular  set  of  specifications  for  at  least 

24  hours.  60,  78-82 

4.  Ease  of  installation  and  operation.  54-61,66 

5.  Maintainability.  61 

6.  Calibration.  61 

7.  Correct  averaging  over  velocity  and  mass  flow 

profile.  14,25-29 

|  8.  Duct  material  and  precision  requirements.  44 

9.  Subsonic,  supersonic  and  hypersonic  velocities.  9-109 


10.  Effects  on  performance  caused  by: 


a. 

Noise  and  vibration  environment  as 
experienced  in  normal  wind  tunnels. 

30-51 

b. 

Particulate  scatterers  (e.  g.  ,  dust  in 
normal  wind  tunnels). 

17,32 

c. 

Humidity. 

32 

11.  Data  acquisition  and  recording.  62 


Let  us  now  state  the  problem,  in  language  essentially  the  same  as 
in  Annex  "A"  of  the  contract. 

I.  STATEMENT  OF  THE  PROBLEM 


In  day-to-day  aerodynamic  and  engine  testing  in  ground  test  facili¬ 
ties  one  of  the  most  common  requirements  is  for  the  determination  of  gas 
velocity  and  mass  flow  rates  in  ducts.  Conventional  methods  of  adequate 
accuracy,  e.  g.  ,  venturi  or  orifice  plates,  total  and  static  pressure  surveys 
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or  rakes,  hot  wire  or  hot  filament  probes,  etc.  ,  have  been  developed 
and  are  presently  used.  However,  these  conventional  methods  have 
disadvantages  in  certain  applications  which  at  present  are  accepted 
because  of  lack  of  a  better  solution.  To  list  a  few: 

1.  Venturi  and  orifice  plates  cause  a  sizable  pressure  loss  in  the 
flow,  and  require  undisturbed  inflow  conditions  and  space  which  may 
not  always  be  available. 

2.  Total  and  static  pressure  surveys  by  a  movable  single  probe 
are  impractical  for  routine  development  testing.  Total  and  static  fixed 
pressure  rakes  (multiple  probes)  disturb  the  flow,  normally  to  an  tin- 
known  degree. 

3.  Movable  hot  wire  or  filament  probe  surveys  are  too  cumbersome 
and  time  consuming  for  development  testing  in  large  facilities. 

In  recent  years  ultrasonic  methods  have  been  introduced  as  mea¬ 
suring  techniques  for  liquids  and  gases.  Nonintrusive  methods  have  been 
described  for  the  measurement  of  volume  and  mass  flow  rate,  tempera¬ 
tures,  speed  of  sound  and  density.  Flowmeters  have  been' developed  for 
liquids,  e.  g.  ,  water  or  fuel  flow  and  natural  gas  in  large  supply  lines. 
Although  there  has  been  progress  and  success  in  applying  the  ultrasonic 
technique  to  measure  flow  parameters,  these  developments  and  applica¬ 
tions  have  been  limited  and  actual  potential  of  acoustic  diagnostics  methods 
in  aerodynamics  testing  in  comparison  to  conventional  procedures  has  not 
been  fully  established. 

Some  of  the  expected  basic  advantages  of  ultrasonic  methods  are 
as  follows: 

1.  No  interference  with  flow  (transducers  recessed  or  on  outside 
wall). 

2.  Ease  of  installation  and  operation. 

3.  Fast  response. 

4.  Inherent  independence  of  gas  type  or  composition. 

5.  Minimum  space  requirement. 

6.  Potential  of  combining  measurements  of  several  flow  parameters 
such  as:  density,  temperature,  speed  of  sound,  and  volume  or 
mass  flow  rate. 
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The  purpose  of  the  present  effort  was  to  assess  the  feasibility, 
potential  practicability  and  performance  of  ultrasonic  techniques  to 
measure  velocity  and  mass  flow  rates  and  related  parameters  of  gases 
for  typical  applications  at  AEDC. 

II.  BACKGROUND 


Testing  in  the  Arnold  Engineering  Development  Center's  numer¬ 
ous  facilities  requires  in  many  instances  and  voider  different  conditions 
determination  of  gas  mass  flow  rates.  Present  techniques,  installations 
and  requirements  vary  considerably  over  large  ranges  of  mass  flows  to 
be  measured  at  permanently  installed  measuring  stations  or  at  individ¬ 
ually  designed  stations  inside  or  outside  a  test  model  for  specific 
inve  stigations . 

Examples  of  such  requirements  are: 


1.  Direct-Connect  Engine  Testing 
Duct  Mach  number:  0.  1  to  0.  3 


Duct  static  pressure: 
Temperature: 

Mass  flow  rate: 
Response  time: 

Duct  diameter: 
Typical  duct  length: 
Accuracy  desired: 
Present  technique: 


0.  5  to  400  psia 
-75  to  750°F 
200  to  2200  lbs/ sec 
1  sec 
5  to  30  ft 
30  ft 
0.  5% 

single  or  multiple  venturi  or 
orifice 


2.  Typical  Airplane  Model  with  Flow  Through  Duct 
Duct  Mach  number:  0.  2  to  0.  9 


Duct  static  pressure:  1.  5  to  15  psia 

Temperature:  100  to  130°F 
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Response  time:  1  sec 

Duct  diameter:  2  to  20  inches 

Length  of  duct:  5  to  20  diameters 

Accuracy  desired:  0.  5% 

Present  technique:  Total  and  static  pressure  rates 

3.  Supersonic  and  Hypersonic  Flow:  Other  needs  are  to  measure 
supersonic  and  hypersonic  velocities  without  disturbing  the  flow  in  test 
facilities  up  to  approximately  20,  000  feet  per  second.  Conventional 
probing  techniques  in  continuous  and  impulse  facilities  are  extremely 
difficult  for  these  conditions.  Not  all  these  conditions  will  occur  simul¬ 
taneously  and  are  given  below  only  to  indicate  the  range  of  conditions: 

Mach  number:  l.  5  to  20 

Velocity:  1500  to  20,  000  fps 

Temperature:  500  to  1 0,  000°R 

Pressure:  0.  1  to  15  psia 

Present  technique:  Total  and  static  pressure  rakes 

These  are  only  examples  and  general  application  at  AEDC  will 
require  other  conditions  of  velocity,  temperature,  pressure,  duct  Mach 
number,  etc.  ,  therefore,  these  are  not  presented  as  limits  to  the  region 
of  interest.  The  limitations  of  the  method  were  a  subject  of  the  study. 

Application  of  ultrasonic  measurements  for  conditions  approximately 
described  in  II.  2  was  of  immediate  interest  and  had  the  higher  priority, 
since  presently  no  completely  satisfactory  method  is  available.  Require¬ 
ments  similar  to  those  described  in  II.  1  may  develop  in  the  future  under 
conditions  where  venturis  or  orifices  may  no  longer  be  acceptable  (because 
of  pressure  loss). 

It  was,  therefore,  the  intent  of  the  proposed  effort  to  concentrate 
on  conditions  encompassing  the  requirements  of  II.  2. 
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m.  APPROACH 


The  contractor  was  to  analytically  evaluate  the  potential  and 
performance  of  different  nonintrusive  (noninterfering)  ultrasonic 
techniques  to  measure  gas  velocity  and  mass  flow  rates  and  related 
parameters  within  the  envelope  of  H.  2.  This  evaluation  could  extend 
beyond  the  described  envelope  to  determine  the  actual  boundaries 
within  which  the  essential  requirements  can  be  fulfilled.  The  initial 
studies  were  to  be  based  on  simple  straight  ducts  of  circular  or 
square  cross  section. 

Panametrics  was  to  describe  the  different  approaches  studied, 
presenting  advantages  and  disadvantages  of  each  and  recommend  the 
best  system  for  possible  experimental  verification.  ' 

The  analysis  was  to  cover,  but  not  be  limited,  to,  the  items  in 
Table  1.  During  the  analytical  phase  we  decided  to  conduct  certain 
small  scale  bench  tests  of  system  components  as  required  to  supple¬ 
ment  and  support  the  analytical  study. 

If  in  the  analysis  a  system  was  defined  which  could  be  expected 
to  meet  the  requirements  in  II.  2,  the  contractor  was  to  briefly  explore  the 
possible  extension  of  such  system  to  the  more  severe  conditions  in  II.  1 
and  H.  3.  Regarding  this  cited  Requirement  2,  it  was  agreed  early  in 
the  program,  during  a  visit  to  AEDC,  to  conduct  "bench  tests"  on  a 
simple  duct  under  these  conditions  including  as  much  as  possible  of  the 
ranges  in  Requirement  2.  This  turned  out  to  involve  design,  construc¬ 
tion  and  testing  in  a  square  duct  ~25  cm  x  25  cm,  operated  essentially 
at  atmospheric  temperature  and  pressure,  at  a  Mach  No.  of  ~0.  15. 

During  this  same  visit  to  AEDC,  it  was  emphasized  by  AEDC  and 
ARO  personnel,  that  the  pressure  drop  associated  with  existing  venturi 
plates  was  undesirable  from  both  technical  and  energy- consumption  points 
of  view,  especially  in  large-diameter  ducts.  Some  of  these  ducts  are  sev¬ 
eral  meters  in  diameter,  which  may  be  too  large  to  reliably  propagate 
across,  with  adequate  signal/ noise  to  achieve  the  desired  accuracy  of  1/2%. 
For  such  large  ducts,  we  suggested  inserting  an  assembly  of  short  tubes, 
or  a  honeycomb,  to  divide  the  area  into  cells  of  diameter  small  enough  to 
be  measured  ultrasonic  ally.  The  cell  walls  could  be  rather  thin,  ~1  mm, 
plus  some  allowance  for  transducer  bodies  and  cabling.  It  would  therefore 
appear  that  even  a  large  number  of  cells,  ~100,  might  be  installed,  and 
yet  avoid  pressure  drops  larger  than  a  few  percent.  Admittedly,  this 
"ultrasonic  honeycomb"  concept  is  not  nonintrusive.  But  it  appears  to 
offer  a  profiling  and/or  averaging  capability,  with  minimal  penalty  due 
to  pressure  drop.  It  also  offers  the  prospect  of  applying  a  design  which 
can  be  tested  and  verified  in  a  small  flow  facility,  to  much  larger  ducts, 
simply  by  connecting  additional  cells,  to  map  out  the  full  duct  area  of  interest. 
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ULTRASONIC  FLOWMETER  MEASUREMENT  APPROACHES 

Determination  of  the  total  mass  flow  rate  M  through  a  duct  gen¬ 
erally  involves  integrating  the  local  p  v  products.  Let  us  first  consider 
ultrasonic  principles  related  to  measuring  v,  and  then  ultrasonic  and 
other  principles  related  to  p  . 

ULTRASONIC  PRINCIPLES  UNDERLYING  v  DETERMINATION 


We  may  distinguish  six  (6)  different  ways  of  ultra sonically  sens¬ 
ing  or  interacting  with  a  flowing  gas: 


•  Transmission 

•  Reflection/ Doppler  Scatter 

•  Beam  Drift 

•  Correlation 

•  Noise 

•  Vortex  Shedding 


Principle  Summarized  in  Equation  Form 
Af  ~  v;  At  ~  v/c^;  At/ Et  ~  v/c 


Af  ~ 

v;  At 

fD~ 

1  v/c 

x/D 

=  v/c 

x/r 

=  V 

n  ~ 

vP,  p  1 

vP,  p  «  3  to  8,  depending  on  v 
f  ~  v 


Nonintrusive  aspects  of  the  first  four  of  these  have  been  discussed 
in  the  paper  (Lynnworth  et  al.  ,  1971)  presented  in  May  1971  at  the  ISA 
Symposium  on  Flow,  Noninvasive  Ultrasonic  Measurement  of  Flow  Velocity 
and  Mass  Flow  Rate.  With  the  permission  of  the  ISA,  a  copy  of  this  paper 
is  appended. 

The  principle  of  passively  "listening"  to  the  flow  rate  is  the  basis 
for  at  least  one  commercially  available  ultrasonic  flowmeter  (abstract  349, 
p.  26,  in  Dowden,  1972).  The  principle  of  measuring  ultras onically,  the 
frequency  at  which  vortices  are  shed  by  a  strut,  is  the  basis  for  other  com¬ 
mercially  available  ultrasonic  flowmeters  for  liquids  and  gases.  However, 
the  intrusivenesB  of  the  strut  generally  puts  vortex  shedding  devices  beyond 
the  scope  of  extensive  consideration  in  this  report.  Likewise,  the  intrusive¬ 
ness  of  fluidic  flowmeters  (abstracts  1672-1676,  p.  140,  in  Dowden)  rule  oul 
further  consideration  of  this  type  flowmeter  in  this  report,  even  though  the 
fluidic  oscillator  frequency  may  lie  in  the  ultrasonic  range. 


For  general  background  and  historical  perspectives,  interesting 
reviews  of  sound  propagation  in  air  are  given  by  Allen  et  al.  (1950)  and 
later  by  Nyborg  and  Mintzer  (1955).  Sonic  anemometry  was  reviewed 
briefly  by  Foxx,  1968.  Let  us  now  briefly  identify  the  six  above-listed  ap¬ 
proaches.  We  will  then  concentrate  on  those  principles  most  relevant  to 
AEDC's  Requirements. 
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T  ransmission 


Principle:  Sound  travels  faster  with  the  wind  than  against  it. 


Typical  Configurations: 


■FLOW 


Fig.  la.  Transmission  approach  wherein  upstream  and  downstream 
waves  propagate  over  the  same  path. 


FLOW 

...  ■ 

JL-L 

c\ 


TRANSDUCER 


SOUND  SOURCE 


TRANSDUCER 


RECTANGULAR  DUCT 
19  x  21  mm 
(0.75x0.825  in.) 


Fig.  lb.  Transmission  approach  wherein  a  single  sound  source 
launches  upstream  and  downstream  waves  which  propagate  over  dif¬ 
ferent  paths  (Ingaard  and  Singhal,  1973). 


Mathematics:  If  time  intervals  t^  and  t ^  are  measured  upstream  and 
downstream,  their  difference  is  At  =  2  l  v/c2.  A  similar  relation  holds, 
if  phase  shifts  are  measured.  If  reciprocals  of  time  or  phase  are  mea¬ 
sured,  their  difference  Af  is  directly  proportional  to  v,  independent  of  c. 
See  Figs,  la,  lb. 

Panametrics'  Application  in  Liquids:  Aviation  fuel  mass  flow  rate  M 
(Eustis  Directorate  programs;  reference:  Lynnworth  et  al.  ,  1972); 
v  in  water  and  cryogenic  liquids  (NASA-Langley  programs;  reference: 
Pedersen  et  al.  ,  1972,  1973).  Summaries  of  these  transmission  ap¬ 
proaches  are  appended.  Technique  is  widely  used  in  commercial  ultra¬ 
sonic  flowmeters  for  liquids  and  gases,  is  seen  in  ultrasonic  flowmeters 
currently  manufactured  by  NUSonics,  Saratoga  Systems,  Scarpa  Labora¬ 
tories,  and  Westinghouse.  This  principle,  used  at  MIT  to  a  Mach  No.  of 
~0.  5,  in  a  small  duct,  was  reported  recently,  and  is  sketched  in  Fig.  lb 
(Ingaard  and  Singhal,  1973).  In  this  MIT  work,  the  sound  source  and  re¬ 
ceiving  transducers  were  isolated  from  the  duct  wall,  and  were  coupled 
directly  to  the  air  flowing  in  the  duct. 
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Principle:  The  frequency  of  a  wave  is  shifted  when  the  wave  is  reflected 
from  a  moving  scatterer. 


Typical  Configuration: 

Fig.  2.  In  doppler  flowmeters,  the  trans¬ 
mitted  and  scattered  beams  of  interest 
may  be  orthogonal,  as  shown  here,  or 
oblique,  or  parallel  to  one  another.  As 
long  as  at  least  one  of  the  beams  propa¬ 
gates  not  normal  to  the  flow  vector,  the 
scattered  wave  will  be  doppler- shifted. 


Mathematics:  The  difference  f^  between  the  transmitted  frequency  fQ  and 
the  received  frequency  f]_,  is  proportional  to  v/c.  The  "constant"  of  pro¬ 
portionality  depends  partly  on  geometry.  It  ranges  between  0  and  2  fQ. 

This  approach  is  mainly  used  in  biomedical  studies  of  blood  flow 
and  tissue  movement  (Ultrasonics  11  (2),  52-53,  March,  1973).  In  indus¬ 
trial  applications,  this  approach  has  been  considered  when  the  fluid  already 
contains  scatterers.  Examples  include  contaminants  in  aviation  fuel,  vapor 
bubbles  in  prepolymers,  and  the  as-yet-undetermined  scatterers  in  rela¬ 
tively  pure  liquid  metals  (Fowlis,  1 9 7 3 ) . The  doppler  method  offers  a  the¬ 
oretical  immunity  to  the  acoustic  short  circuit  problem  which  besets  most 
transmission  approaches.  However,  it  suffers  from  phase  shifts  due  to 
temperature  gradients.  It  provides  correct  area-averaging  only  under  con¬ 
ditions  of  uniform  intensity  and  uniform  distribution  of  scatterers,  which 
conditions  would  be  rather  difficult  to  maintain  in  flowing  gases.  Lastly, 
in  order  to  equate  the  motion  of  scatterers  with  the  motion  of  the  gas,  one 
must  solve  the  dilemma  of  having  a  particle  small  enough  to  travel  precisely 
with  the  gas,  yet  large  enough  to  scatter  efficiently.  Such  particles  gener¬ 
ally  are  not  to  be  naturally  expected  in  a  flow  facility,  nor  is  their  injection 
into  the  gas  considered  welcome  or  practical.  Hence,  we  terminate  con¬ 
sideration  of  the  doppler  approach. 

*Lynnworth  has  suggested  that  strong  velocity  gradients,  or  turbulence 
might  be  a  source  of  backs catte ring  in  an  otherwise  homogeneous  uniform- 
density  liquid.  Apart  from  the  well-known  ray-bending  or  forward¬ 
scattering  mechanisms,  small  vortices  might  cause  local,  perhaps  short¬ 
term  departures  from  c,  in  the  effective  sound  propagation  speed.  Thus, 
backs  catte  ring  could  result  due  to  local  pc  mismatch.  This  is  somewhat 
analogous  to  backs cattering  from  anisotropic  grains  in  a  polycrystalline 
solid.  (Earlier  (February  1973)  Pedersen  identified  the  short-term  radial 
components  of  vortices  as  having  a  significant  effect  on  dynamic,  fast- 
response  ultrasonic  flow  measurements  where  the  acoustic  path  includes 
a  radial  component  of  the  eddies.  )  See  also:  D.  W.  Beran  et  al.  ,  J.  Acoust. 
Soc.  Amer.  55,  2,  334-338  (Feb.  1974). 
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Beam  Drift 


Principle:  Sound  beam  is  carried  downstream  by  wind. 
Typical  Configuration: 


NO  FLOW 

FLOW 


Fig.  3.  The  beam  drift  approach  to  measuring  gas  flow  velocity  is 
probably  the  simplest  ultrasonic  flowmetering  method  in  use.  It  is 
applicable,  depending  on  geometry, 'to  flows  below  Mach  No.  0.1,  up 
to  Mach  1  and  beyond.  Not  unexpectedly,  however,  the  advantages  of 
simplicity  and  rangeability  are  offset  by  accuracy  limitations,  esti¬ 
mated  as  ~5%,  in  one  of  the  latest  systems  under  evaluation. 

Mathematics:  If  the  normal  to  the  wave  fronts  is  orthogonal  to  the  flow 
vector,  the  drift  angle  0  is  given  by  tan  0  =  v/c.  The  downstream  drift  x 
amounts  to  D  v/c.  The  amplitude  and  phase  of  the  wave  at  a  particular 
receiving  point  will  depend  on  v/c’ over  the  path. 

Application  in  Air:  National  Physical  Laboratory,  England.  Reported 
in  journal  Ultrasonics,  pp.  197-199,  September  1972.  Flow  was  mea¬ 
sured  at  1  MHz  in  ducts  up  to  1  ft  diameter.  Panametrics  used  a  beam 
drift  method  to  measure  the  Mach  No.  of  combustor  exhaust  gas  (M  «  .  25) 
in  a  program  conducted  for  GE-Evendale,  Ohio  (jet  engine  turbine  inlet 
simulator).  We  have  learned  from  G,  G.  Twidle  at  NPL,  that  in  the  beam 
drift  approach  under  investigation  in  England's  National  Physical  Labora¬ 
tory  (NPL),  the  present  NPL  system  uses  discrete  transducers,  not  a 
continuous  transducer.  The  discrete  transducers  would  be  expected  to 
yield  accuracy  of  '"’5%,  based  on  processing  of  amplitude  information  only, 
as  a  measure  of  the  beam  drift,  and  hence  a  measure  of  the  Mach  No.  The 
NPL  system  does  not  process  phase  information,  which  if  used,  could  in 
principle  improve  the  accuracy.  Delivery  of  the  final  components  for  the 
NPL  Flowmeter  was  scheduled  to  be  completed  in  1973.  At  the  present 
time,  no  test  results  have  been  made  available  to  us.  At  the  present  time, 
therefore,  we  have  no  evidence  that  this  technique,  as  practiced,  can  yield 
accuracy  better  than  5%.  Whether  an  order- of-magnitude  improvement  is 
possible  based  on  phase  measurements  instead  of  or  in  addition  to  ampli¬ 
tude  measurements,  is  not  yet  known. 
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Correlation 


Principle:  Tags  "attached"  to  the  flowing  air  will  successively  modulate 
two  interrogating  beams  or  two  recordings  from  spaced  microphones  at  a 
time  interval  T  inversely  proportional  to  v,  and  essentially  independent 
of  c. 


Fig.  4b. 


PULSE-  ECHO 


THROUGH - 
TRANSMISSION 


Fig.  4.  Correlation  approach  to  measuring  gas  flow  velocity.  In  (a)  a 
tag,  which  may  be  a  naturally- occurring  eddy,  or  an  introduced  bubble,  is 
indicated  as  a  possible  means  of  providing  information  on  flow  velocity  v 
along  the  surface  of  a  model.  In  (b)  a  pulse- echo  version  is  shown,  other¬ 
wise  similar  to  the  through-transmission  version  practiced  by  Coulthard 
(1972)  and  shown  in  the  right-half  of  (c).  At  the  left-half  of  (c)  we  also 
show  an  optional  ultrasonic  source  of  pulsed  high-intensity  waves,  to  in¬ 
troduce  a  tag  nonintrusively. 
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Mathematics:  The  correlation  time  T  =  x/v  where  x  =  distance  between 
beams. 

Application:  Normally-occurring  thermal  "tags"  have  been  exploited  in 
order  to  measure  flow  velocity  profile  in  liquid  metals  by  Wheelock  et  al. 
of  the  Liquid  Metals  Engineering  Center  in  California.  Thermocouple 
rakes  detected  the  thermal  inhomogeneities.  We  and  others,  have  sug¬ 
gested  using  nonintrusive  ultrasonic  beams  to  perform  an  analogous 
ultrasonic  measurement  (Fig.  4). 

In  Fig.  4a  we  sketch  a  possible  method  for  remotely  determining 
the  value  of  v  next  to  a  model  surface.  In  Fig.  4b  an  average  flow  is 
sensed.  Neglecting  attenuation,  it  is  interesting  to  note  that  if  the  trans¬ 
ducers  in  Fig.  4b  are  applied  circumferentially  so  as  to  interrogate  the 
entire  area  within  the  (assumed)  circular  duct  or  pipe,  a  measure  of  the 
area-averaged  flow  velocity  could  be  obtained,  independent  of  flow  pro¬ 
file.  On  the  contrary,  the  same  transducer  configuration,  utilized  in  a 
transmission  measurement,  would  not  generally  yield  the  area-average, 
but  the  diameter-average."' 

Figure  4c  represents  a  cross -correlation  flowmeter  developed  at 
Tees  side  Polytechnic  in  England  (Coulthard,  1973).  Consider  a  disturb¬ 
ance  detected  at  A  and  then  later  at  B.  If  the  signal  recorded  at  point  A 
is  x(t)  whose  value  is  directly  related  to  the  instantaneous  flow  at  point  B, 
then  the  cross -correlation  function  expressing  x(t)  and  y(t)  in  terms  of  a 
time  variable  T  is  given  by  the  expression 

T 

( T)  =  lim  ^  f  x(t )  y(t  +  T)  dt. 

*Y  T— oo  1  J0 

If  all  components  of  the  flow  can  be  assumed  to  move  at  a  velocity 
v,  then,  for  a  transducer  spacing  x  the  flow  velocity  v  is  given  by  the 
expression 

x 

V  =  T 
m 

where  Tm  is  the  value  of  T  corresponding  to  the  peak  value  of  ^Xy(T). 

It  may  be  argued  that  the  value  of  T  corresponding  to  the  center 
of  gravity  of  the  cross -cor  relation  peak  should  be  taken  as  a  measure  of 

^Optionally,  an  upstream  circumferential  transducer  periodically  could 
be  energized  impulsively  to  yield  a  brief,  high- intensity  beam  directed 
radially  inward  (Fig.  4c).  This  would  nonintrusively  produce  a  thermal 
p\ils e  or  tag  in  the  gas.  As  the  tag  was  carried  downstream,  its  velocity 
could  be  measured  (see  also  Fig.  6,  p.  24). 
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the  time  of  flight  of  the  fluid  between  A  and  B,  particularly  if  the  peak 
exhibits  any  degree  of  asymmetry.  Coulthard,  however,  took  the  peak 
value  of  tf'xytT)  as  the  measure  of  flow  since  this  lends  itself  more 
readily  to  electronic  detection  in  a  practical  flowmeter.  Departures 
from  linearity  achieved  in  initial  trials  were  reported  to  be  ~5%,  for 
these  experimental  conditions;  air  path,  150  mm  (6  in.  );  downstream 
separation  between  transducers,  250  mm  (10  in.  );  v  =  10  to  35  m/s. 
Further  consideration  of  the  correlation  flowmeter  is  postponed  until 
such  time  as  theory  predicts  or  experiments  confirm  an  accuracy  of 
~0.  5  to  1%. 

Noise 


In  a  technical  historical  paper  on  advances  in  aero-acoustics 
(Powell,  1965),  it  was  shown  that  several  methods  of  expressing  aero- 
dynamically  generated  sound,  in  the  absence  of  solid  surfaces,  yielded 
the  "eighth  power  law,  "  i.  e.  ,  the  sound  power  varies  as  the  eighth 
power  of  some  characteristic  flow  velocity,  on  corresponding  assump¬ 
tions  of  similarity.  In  this  program,  it  is  recognized  that  the  existence 
of  a  law  for  sound  power  might  provide  a  basis  for  measuring  flow  ve¬ 
locity.  Initially,  however,  the  principal  subject  for  application  of  such 
a  theory  was  the  noise  caused  by  the  exhaust  flow  of  jet  engines.  The 
eighth  power  relationship  was  found,  but  the  sound  field  was  heavily 
biased  at  an  acute  angle  to  the  jet  direction.  Fig.  5a.  A  transformation 
to  moving  axes  explained  how  convection  effects  produced  a  bias  of  the 
observed  type,  but  this  simultaneously  caused  the  power  to  increase 
more  rapidly  than  according  to  the  aforementioned  eighth  power  law  re¬ 
lationship  and  markedly  so  near  a  convection  speed  approaching  the  speed 
of  sound  in  the  atmosphere.  It  is  noted  that  experimental  data  showed 
that  this  relationship  held  remarkably  well  to  beyond  such  a  speed,  as 
respectively  indicated  in  the  uppermost  and  lowest  curves  of  Fig.  5b,  in 
which  the  convection  velocity  is  supposed  to  be  well  represented  by  half 
the  jet  exit  velocity. 

As  Powell  explains,  if  one  accepted  the  well-established  fact  that 
the  simple  dimensional  analysis  was  consistent  with  experimental  data  so 
far  as  the  total  sound  power  is  concerned  (and  took  the  view  that  the  con¬ 
vection  effects  were  yet  to  be  fully  explained),  then  it  seemed  plausible 
that  various  given  volumes  of  the  jet  would,  more  or  less,  also  follow  the 
relationships  of  the  simple  similarity.  Thus  by  assuming  the  flow  charac¬ 
teristics  to  be  similar  in  various  slices  of  the  jet,  e.  g.  ,  that  the  Strouhal 
numbers  of  the  turbulence  spectral  peaks  were  all  the  same,  then  one  could 
discuss  the  acoustic  power  produced  per  unit  length  of  the  jet.  Like  a 
boundary  layer,  the  acoustic  power  was  determined  to  depend  only  on  the 
area  of  the  constant- shear  layer,  because  as  the  mixing  region  increases 
in  thickness,  the  lesser  gradient  of  the  shear  and  the  lesser  frequency 
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turn  out  to  be  just  compensated  by  the  larger  eddy  size  and  thickness 
of  the  mixing  region  occupied  by  them.  However,  further  downstream, 
where  the  mean  velocity  on  the  jet  axis  falls  off  rather  rapidly,  the 
assumed  corresponding  reduction  in  turbulent  velocity  and  frequency 
fluctuations  completely  overwhelm  all  other  factors,  so  that  the  sound 
power  generated  per  unit  length  then  falls  off  with  the  seventh  power 
of  the  distance  from  the  orifice. 

Continuing,  Powell  traces  the  later  improvement  of  the  convec¬ 
tion  transformation  itself.  There  were  two  quite  independent  approaches 
which  yielded  very  similar  forms  of  results.  These  show  how  the  eighth 
power  law  ultimately  gives  way  to  a  lower  one  at  supersonic  speeds,  as 
indeed  it  must  sooner  or  later  since  the  total  acoustic  power  can  never 
exceed  the  total  jet  power,  in  practice  being  about  one  percent  or  so  of 
it.  No  convection  limitations  {and  consequently  no  diffractive  effects)  are 
incorporated  in  these  theories,  and  this  possibly  goes  towards  accounting 
for  the  fact  that  the  theoretical  curve  still  rises  appreciably  beyond  the 
simple  law,  as  the  middle  curve  of  Fig.  5b  suggests.  At  high  speeds  the 
theories  suggest  a  third  power  law,  if  the  effects  of  the  change  in  jet 
structure  and  density  are  neglected  (Fig.  5c). 

These  later  developments  have  gone  far  to  provide  explanations 
of  the  various  features  of  jet  noise,  about  which  much  of  the  interest  in 
aerodynamic  noise  has  been. 

At  present,  one  type  of  acoustic  noise  flowmeter  is  available 
commercially,  manufactured  by  Scarpa  Laboratories,  for  nonintrusive 
use  on  pipes  in  which  a  liquid  is  flowing.  Accuracy  is  claimed  to  be  2%. 
However,  to  relate  noise  due  to  a  flowing  gas,  to  the  flow  velocity,  in¬ 
dependent  of  the  local  solid  surfaces,  appears  to  be  a  rather  formidable 
task* 


Consideration  of  a  noise  flowmeter  for  gases  appears  to  be  of 
academic  interest,  and  perhaps  of  limited  experimental  interest,  where, 
for  gas  flow,  v  estimates  only  somewhat  better  than  order-of-magnitude 
would  be  adequate.  The  possibility  of  measuring  v  to  ~1%  appears  remote, 
based  on  the  nature  of  the  variables,  as  summarized  by  Powell.  Accord¬ 
ingly,  the  noise  flowmeter  approach  will  not  be  pursued  further. 

*See,  however,  abstract  349,  p.  26,  in  Dowden,  1972.  It  was  reported 
in  1969  that  the  power  level  of  the  ultrasonic  noise  was  directly  propor¬ 
tional  to  the  volume  flow  rate  of  the  moving  gas,  liquid,  slurry  or  powder 
(New  Scientist  44  (679),  558,  December  11,  1969).  Theoretically,  of 
course,  it  is  clear  that  "direct  proportionality”  could  be  expected  as  an 
approximation  only  over  a  limited  flow  range.  A  flow/no  flow  acoustic 
noise  sensor,  model  SF-14A,  is  manufactured  by  Sensor  Technology 
Company  of  Seattle. 
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Fig.  5a.  Direction  distribution  of  the  far  field  sound  radiation  for  turbulent 
jets  of  three  different  speeds. 


CONVECTION  MACH  NUMBER 


Fig.  5b.  The  straight  line  is  for  the  "eighth  power  law".  The  highest 
curve  is  that  when  the  original  convection  theory  was  applied,  middle  curve 
being  for  the  improved  convection  theory.  The  convection  curves  are  approxi 
mations  only,  based  upon  the  average  effect  for  various  types  of  quadrupoles. 


Fig.  5c.  The  sound  power  "P"  of  various  classes  of  jets  after  the  effect 
of  the  diameter  "D"  has  been  allowed. 
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Vortex-Shedding 

Principle:  The  rate  or  frequency  at  which  vortices  or  eddies  are  shed 
from  a  body  effectively  depends,  in  the  simplest  case,  on  flow  velocity 
v.  The  shed  vortices  can  modulate  an  ultrasonic  beam  aimed  across 
the  wake. 

Typical  Configuration: 


Fig.  6.  Vortex- shedding  ultrasonic  flowmeter. 


Mathematics:  f  ~  v.  The  mathematics  is  more  complicated  for  non¬ 
circular  bodies  (Mujumdar  et  al.  ,  1970)  for  turbulent  flow  and  for  flow 
gradients. 

Application:  Airspeed,  angle  of  attack  measurements.  These  and  other 
applications  are  being  developed  by  J- Tec  Associates  (Klass,  1972).  The 
basic  technique  is  described  in  U.  S.  Patent  No.  3,680,  375.  Accuracy 
is  1%  of  full-scale. 

Comment:  This  principle  is  considered  generally  not  applicable  to  the 
present  requirements,  because: 

•  Body  which  sheds  vortices  is  inherently  intrusive  (exception:  shedding, 
off  parts  of  model  in  duct). 

•  Shedding  frequency  depends  on  average  Reynolds  number,  or  average 
v,  over  an  extended  region.  However,  a  strut  of  nonuniform  diameter, 
or  a  group  of  parallel  struts,  might  be  used  to  obtain  profile  information. 
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PRINCIPLES  UNDERLYING  DENSITY  DETERMINATION 
Ultrasonic  p  orpc  Determinations 

There  are  several  ultrasonic  approaches  to  measuring  fluid 
density  or  impedance: 

•  Measure  resonant  frequency  of  vane  (November,  1972); 
commercially  available  from  ITT-Barton.  Yields  p  near  vane. 

•  Measure  electrical  input  impedance  of  transducer  (e.  g.  ,  quartz), 
radiating  into  fluid  (Kritz,  1955).  Yields  p  c  near  transducer. 

•  Compare  reflection  from  rod  whose  end  is  in  contact  with  fluid,  to 
reflection  from  another  portion  of  rod,  or  from  another  similar 
rod,  which  is  not  radiating  into  the  fluid  (Lynnworth  et  al.  ,  1971, 

1972).  Yields  p  c  at  rod  end. 

•  Measure  sound  speed  in  fluid.  Yields  average  c  or  c  profile  across 
duct. 

•  Measure  signal  amplitude  after  transmission  across  duct  or  across 
other  sampled  volume.  Result  depends  on  p  c  and  attenuation,  not 
just  p . 

If  the  gas  density  is  uniform  or  of  predictable  profile,  nonultra- 
sonic  methods  are  probably  easiest  to  use.  If  the  profile  is  not  predictable, 
then  the  resonant  vane  is  ruled  out  by  virtue  of  the  intrusiveness  required 
to  determine  the  profile,  and  the  next  two  "loading"  methods  are  ruled  out 
because  they  sample  only  the  fluid  close  to  the  loaded  or  radiating  surface. 
Although  they  appear  to  be  moderately  (<,  0.  3%)  sensitive  for  detecting 
liquid  p  c,  they  would  be  essentially  insensitive  to  changes  in  gas  pc  over 
the  conditions  of  this  program.  The  last  method  is  not  likely  to  provide 
high  accuracy  because  of  amplitude  fluctuations  due  to  turbulence  and 
attenuation.  Thus,  the  only  nonintrusive  ultrasonic  densitometry  approach 
holding  potential  merit  is  based  on  measuring  the  speed  of  sound  c  in  the 
gas. 


One  application  of  c  to  p  profile  determination  is  as  follows.  To 
the  extent  that  the  gas  is  isobaric,  the  density  profile  is  determined  by 
the  temperature  profile,  and  a  knowledge  of  the  static  pressure  at  one 
point,  say,  at  the  wall.  Now  the  temperature  (T)  profile  can  be  deter¬ 
mined  from  multipath  measurements  of  c  (which  yield  the  c  profile).  That 
is,  when  p  =  constant,  c — -T — ►  p  . 
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A  check  on  the  isobaric  state  may  be  obtainable  by  using  a  mo¬ 
mentum  equation.  This  would  utilize  the  v  profile,  possibly  measured 
with  the  same  transducers  used  to  get  the  c  profile. 

It  should  be  noted  that  ultrasonic  impedance-sensing  probes, 
which  yield  pc,  rather  than  p  alone,  are  not  necessarily  undesirable 
for  that  reason.  For  it  is  possible  to  combine  transmission  measure¬ 
ments  upstream  and  downstream,  or  across  the  stream  and  oblique  or 
along  the  stream,  to  obtain  the  Mach  No.  v/c,  which,  when  multiplied 
by  p  c,  results  in  a  product  proportional  to  M,  the  mass  flow  rate. 

Given  the  nonintrusive  constraint,  the  only  instrumentation 
(for  p  )  techniques  in  common  use  today  involve  nuclear  measurements. 
Recognizing  that  objections  have  been  raised  in  the  past  about  this  ap¬ 
proach,  and  following  AEDC's  recommendation  not  to  pursue  this  approach 
within  the  present  contract,  we  have  no  new  results  to  report  in  this  area. 
However,  to  provide  some  perspective  on  the  present  work,  and  to  provide 
a  basis  for  discussion  concerning  possible  future  requirements  for  densi¬ 
tometry  not  adequately  covered  by  existing  nonnuclear  technologies,  we 
include  next,  a  brief  discussion  of  nuclear  measurements  of  air  density 
profile,  prepared  by  F.  B.  Sellers. 

NUCLEAR  MEASUREMENT  OF  AIR  DENSITY  PROFILE 


a.  Two  Dimensional  Flow 

The  measurement  of  density  profiles  in  the  vicinity  of  an  object, 
such  as  a  wing  section,  which  produces  essentially  a  two  dimensional  flow 
pattern  is  straightforward.  In  such  instances  the  symmetry  of  the  flow 
allows  measurement  of  the  attenuation  of  a  beam  of  nuclear  radiation  (x  or 
beta  rays  for  densities  of  interest  here)  to  define  the  localized  density.  The 
quantity  measured  in  this  type  gage  is  the  integral  of  the  density  p  (g/crn^) 
along  the  radiation  path  x(cm) 

f 

d(g/cm2)  =  f  p(x)dx  (1) 

*/0 

from  the  source  location  at  x  =  0  to  the  detector  location  at  1  .  The  aver¬ 
age  density  is  then 


1 


If  p  (x)  is  constant,  independent  of  x,  then  p  =  p  .  which  is  the  case  for 
two  dimensional  flow.  This  is  because  the  beam  can  be  directed  such 

*Various  electromagnetic  techniques,  which  are  based  on  index  or  di¬ 
electric  constant  changes  with  p  ,  are  not  considered  in  this  report. 
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that  it  is  perpendicular  to  the  flow  pattern  at  any  desired  location  in 
the  pattern,  and  the  density  at  every  point  within  the  beam  path  is  the 
same. 


Individual  measurements  with  a  density  accuracy  of  better  than 
1%  with  a  time  constant  on  the  order  of  a  second  are  achievable  for  the 
attenuation  method. 

b.  Three  Dimensional  Flow 

Three  dimensional  flow  patterns  present  a  more  difficult  prob¬ 
lem,  because  the  nuclear  radiation  beam  cannot  normally  be  directed 
such  that  the  density  at  all  points  along  its  path  is  the  same.  An  atten¬ 
uation  measurement  can  still  be  employed  usefully  in  those  instances  in 
which  only  the  total  mass  per  unit  area  d(g/cm2)  along  the  path  of  the 
radiation  is  of  interest,  but  in  this  case  the  average  density  p”  defined 
from  (2)  is  not  equal  to  the  actual  density  along  the  radiation  path  (except 
by  chance).  However,  measurement  of  the  density  within  an  arbitrary 
localized  volume  is  not  possible  by  attenuation. 

In  such  a  case  it  is  possible  to  use  a  scattering  technique,  which 
is  illustrated  by  the  simple  drawing  in  Fig.  7.  Here  a  collimated  radia¬ 
tion  source  produces  a  beam  of  radiation  which  passes  through  the  skin 
(or  a  small  aperture)  of  the  body  to  be  measured,  and  a  detector  is  col¬ 
limated  in  such  a  manner  that  it  only  observes  radiation  scattered  from 
a  particular,  "sensed",  volume.  In  this  case  the  detector  output  depends 
on  the  density  within  the  localized,  or  sensed,  volume.  If  the  radiation 
is  such  that  it  is  attenuated  along  the  path  to  and  from  the  volume,  then 
the  detector  output  also  depends  on  the  quantity  d  defined  in  (1).  But  if 
the  density  profile  is  measured  starting  at  a  point  sufficiently  close  to 
the  body,  then  the  values  of  density  needed  for  measurements  made  far 
from  the  body  are  available. 

While  the  source  and  detector  are  shown  inside  the  body  in  Fig.  7, 
they  can  also  be  located  outside  if  measurements  are  to  be  made  far  from 
the  body,  or  if  the  body  is  too  small  to  contain  them.  In  any  case,  it  is 
clearly  necessary  to  move  either  the  source,  detector,  or  both  simulta¬ 
neously  to  complete  a  scan  to  determine  a  density  profile. 

A  technique  that  may  be  of  some  use  for  this  method  is  introduc¬ 
tion  into  the  flow  pattern  of  powdered  particles  of  materials  (such  as 
metals)  whose  scattering  cross  section  is  larger  than  that  of  air.  Pro¬ 
vided  the  density  ratio  of  added  material  to  air  is  independent  of  position 

*The  above  mentioned  objection  to  seeding  still  would  need  to  be  resolved. 
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Fig.  7.  Geometry  for  determination  of  density  profiles  by  nuclear 
radiation  scattering. 
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within  the  flow  pattern,  this  can  be  of  significant  aid  in  increasing 
the  detector  response  for  a  given  radiation  source.  The  detector 
response  could  also  be  enhanced  if  the  material  chosen  were  such 
that  its  fluorescent  x-ray  were  excited  by  the  source.  This  would 
be  the  case,  for  instance,  if  a  Cd-109  22  keV  x-ray  source  were 
used  and  molybdenum,  whose  K  absorption  edge  is  at  20  keV,  were 
used  as  the  powdered  addition. 

It  would  be  necessary  during  any  program  to  consider  the 
maximum  radioactivity  source  that  could  be  used  under  these  con¬ 
ditions.  However,  provided  the  source  is  sufficiently  intense  it  is 
possible  that  accuracy  better  than  1%  can  be  achieved  with  spatial 
resolution  of  less  than  1  cm^  with  time  constants  of  a  few  seconds. 
Theoretically,  the  number  of  counts  that  must  be  recorded  in  a  gage 
of  this  type  to  obtain  1%  statistical  accuracy  is  10,  000,  and  the  num¬ 
ber  for  1/2%  is  40,  000.  Thus,  it  would  be  easier  to  obtain  a  given 
accuracy  with  a  longer  counting  time,  for  example,  10  sec,  since  then 
the  cotinting  rate  need  not  be  too  high. 

c.  Radiation  Source  Types 

Due  to  the  atmospheric  densities  expected  (~.  15  to  15  psi)  it 
will  be  necessary  to  use  x  or  beta  sources  rather  than  alpha  sources. 
Useful  x-ray  sources  would  be  Cd-109  (22  keV),  Am-241  (60  keV) 
and  Co-57  (~100  keV).  Useful  beta  sources  would  be  Pm -147  (~225 
keV  max.  )  and  Sr-Y-90  (~2  MeV  max.  ).  All  of  these  have  half  lives 
greater  than  one  year  and  are  readily  obtainable. 

d.  Parameters 

While  estimates  of  the  expected  accuracy  have  been  given  above, 
it  will  be  necessary  to  consider  the  actual  accuracy  achievable  by  con¬ 
sidering  the  following  parameters  during  a  particular  program: 

•  source  types 

•  source  strength 

•  allowed  geometrical  configurations 

•  time  constants 

•  dimensions  of  sensed  volume 

The  tradeoff  between  these  various  parameters  will  be  determined  to 
define  the  optimum  system  and  accuracy  for  the  required  models. 
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SIGNAL- TO-NOISE  RATIO 
SMALL  SIGNAL  CONSIDERATIONS 


Usually  the  type  of  wave  energy  used  for  probing  a  medium  is 
chosen  to  make  the  medium  as  transparent  as  possible  (radar,  sonar) 
to  the  waves.  Sound  is  more  sensitive  to  flow  than  light  is,  because  the 
interaction  of  sound  and  flow  is  strong.  The  strong  interaction  of  sound 
and  flow  also  makes  ultrasonic  flow  measuring  devices  sensitive  to  flow 
noise  and  leads  to  larger  attenuation.  Therefore,  the  range  over  which 
a  successful  flow  measurement  can  be  made  is  determined  by  the  signal- 
to- noise  ratio  which  can  be  obtained.  The  maximum  signal  strength  one 
can  transmit  through  the  gas  in  acoustics  is  limited  in  principle  by  non¬ 
linear  distortion  of  the  signal.  In  fact,  the  peak  power  one  can  obtain 
from  a  transducer  which  fits  the  other  system  requirements  is  more 
often  the  limit.  The  signal  undergoes  considerable  distortion  and  is  dif¬ 
ficult  to  interpret  (in  terms  of  flow,  for  instance)  at  power  limits  below 
the  maximum  power  one  can  transmit.  The  problems  encountered  using 
the  maximum  transmitted  power  are  briefly  discussed  in  a  later  section. 
The  maximum  power  consistent  with  linear  propagation  will  be  considered 
in  the  small  signal  section. 

Small  signal  is  defined  here  as  the  maximum  power  which  can  be 
transmitted  across  the  required  gap  (£  )  without  shock  formation.  The 
pressure  level  Pmax  corresponding  to  this  power  level  is  (Beyer,  1965) 

2  c2 

*max  (7  +  1)  ul 

where  w  is  the  angular  frequency,  7  the  specific  heat  ratio,  c  the  sound 
speed,  and  p  the  density.  The  longer  the  path  and  the  higher  the  fre¬ 
quency,  the  lower  the  maximum  pressure  level  becomes. 

The  pressure  levels  for  small  signal  operation  at  various  fre¬ 
quencies  and  pathlengths  are  given  in  Table  2.  It  should  be  noted  that 
a  cw  ultrasonic  wave  with  a  pressure  amplitude  10"^-  atm  puts  10  watts/ 
cm^  into  the  gas.  While  the  power  level  is  acceptable  on  a  peak  pulse 
basis,  power  levels  of  0.  1  watts /cm^  (pmax  ~  10"^  atm)  are  more 
realistic  for  cw  operation.  Note  the  maximum  small  signal  pressure 
drops  directly  as  the  gas  pressure. 

The  first  part  of  the  present  section  is  devoted  to  calculation  of 
the  attenuation  of  small  acoustic  signals  in  the  gas  path.  This  includes 
losses  due  to  thermal  absorption,  turbulent  scattering,  refraction  due 
to  large  scale  flows,  and  beam  shape  factors.  The  turbulent  noise 
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pressure  will  then  be  discussed.  The  final  section  will  be  devoted  to 
putting  the  signal -to- noise  information  together  to  obtain  response 
time  versus  accuracy  for  applications  of  interest. 


Table  2.  Maximum  pressure  at  transmitter  (atm)  for 
small  signal  transmission  (1  atm,  300  K). 


Frequency,  kHz 

10 

50 

100 

500 

&  =  5  cm 

1.  31  x  10"1 

2.  6l  x  10  2 

1.  31  x  10"2 

2.  61  x  10"3 

i  =  25  cm  ' 

2.  62  x  10"2 

5.  24  x  10"2 

2.62  x  10"3 

5.  24  x  10-4 

i  =50  cm 

1.  31  x  10-2 

2.  6l  x  10"3 

1.  31  x  10"3 

2.  61  x  10"4 

Thermal  Absorption 


Thermal  absorption  includes  losses  of  sound  energy  to  heat  by 
viscosity,  thermal  conductivity  and  relaxation  (rotational  and  vibrational 
relaxation).  Thermal  conduction,  viscous  dissipation  and  rotational  re¬ 
laxation  dominate  the  true  sound  absorption  in  dry  air  under  the  temperature- 
pressure  conditions  contemplated  here.  The  Bound  absorption  due  to  thermal 
conductivity  and  viscosity  is  given  by 


a 

c 


?  2* 
2  TT  f 


P  c 


(4) 


Here  f  is  the  ultrasonic  frequency,  c  the  sound  speed,  p  the  density, 
7  the  specific  heat  ratio,  and  Pr  the  Prandtl  number.  The  total  loss 
including  rotation  relaxation  (Zmuda,  1951)  is 


a  -  a  (1  +  0.  067  Z  )  . 
c  K 


(5) 


The  rotational  collision  number  (Z^)  for  air  is  close  to  that  for  nitrogen 
and  oxygen  (Carnevale,  1967;  Brau,  1970;  Lordi,  1970).  The  value  of 
Zj^  at  300  K  in  air  is  Zj^  ~  3.  8  collisions  and  increases  slowly  with 
temperature. 


^Below  ~100  atm  p  c2 


=  7  p  is  a  useful  simplification. 
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The  attenuation  times  the  pathlength  is  plotted  versus  pathlength 
in  Fig.  8  for  air  at  300  K,  1  atm  for  frequencies  of  50,  100,  200  and  400 
kHz.  The  ratio  of  received  amplitude  to  transmitted  amplitude  is  also 

plotted  on  the  same  graph.  Note  that  there  is  very  little  (tH 

absorption  in  a  20  cm  path  at  100  kHz  and  below  due  to  thermal  effects. 
The  attenuation  at  any  temperature  and  pressure  can  be  scaled  from  the 
results  in  Fig.  8  using  the  relation 


a  L  (T,  p,  f) 


(6) 


where  f  is  the  operating  frequency  in  kHz,  T  the  temperature  in  °K,  and 
p  the  pressure  in  atm. 


Vibrational  relaxation  also  occurs  in  air.  However,  at  the  tem¬ 
peratures  and  pressures  of  interest  here,  in  dry  air,  the  relaxation  fre¬ 
quency  and  the  relaxation  strength  are  both  too  small  to  cause  any  effect. 
The  addition  of  water  vapor  increases  the  efficiency  of  collisional  transfer 
in  N£  and  C^. 

Relaxation-  Humidity 

Relaxation  of  air  is  very  complex  (Bass,  1972)  due  to  the  many 
molecules  which  make  it  up.  Water  vapor  has  catalytic  effects  on  the 
vibrational  relaxation  of  oxygen,  nitrogen  and  CO2.  The  addition  of  water 
vapor  can  shift  the  relaxation  time  of  20  C  air  from  less  than  100  Hz  to 
just  over  10^  Hz  at  100%  humidity  (Evans,  1974).  The  absorption  triples 
at  50  kHz  as  the  relative  humidity  of  the  air  goes  from  0-100%.  This 
multiplier  may  be  used  on  the  results  in  Fig.  8  to  include  the  effects  of 
relaxations  at  high  relative  humidities.  A  plot  of  total  absorption  for 
different  relative  humidities  taken  from  Evans,  1974,  is  given  in  Fig.  9. 
There  is  also  a  small  effect  on  the  sound  speed.  An  analytic  expression 
for  sound  speed  and  absorption  is  given  in  Bass,  1972.  As  pointed  out  by 
the  authors,  the  lack  of  good  relaxation  times  prevents  accurate  calcula¬ 
tions  of  sound  absorption  beyond  about  3  of  20  C.  These  problems  do 
not  affect  the  flow  measurement  since  sound  speed  is  measured,  not  cal¬ 
culated,  in  the  flowmeter  designs  contemplated  in  the  present  work. 


Turbulent  Scattering 

Scattering  of  sound  by  the  turbulent  flow  fields  can  cause  significant 
signal  attenuation  in  many  flows  of  aerodynamic  interest.  Theories  for 
scattering  of  sound  by  isotropic  turbulence  have  been  developed  and  tested 
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Fig.  8.  Thermal  sound  attenuation,  in  nepers, and  amplitude  ratios  in  air 
at  several  frequencies,  versus  path  length  L,.  Knowing  L  and  f,  one  uses 
this  graph  to  find  a  L  and  A/A  . 
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Fig.  9.  Total  absorption  of  sound  in  air  as  a  function  of  frequency 
(after  Evans,  1974). 
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for  oceanic  and  atmospheric  sounding.  However,  for  nonisotropic 
turbulence  such  as  is  the  case  for  boundary  layer  and  pipe  flow,  the 
best  one  can  do  is  appeal  to  the  notion  of  locally  isotropic  turbulence 
(Tartarsky,  1961,  pp.  66-70).  Even  when  one  is  willing  to  assume 
local  isotropy  the  required  correlation  functions  for  turbulent  tem¬ 
perature  and  flow  velocity  fluctuation  fields  are  not  usually  available 
for  the  flows  of  interest.  Below,  the  scattering  due  to  turbulence  is 
calculated  in  terms  of  quantities  that  are  readily  available  in  standard 
texts  on  turbulent  flow. 


The  attenuation  of  sound  from  turbulence  is  much  less  than  that 
due  to  viscous  and  thermal  losses  when  the  wavelength  is  much  longer 
than  the  scale  of  the  turbulence.  The  reason  for  the  small  scattering 
in  the  long  wavelength  region  can  be  seen  by  considering  the  size  of  the 
terms  in  the  absorption  expressions  one  finds  in  the  literature,  i.  e.  , 
(Chernov,  I960,  p.  55) 


~  2,43 

a  =  (jl  k  a 


(7) 


where  |i  is  the  mean  square  Mach  number  fluctuation.  The  mean 

2  -3 

square  Mach  number  fluctuation  fj.  is  typically  1  0  .  Since  the  wave 

number  turbulent  scale  product  (ka)  is  less  than  1,  the  absorption  never 
reaches  0.  01  neper/cm. 


The  scattering  for  ka  >>  1  can  easily  be  greater  than  the  thermal 
absorption.  However,  the  magnitude  of  the  calculated  absorption  can  vary 
by  factors  of  4,  depending  on  the  assumptions  made  about  the  structure 
of  the  turbulence.  The  calculation  below  is  based  on  formulas  given  by 
(Tartarsky,  1961 )  and  corrected  by  (Monin,  1962).  Assuming  that  the 
correlation  functions  are  exponential 


B 

rr 


for  velocity,  and 


B 


T 


(8) 


(9) 


for  thermal  fluctuations  the  absorption  due  to  turbulence  is  given  by 


a  = 


3  TT 
8  a 


A  u 


2  AT 


,2  2 

k  a 


(10) 
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2 

The  correlation  length  (ia),  the  mean  square  velocity  (Au  ),  and 
2 

temperature  (A  T  )  fluctuations  must  be  known  to  use  Eq.  (10).  The 
sound  speed  c,  the  mean  temperature  (T),  and  the  ultrasonic  wave 
number  (k)  are  all  known. 


The  correlation  length  was  taken  to  be  the  mixing  length 
(Schlichting,  I960,  p.  510)  in  the  absence  of  measured  correlation 
lengths.  The  mixing  length  for  a  pipe  varies  from  zero  at  the  wall 
to  about  14%  of  the  pipe  radius  at  the  centerline.  The  mixing  length 
for  a  turbulent  boundary  layer  (Rotta,  1972,  p.  350)  rises  from  zero 

to  10%  of  the  displacement  thickness  (6  )  at  g-j  ~  0.  2  and  remains 


y  /\  n 

constant  almost  to  6  .  The  turbulent  intensity  — —  is  taken  as 

constant  and  equal  to  0.  05  for  both  pipe  flows  and  boundary  layer  flows. 


Temperature  changes  in  turbulent  flow  will  usually  vary  about 
the  mean  temperature  by  the  same  factional  amount  as  the  flow  variations 


(11) 


Thus,  for  order  of  magnitude  calculations  such  as  are  contemplated  here, 
Eq.  (10)  becomes 


a 


2.  59  x  10 


(k  a)2  . 


(12) 


It  should  be  noted  that  where  aerodynamic  heating  or  cooling  are  involved 


can  easily  be  greater  than 


especially  at  low  Mach  numbers. 


For  instance,  if  a  turbulent  pipe  flow  with  a  3°C  temperature  difference 
between  the  walls  and  the  free  stream  flows  in  a  pipe  with  a  bend,  the  re¬ 
sulting  secondary  flows  convect  1%  temperature  differences  to  the  cooler 

A  -4  A  u'^  -6 

walls.  At  M  =  0. 1 , - —  is  ~  1 0  while — ~ —  ~1 0"  .  Thus,  at  low 

T  c 

Mach  numbers  Eq.  (12)  may  underestimate  the  absorption  due  to  turbulent 
scattering. 


36 


AEDCTR-74-77 


Absorption  due  to  turbulence  according  to  Eq.  (12)  has  been 
plotted  in  Fig.  10  as  a  function  of  pathlength  and  frequency  for  a 
M  =  0.  5  fully  developed  pifce  flow.  The  absorption  due  to  turbulence 
is  greater  than  the  thermal  absorption  (Fig.  8)  under  the  conditions 
calculated.  Numbers  can  be  scaled  from  Fig.  10  using  Eq.  (1Z). 
Note  that  the  absorption  due  to  turbulence  does  not  depend  on  pres¬ 
sure.  Thus,  at  sufficiently  low  pressures  thermal  absorption 
dominates  while  at  elevated  pressures  turbulent  losses  dominate. 

Harmonic  Generation 


One  might  expect  that  the  signal  attenuation  due  to  turbulence 
and  thermal  loss  can  be  added  to  give  a  total  loss  in  the  path.  How¬ 
ever,  we  defined  the  small  signal  power  level  so  that  the  path  is  one 
discontinuity  distance  (p.  Z2).  There  is  considerable  loss  at  these 
power  levels  to  harmonics  generated  by  the  nonlinear  mixing  of  the 
signal  with  itself  (Beyer,  1965).  The  effective  attenuation,  including 
losses  to  harmonic  generation,  are  as  much  as  four  times  the  absorp¬ 
tion  calculated  from  Figs.  8  and  10.  -A  plot  of  the  effective  attenuation 
divided  by  thermal  and  turbulent  absorption  {a  Q)  versus  normalized 
pathlength,  is  given  in  Fig.  11.  Notice  the  multiplier  on  attenuation 
is  largest  for  the  lowest  path  loss.  This  is  because  the  signal  strength 
and  therefore  the  nonlinear  losses  are  great  over  more  of  the  path  than 
is  the  case  for  a  large  a  Q.  Nonlinear  mixing  leads  to  diminishing  re¬ 
turns  on  increasing  driver  amplitude. 


Diffraction- Beam  Spread  (Near  Field,  Far  Field,  Steered  Beam) 


A  major  source  of  signal  loss  is  due  to  beam  spreading  and  mis¬ 
alignment  of  transmitter  and  receiver.  The  general  shape  of  the  field  of 
a  piston  radiator  is  shown  in  Fig.  12.  The  beam  is  highly  collimated  in 


the  near  field  zone  between  Z 


0  and  Z  = 


4 


Here,  d  is  the  diameter 


of  the  transducer,  X  the  wavelength  of  sound,  and  Z  the  distance  from  the 
transducer  face.  The  length  of  the  near  field  increases  as  the  frequency 
of  the  ultrasonic  wave  is  increased  (X  decreased),  and  as  the  diameter  of 
the  transducer  is  increased  the  near  field  extends  further  from  the  transmitter. 


The  near  field  transducer  with  its  wave  fronts  perpendicular  to  the 
flow  reaches  the  receiver  with  a  minimum  loss.  The  advantage  of  the  beam 
drift"1'  configuration  is  that  it  minimizes  alignment  losses.  A  near  field 

*The  transducers  positioned  so  that  the  receiver  is  just  downstream  enough 
to  catch  the  drift  of  the  beam  due  to  flow  is  called  the  beam  drift  configuration. 
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PATH  LENGTH,  CM 


Fig,  10.  Attenuation  in  air  due  to  turbulence,  in  nepers,  and  amplitude 
ratio,  versus  path  length. 
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Fig.  12.  Beam  characteristics  of  circular  piston.  Dimensions  shown  are  approximations 
suitable  for  estimating  purposes. 
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transducer  is  very  sensitive  to  tilt  in  the  wave  front.  To  maintain  a 
near  field  out  to  50  cm  at  100  kHz,  requires  a  transducer  with  a 
diameter  of  about  24  wavelengths  (8  cm).  A  misalignment  of  ~0.  13 
radians  causes  a  substantial  phase  shift  across  the  receiving  trans¬ 
ducer  face,  resulting  in  large  phase  jitter.  However,  the  average 
amplitude  is  preserved. 

The  other  extreme  is  in  the  wide  beam  transducer.  Referring 
again  to  Fig.  11,  a  wide  beam  is  obtained  by  decreasing  d.  Any  prac¬ 
tical  phase  measuring  system  for  ultrasonic  flow  measurement  requires 
upstream  and  downstream  signal  paths.  A  50  kHz  system  suitable  for 
measuring  flows  between  M  =  0.  2  and  M  =  0.  7  must  have  a  beam  width 
of  45  and  a  transducer  diameter  roughly  equal  to  one  wavelength  (X.  ) 
is  ~0.  66  cm.  Only  a  small  solid  angle  of  the  signal  is  subtended  by  the 
receiver  at  each  flow  Mach  number.  Thus,  the  pressure  loss  due  to 
spherical  spreading  in  a  25  cm  path  is  a  factor  of  400  or  -63  dB.  In 
fact,  the  loss  will  be  about  12  dB  greater  due  to  the  angle  between  the 
signal  wave  fronts  and  the  receive  transducer. 

Finally,  there  is  the  loss  of  signal  due  to  spreading  of  the  ray 
trajectories  for  a  realistic  flow  as  indicated  in  Fig.  12.  Note  that  three 
sound  sources  are  required  to  obtain  an  upstream  and  a  downstream  ray 
in  many  cases.  When  the  flow  velocity  gradients  are  sufficiently  sharp 
a  downstream  ray  will  return  to  the  wall  due  to  refraction.  This  phe¬ 
nomenon  called  shadow  zone  propagation  is  aggravated  by  flows  with  a 
hot  free  stream.  A  set  of  typical  rays  is  shown  in  Fig.  13  for  various 
launch  angles,  a Q.  In  addition,  Fig.  13  shows  downstream  rays  are 
spread  apart  leading  to  an  effective  attenuation  which  depends  on  the 
flow  and  temperature  gradients.  Ray  spreading  can  lead  to  an  additional 
20  dB  loss  for  downstream  rays. 

NOISE 


Source  of  Noise 


Acoustic  noise  generated  by  the  flow  itself  may  be  the  limiting  factor 
on  the  accuracy  or  feasibility  of  many  acoustic  measurements.  The  pres¬ 
sure  fluctuations  on  the  face  of  the  transducer  pn  and  the  signal  ps  enter 
the  detection  system  together.  Thus,  there  is  nothing  that  can  be  done 
to  improve  the  signal-to-noise  ratio  by  electrical  or  acoustic  nonresonance 
impedance  matching  on  the  receive  side.  However,  increasing  the  trans¬ 
ducer  size,  increasing  operating  frequency  and  narrow  band  filtering, 
discriminate  against  the  turbulent  noise.  The  most  noise  reduction  one 
can  expect  from  transducer  size  is  a  factor  of  10^.  Any  attempt  to  broaden 
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Fig.  13.  Schematic  of  pulse  trajectory  and  wave  fronts  in  a  real  flow. 
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or  steer  the  beam  forces  the  use  of  smaller  transducers,  thus  increas¬ 
ing  the  noise.  A  narrow  band  system  can  reduce  the  noise  at  the  expense 
of  response  time.  As  discussed  below,  there  are  practical  limits  to  how 
much  the  bandwidth  reductions  can  improve  the  signal-to-noise  ratio. 
However,  a  factor  of  10^  relative  to  the  total  wide  band  noise  pressure 
fluctuation  can  probably  be  obtained  with  cw  acoustic  transmission.  Under 
turbulent  flow  conditions  in  air  at  Ms  =  0.  5,  the  turbulent  wall  pressure 
noise  ranges  between  1.  78  x  10^  dyne/cm^  to  ~3.  46  x  10^  dyne/cm^.  Using 
optimum  filtering  and  transducer  size  an  equivalent  acoustic  pressure  of 
1.  78  x  10-5-8.  46  x  10"1  dyne/cm^  enters  the  transducer  with  the  signal.* 

The  noise  generated  by  turbulence  depends  to  a  large  extent  on  the 
type  of  flow  involved  (Sidden,  1973).  The  fraction  of  the  dynamic  head 
(p  U^/2)  converted  into  noise  pressure  fluctuations  in  different  types  of 
flow  varies  over  a  three  order  of  magnitude  range  as  shown  in  Table  3. 

True  sound  radiated  from  distant  turbulent  sources  makes  only  a  small 
contribution  to  the  total  flow  noise.  The  true  sound  level  is  representative 
of  the  noise  level  one  should  expect  from  a  transducer  mounted  on  a  bound¬ 
ary  layer  stripper  or  on  a  model  in  a  low  turbulence  tunnel. 


Table  3.  Maximum  percent  of  dynamic  head  which  appears  as 
wall  pressure  fluctuations  in  various  flow  situations. 


True 

Acoustic 

Radiation 

P  s  eudo-Sound 
Under  Turbulent 
Boundary  Layer 

Separated  Flow  Near 
Reattachment  Zone 

Lift  and  Drag 
Fluctuations  Due 

To  Freestream 
Turbulence  Over 
Thin  Boundary  Layer 

0.  01 

0.  6 

6 

10 

The  noise  under  a  thick  (~0.  5")  turbulent  boundary  layer  or  in 
fully  developed  pipe  flow  is  due  to  pressure  fluctuations  on  the  wall  itself 
(Haddle,  1969  and  I960).  The  rms  pressure  fluctuations  are  typically  0.  6% 
of  the  dynamic  pressure.  Separation  (Frick,  1968,  1970,  and  Siddon,  1973) 
of  the  boundary  layer  leads  to  a  factor  of  ten  increase  in  the  wall  pressure 
near  reattachment. 


*This  should  be  compared  with  the  signal  levels  of  10“2-10“9  dyne/cm^ 
which  are  likely  to  be  received  in  practice  as  indicated  in  the  previous 
s  ection. 
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The  fluctuating  lift  and  drag  forces  due  to  a  turbulent  free- 
stream  can  cause  very  large  pressure  fluctuations  up  to  20%  of  the 
dynamic  pressure.  Lift  and  drag  pressures  are  important  sources 
of  noise  on  models  (Siddon,  1973)  near  leading  edges  where  the 
boundary  layers  are  thin.  Transducers  which  are  immersed  into 
the  flowing  stings  or  acoustic  wave  guides  develop  large  drag  noise 
signals.  In  general,  the  acoustic  wave  guides  and  transducers  must 
be  carefully  isolated  from  those  portions  of  wall  which  protrude  into 
the  flow. 

The  order  of  magnitude  of  the  acoustic  noise  and  some  general 
constraints  on  the  acoustic  flowmeter  system  can  be  determined  from 
the  general  information  given  above.  Probes  which  protrude  or  are 
recessed  from  the  duct  walls  will  result  in  a  factor  of  30  increase  in 
turbulent  noise  over  that  due  to  the  turbulent  boundary  layer.  Thus, 
flush-mounted  transducers  are  advisable.  Similarly,  the  noise  levels 
on  models  in  a  turbulent  freestream  may  be  larger  than  the  noise  on 
the  duct  walls  by  a  similar  factor.  Once  the  expected  wide  band  noise 
levels  are  calculated  for  the  type  of  flow  involved,  the  dynamic  head, 
and  the  factor  in  Table  3,  the  bandwidth  and  transducer  factors  can  be 
used  to  obtain  a  first  estimate  of  the  final  signal -to -nois  e  level.  Typ¬ 
ically,  for  a  1  second  response  time  a  10^  reduction  in  flow  noise  can 
be  expected  with  cw  operation.  This  improvement  is  reduced  by  the 
duty  cycle  for  pulse  operation.  See  also:  Camp  et  al.  ,  1970. 

More  exact  noise  information  shows  that  the  flow  noise  drops 
off  as  some  high  (~3)  power  of  frequency.  In  addition,  a  more  exact 
calculation  of  spatial  filtering  due  to  transducer  size  allows  the  noise 
penalty  for  small  transducers  to  be  calculated  (Corcos,  1967;  Racine, 
1971;  Petrovskii,  1973;  Smolyakov,  1970).  Many  methods  for  obtaining 
rays  over  a  wide  solid  angle  requires  the  use  of  small  transducers. 
Therefore,  the  remainder  of  this  section  is  devoted  to  more  detailed 
calculation  of  turbulent  noise. 

The  most  intense  source  of  noise  is  due  to  the  boundary  layer 
flow  near  the  wall.  The  noise  from  turbulent  boundary  layers  has  been 
exhaustively  studied  by  aerodynamicists  and  sonar  specialists.  How¬ 
ever,  emphasis  has  been  on  the  low  frequency  end  of  the  noise  spectra. 
Fortunately,  in  the  megahertz  range  the  noise  spectrum  is  dominated 
by  the  balance  of  energy  transfer  to  the  higher  frequencies  via  nonlinear 
mixing  (Ellison,  1962;  Kuznetzov,  1970).  The  turbulent  energy  is 

jL 

With  the  exception  of  low  turbulence  freestreams  where  boundary  layer 
strippers  are  advisable  to  avoid  the  boundary  layer. 
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NORMALIZED  RMS  NOISE  PRESSURE  PER  HERTZ 


Fig.  14.  Normalized  rms  pressure  per  hertz  versus  normalized 
frequency. 
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introduced  to  the  flow  at  low  frequencies  near  the  wall.  Therefore,  the 
high  frequency  noise  power  universally  falls  off  as  f^. 


A  convenient  analytic  form  for  normalized  mean  square  pressure/ 
unit  frequency  (power  spectral  density  at  frequency  u)  vs  normalized  fre¬ 
quency  was  taken  from  Strawderman,  1968.  The  normalization  used  in 
Fig.  14  correlates  boundary  layer  noise  from  many  extremely  different 
types  of  flows  in  liquids  and  gases.  This  graph  is  based  on  the  analytical 
expression  for  mean  square  pressure  given  by 


A 


0.  75  x  10 


2 

a 


(13) 


(14) 


and  p£  is  the  fluid  density,  6  *  the  boundary  layer  displacement  thickness, 
or  the  pipe  radius  for  fully  developed  pipe  flow,  and  a  a  dimensionless 
quantity: 


1.  0  (for  air) 

(15) 

3.  0  (for  water) 

(16) 

and  U0  is  the  freestream  flow  velocity.  The  convection  frequency  is 

U 


o 


The  pressure  fluctuations  on  an  arbitrarily  small  transducer  are 
constant  out  to  a  frequency  w  =  1.  256  wQ.  The  noise  pressure  then  falls 
off  as  Thus,  it  is  advantageous  to  operate  with  a  carrier  frequency 

which  is  as  high  as  the  attenuation  will  allow.  The  normalized  rms  pres¬ 
sure  amplitude  versus  frequency  is  shown  in  Fig.  14.  This  graph  is  uni¬ 
versal.  The  noise  pressure  on  the  crystal  face  can  be  calculated  from  the 
displacement  thickness  (or  pipe  radius,  R),  centerline  flow  velocity  UQ  and 
dynamic  head  (p  UQ/2)  for  any  frequency. 

The  noise  pressure  is  not  coherent  across  the  face  of  the  trans¬ 
ducer.  Thus,  the  noise  which  enters  the  transducer  is  reduced  by  phase 
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cancellation  over  the  transducer  face  (Corcos,  1967,  I960).  The  can¬ 
cellation  effect  is  improved  somewhat  when  the  transducer  diameter  is 
larger  than  the  displacement  thickness  (Smolyakov,  1970).  Cancellation 
is  less  complete  when  transducer  sensitivity  is  not  constant  over  the 
transducer  face  (Petrovskii,  1973).  The  signal  is  more  coherent  than 
the  noise  over  the  transducer  face  and  does  not  suffer  as  much  attenua¬ 
tion  due  to  spatial  filtering. 


A  spatial  filtering  factor  was  taken  from  Recine,  1971,  and  is 

1/2 


given  in  Fig.  15.  The  existing  measurements  end  at 


0)  A 


U 


=  15.  The 


u)  A 


1/2 


data  were  extrapolated  to  — — -  =  200  for  Fig.  15.  The  correlation 


of  pressure  fluctuations  across  the  transducer  face  and  the  spatial  varia¬ 
tion  of  transducer  sensitivity  must  be  known  to  determine  the  effect  of 
transducer  size  exactly.  However,  the  product  of  the  noise  pressure 
from  Fig.  14  times  the  spatial  factor  from  Fig.  15  gives  an  order  of 
magnitude  of  the  noise  to  be  expected  from  typical  transducers  subject 
to  turbulent  boundary  layers.  The  percentage  factors  in  Table  3  should 
be  applied  to  the  turbulent  boundary  layer  predictions  to  arrive  at  the 
frequency  dependent  flow  noise  in  Hie  other  types  of  flow  enumerated  in 
Table  3. 


NOISE  IN  DUCTS  VERSUS  MACH  NUMBER 


The  noise  in  a  25  cm  diameter  duct  with  a  fully  developed  pipe 
flow  is  of  specific  interest  for  the  tests  in  the  PWT  facilities.  The  rms 
noise  pressure  appearing  at  the  electrical  terminals  of  a  2.  5  cm  trans¬ 
ducer  in  the  25  cm  duct  is  shown  in  Fig.  16.  This  transducer  is  the  size 
that  would  be  used  in  a  transmission  (upstream-downstream)  measure¬ 
ment.  The  frequency  at  which  the  noise  falls  off  varies  from  654  Hz  at 
Ms  =  0.  2  to  2800  Hz  at  Ms  =  0.  8.  The  rms  noise  pressure  decreases  as 
the  cube  of  the  flow  velocity  above  thiB  frequency.  Suppose  we  wish  to 
receive  a  signal  after  processing  which  is  10  to  100  times  the  pressure 
level  generated  by  turbulence  at  the  highest  Mach  number  of  interest.  This 
means  that  a  system  designed  to  operate  up  to  Mach  0.  8  with  a  bandwidth  of 
10  Hz  at  ~100  kHz  must  receive  a  signal  level  of  10~^-10~*  dynes/cm^.  The 
conditions  under  which  this  signal  level  can  be  achieved  are  deduced  from 
the  signal  level  discussion  on  p.  41. 


The  noise  pressure  in  other  duct  diameters  can  be  scaled  from 
the  values  in  Fig.  16  using  Eq.  (17).  Thus,  the  rms  noise  pressure  p^ 
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Fig.  15.  Spatial  filtering  factor  for  finite  transducers  of  diameter  (d). 
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1/2 

Fig.  16.  Turbulent  RMS  noise  pressure/(Hz)  ,  calculated  for 
empirical  data  for  a  fully  developed  pipe  flow  in  a  25  cm  diameter 
duct  with  a  2.  54  cm  diameter  transducer. 
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in  a  fully  developed  pipe  flow  in  a  duct  of  diameter  D  is 


Here  pQ  is  the  noise  pressure  from  Fig.  l6. 

Similarly,  to  obtain  the  noise  in  ducts  where  the  flow  is  not 
fully  developed  the  displacement  thickness  2  6*is  put  into  Eq.  (17)  in 
place  of  the  duct  diameter  D.  Noise  pressures  for  the  small  diameter 
~5  cm  ducts  in  the  "ultrasonic  venturi"  would  be  determined  by  the 
displacement  thickness  of  the  ducts.  It  is  important  in  the  "ultrasonic 
venturi"  application  that  the  ultrasonic  transducers  be  located  at  least 
5  cm  downstream  from  the  leading  edge  of  the  duct  to  avoid  lift  and 
drag  noise.  In  addition,  the  operating  frequency  can  proportionately 
be  increased  in  smaller  ducts.  For  instance,  if  50  kHz  is  appropriate 
for  a  25  cm  duct,  250  kHz  can  be  used  in  a  5  cm  duct.  Looking  at  the 
noise  fall-off  with  frequency  in  Fig.  16,  a  factor  of  25  noise  reduction 
is  obtained  by  increasing  the  frequency  of  the  ultrasonic  wave  by  this 
factor  of  5.  This  makes  the  ultrasonic  venturi  a  very  attractive 
prospect. 

SMALL  SIGNAL- TO-NOISE  RATIO 


The  signal  and  noise  can  now  be  put  together  to  see  if  sufficient 
signal-to-noise  ratio  can  be  achieved  to  perform  the  desired  measure¬ 
ment.  The  case  of  fully  developed  pipe  flow  in  a  25  cm  duct  with  0.  6 
cm  spherical  radiators  operating  at  several  frequencies  were  calculated. 

The  signal-to-noise  ratios  are  shown  as  a  function  of  Mach  number  in 
Fig.  17  for  a  beam  drift  trajectory.  The  signal-to-noise  ratio  is 
optimum  at  100  kHz.  However,  sufficient  signal-to-noise  level  is  ob¬ 
tained  at  50  kHz  to  do  the  measurement.  It  is  more  convenient  to  gen¬ 
erate  a  50  kHz  sound  wave  at  the  small  signal  maximum  pressure  than 
it  is  at  100  kHz. 

The  response  time  required  to  make  a  measurement  to  a  given 
precision  in  the  presence  of  uncorrelated  noise  decreases  with  the  signal- 
to-noise  ratio.  The  exact  relation  between  response  time,  precision,  and 
signal-to-noise  ratio  depends  on  details  on  the  quantity  being  measured 
and  the  measurement  system.  However,  the  typical  precision  one  should 
expect  is  one  minus  the  amplitude  signal-to-noise  ratio.  The  response  time 
required  to  make  the  measurement  to  this  precision  is  the  one  over  the  band¬ 
width  of  the  receiver.  Figure  17  is  calculated  for  a  constant  bandwidth  of 
10  Hz.  Thus,  the  greatest  precision  we  should  expect  from  a  0.  1  second 
response  time  system  at  the  frequencies  and  Mach  numbers  given  in  Fig.  17 
is  %  precision  =  /l  —  _  \  where  —  is  the  signal-to-noise  from  Fig.  17. 
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MACH  NUMBER 

Fig.  17.  Signal-to-noise  ratio  in  a  25  cm  duct  with  a  2.  5  cm  diameter 
transducer. 
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OTHER  ANALYTIC  OR  DESIGN  CONSIDERATIONS 
GAS  PRESSURE  EFFECTS  ON  SIGNAL  LOSSES 


The  gas  pressure  p  has  two  principal  effects  on  classical  sound 
propagation: 

(1)  The  sound  energy  transmission  coefficient  increases  as  p 
increases,  as  the  probe/ gas  mismatch  is  lessened; 

(2)  The  classical  absorption  coefficient  decreases  as  p 
increases. 

Other  effects  of  p,  relative  to  maximum  signal  intensities  that 
to  not  yet  introduce  nonlinear  or  "finite  amplitude"  effects,  and  relative 
to  noise,  are  not  considered  here.  The  effect  of  temperature  T  is  to 
increase  slightly,  the  losses  Lz  and  La  ,  as  shown  in  "Table  II",  p.  719, 
of  the  cited  paper  on  ultrasonic  temperature  measurement  (Lynnworth 
and  Carnevale,  1972). 

To  graphically  show  the  two  principal  effects  listed  above,  in 
Fig.  18  we  have  plotted  lines  of  constant  loss  due  to  classical  absorption, 
La,  and  a  line  of  loss  2LZ  due  to  impedance  mismatches  at  transmitter 
and  receiver  sides.  Lz  is  independent  of  pathlength  and  frequency.  Lff 
increases  in  proportion  to  pathlength  and  in  proportion  to  frequency 
squared.  The  graph  shows  that  for  most  conditions  of  interest  in  Re¬ 
quirement  2,  Lff  is  very  large  at  1  MHz,  but  reasonably  smaller,  by  a 
factor  of  100,  at  0.  1  MHz.  (Losses  due  to  turbulence  have  not  been 
calculated.  ) 

The  additional  diffraction  (beam  spread)  loss  at  0.  1  and  1  MHz 
is  readily  calculated,  for  an  assumed  transducer  diameter  2  a  =  3  cm,  a 
sound  speed  in  the  gas  of  c  =  300  m/s,  and  a  path  z  =  30  cm  (»  1  ft).  The 
dimensionless  Seki  parameter  S  =  z\/a^  =  z  c/a^f  =  4  at  0.  1  MHz,  and 
0.  4  at  1  MHz.  The  corresponding  diffraction  losses  are  3.  85  dB  and  1.  2 
dB.  These  values  are  negligibly  small  relative  to  Lz,  suggesting  that  low 
frequencies  (~0.  1  MHz)  can  be  used,  without  significant  penalty  due  to 
beam  spread  losses. 

COMPARISON  OF  ULTRASONIC  FLOWMETERS  FOR 
GASES  VERSUS  LIQUIDS 


On  page  15  we  identified  six  (6)  basic  types  of  ultrasonic  flowmeters. 
All  but  the  correlation  type  have  been  manufactured  as  commercial  instru¬ 
ments  for  application  in  gases  and  liquids.  Based  on  data  in  the  literature. 
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Fig.  18.  Calculated  losses  due  to  acoustic  impedance  mismatch,  2  Lz 
and  classical  absorption,  h  ,  at  1  MHz  (and  at  0.  1  MHz)  in  nitrogen 
at  300K,  as  functions  of  pathlength  and  gas  pressure.  Calculations 
assume  a  solid  buffer  impedance  of  1.  3  g/cm^-ps,  and  neglect  beam 
spread. 
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and  manufacturer'  product  information  sheets,  and  our  analysis,  we 
can  compare  some  general  characteristics  of  ultrasonic  flowmeters. 
The  purpose  of  Table  4,  p.  55,  is  to  afford  some  perspective  on 
approaches  and  on  the' practicality  of  these  devices,  and  provide  an 
estimate  of  accuracy. 

ELECTROSTATIC  TRANSDUCERS 


Several  transducers  using  flexure  to  match  intense  acoustic 
signals  to  air  in  the  tens  of  kHz  range  (Barone,  1971;  Massa,  1965), 
have  been  used.  Renewed  interest  in  electrostatic  transducers  for 
this  purpose  has  come  into  use  recently  (Kuhl,  1954;  Legros  and 
Lewiner,  1973).* 

The  electrostatic  transducer  has  additional  advantages  for  ease 
of  installation  and  operation.  An  electrostatic  transducer  schematic  is 
shown  in  Fig.  19.  The  transducer  is  made  of  inert  flexible  films  which 
can  easily  be  installed  flush  on  curved  surfaces.  In  addition,  the  stiffness 
of  the  transducer  is  the  air  cushion  behind  the  plastic  films.  Thus,  the 
impedance  match  is  better  than  is  obtained  with  other  transducers. 

The  initial  design  for  the  flow  transducer  is  a  simple  broad  beam 
radiator.  However,  a  phased  array  offers  the  possibility  of  recovering 
most  of  the  80  dB  signal  loss  due  to  spherical  spreading.  Two  such  electro¬ 
static  arrays  are  shown  in  Fig.  19.  The  array  in  Fig.  19  shows  slots  in 
the  transducer  backup  plate  which  vary  across  the  transducer  face.  The 
first  groove  is  cut  to  k  /4  at  50  kHz,  k  /4  ~  0.  17  cm.  The  other  grooves 
vary  in  even  steps  down  to  0.  03  cm  deep.  The  k  /4  groove  will  cause  res¬ 
onant  phase  shift  (90°)  in  the  motion  of  the  portion  of  the  membrane  just 
above  it  (Kuhl,  1954).  The  phase  shift  of  each  succeeding  portion  of  mem¬ 
brane  will  drop  off  towards  45  .  Thus,  the  transducer  when  driven  at  the 
frequency  for  which  the  first  groove  is  k  /4  will  radiate  a  pencil  beam  at 
45°.  For  lower  frequencies  the  array  will  sweep  its  beam  towards  the 
normal.  A  focused  transducer  realization  is  shown  in  Fig.  19  which  op¬ 
erates  on  the  small  principle. 

TRANSDUCER  MOUNTING  -  DESIGN  CONFIGURATIONS 


The  schematics  in  Fig.  20  illustrate  nonrotated,  small  signal  plane 
wave  methods  for  interrogating  on  a  gas  path  along  the  three  directions 
which,  relative  to  the  duct  axis,  are:  parallel,  oblique  and  orthogonal. 

In  one  of  the  oblique  beam  drift  cases,  the  triangular  transducer  mosaic 
is  shown  connected  differentially  to  yield  an  output  proportional  to  the  Mach 
No.  v/c.  This  mosaic  is  an  "analog"  or  continuous  version  of  the  "digital" 
or  discrete  transducer  arrangement  indicated  in  the  British  NPL  work. 

*Electret  transducers  subject  to  large  electrode  displacements  are  analyzed 
by  G.  M.  Sessler,  J.  Acoust.  Soc.  Amer.  55,  2,  345-349  (Feb.  1974). 
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Table  4.  Comparison  of  typical  parameters,  relative  to 
ultrasonic  flowmetry. 


Parameter 

Gas 

Liquid 

Usual  approach 

1.  3* 

1,  2 

Common  path  up 

Along  axis,  M  <  1, 

Yes,  since  M  <<  1 

and  down 

or  oblique,  M  <<  1 

Sound  speed  c,  m/s 

~300 

~1000 

3 

Density,  g/cm 

10“4  at  NTP 

1 

Transmission  loss,  dB 

>  100 

~20 

Mach  No.  M 

~1  0  to  ~1 0* 

~io"3 

Gradients 

Large 

Small 

Interaction  of  sound  beam 
and  fluid 

Large 

Small 

Focusing  by  circular  pipe 

On- axis 

Off-axis 

Fluctuations  in  Ac/c  vs 
those  in  Av/v 

Small 

Large 

Ratio  of  c 

solid  fluid 

20:1 

5:1 

Response  time  required, 
seconds 

-1  4 

10  to  10 

10"Z  to  102 

Rangeability  required 

0  to  hypersonic 

100:1 

Accuracy,  %  of  range, 
standard  models 

0.  5% 

0.  5% 

*(1)  Upstream  minus  downstream;  (2)  Doppler  shift  off  scatterers; 
(3)  Beam  drift;  (4)  Correlation;  (5)  Noise;  (6)  Vortex  shedder. 
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Fig.  19.  Electrostatic  realization  o£  arrays. 
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SCHEMATIC 


FLOW  INLET  OUTLET 
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Relative  to  Duct  Axis 


Offset) 

orthogonal 
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parallel 


Parallel 
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Parallel 
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v/e  »0 — ' 


Drift  and 
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Oblique 
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Fig.  20.  Nonrotated  small  signal  plane  wave  methods. 
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PIEZOELECTRIC  CRYSTALS  AND  CERAMICS 


Since  several  books  have  been  devoted  to  this  topic  (Mason,  1950; 
Mattiat,  1970)  we  will  only  summarize  a  few  important  considerations. 
The  choice  of  piezoelectric  material  depends  on  what  characteristics  are 
most  important.  For  example,  the  short  list  below  could  be  used  to 
select  a  piezoelectric  material  to  emphasize  activity,  linearity,  high 
temperature  or  broadband  characteristics: 


MATERIAL  SELECTED 


TO  EMPHASIZE 


PZT 

Quartz 

Lithium  Niobate 
Lead  Metaniobate 


Activity 
Linearity 
High  Temperature 
Broadband 


Practical  applications  will  usually  involve  optimizing  the  response 
to  a  number  of  different  requirements.  These  may  typically  include,  in 
addition  to  the  characteristics  listed,  electrical  and  acoustical  impedance, 
mode  purity,  ^  ^eiec*  temperature  coefficients,  sensitivity,  size, 

aging,  availability,  compatibility  with  environmental  factors  other  than 
temperature,  such  as  humidity,  vibration,  pressure,  and,  occasionally, 
nuclear  radiation. 


As  far  as  we  are  aware,  all  present  commercial  ultrasonic  flow¬ 
meters  employ  transducers  of  the  piezoelectric  type.  To  illustrate  their 
use,  we  obtained  several  oscillograms  of  transmission  across  stationary 
(no-flow)  air,  up  to  1500  psig,  room  temperature  and  at  slightly  elevated 
temperature,  and  a  pathlength  of  69  cm  (27  in.  ).  These  transducers  were 

designed  to  have  low  Q  ,  and  were  operated  under  spike  and  also  rf 

mech 

excitation  at  0.  5  MHz.  Figure  28,  p.  77,  showing  one  bench  test  result 
with  these  piezoelectric  transducers  is  contained  in  the  section  on  Experi¬ 
mental  Results.  We  have  also  made  extensive  use  of  piezoelectric  trans¬ 
ducers  in  ultrasonic  measurements  in  shock  tubes  and  in  other  high  tem¬ 
perature  sources. 

As  the  operating  frequency  is  reduced  below  ~1  MHz,  however, 
say,  to  the  decades  10-100  kHz  or  100-1000  kHz,  magnetostrictive  and 
electromagnetic  transduction  become  increasingly  effective.  This  is 
because  the  field  of  the  exciting  wave  can  penetrate  deeper  and  more 
efficiently  into  the  transducer  material  (skin  depth  =  I/'s/  tt  f  jjl  or )  as 
frequency  f  decreases.  On  the  other  hand,  the  capacitive  reactance 
Xc  =  1/2  ir  f  C  of  a  piezoelectric  transducer  increases  as  f  decreases. 
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which  can  aggravate  matching  and  cable  loading  problems.  Therefore, 
we  briefly  examine  magnetostrictive  and  electromagnetic  transducers, 
in  the  following  paragraphs. 

MAGNETOSTRICTIVE  METALS 


A  variety  of  metals  are  commonly  used  as  magnetostrictive  trans¬ 
ducers  (Kikuchi,  pp.  44,  45,  in:  Mattiat,  1970).  We  have  mainly  used 
remendur,  an  alloy  principally  of  48%  Co,  48%  Fe,  4%  Va.  This  alloy 
has  relatively  high  coupling  and  high  curie  point,  nearly  1000  C.  We 
have  conducted  bench  tests  on  several  disc  resonators,  excited  axially, 
using  a  rod  of  remendur,  1. 4  mm  (0.  056  in.  )  diameter  by  >  15  mm  long, 
to  propagate  across  paths  several  hundred  mm  long,  at  f  <,  100  kHz.  The 
discs  may  be  excited  at  their  center  by  one  remendur  rod,  or  at  a  number 
of  points  essentially  along  or  near  a  nodal  path  such  as  a  mid- radius  circle. 
The  beam  patterns  thereby  generated  have  been  analyzed  to  some  degree, 
for  the  discs  themselves  (Barone  and  Juarez,  1972).  The  patterns  may  be 
controlled  to  peak  at  preferred  off-axis  directions  such  as  45°,  for  ex¬ 
ample.  Thus,  they  may  be  useful  in  beam  steering.  They  may  also  be 
installed  in  a  manner  that  hardly  perturbs  a  given  duct.  Installation  can 
even  take  advantage  of  an  all-metal  seal,  since  the  energizing  field  is  of 
low  enough  frequency  to  penetrate  a  tubular  sheath  such  as  stainless  steel, 
of  wall  thickness  ~0.  2-0.  5  mm  thick.  This  construction  is  illustrated 
below,  and  has  been  found  useful,  at  least  as  far  as  sealing  and  trans¬ 
duction  is  concerned,  in  nuclear  thermometry  applications  (Lynnworth,  1971). 


Fig,  21.  Magnetostrictive  transducer  used  with  all-metal 
seal,  with  nonmetal  gasket,  and  in  direct  coupling  to  duct. 
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ELECTROMAGNETIC  TRANSDUCTION 


Thorough  accounts  of  this  type  of  transduction  are  now  avail¬ 
able  in  the  literature  (Gordon,  1969;  Thompson,  1973;  Dobbs,  pp. 

127-191,  in  Mason  and  Thurston,  1973).  Briefly,  in  electrically 
conductive  materials,  through  the  Lorentz  force,  ultrasonic  waves 
can  be  generated  or  detected  through  the  interaction  of  the  ultrasonic- 
frequency  oscillation  and  a  dc  biasing  magnetic  field. 

The  effect  is  generally  much  weaker  than  for  piezoelectric  or 
magnetostrictive  transduction,  and  so  it  has  seen  little  industrial  use 
so  far.  However,  through  improved  signal  processing,  whereby  ade¬ 
quate  S/N  can  be  achieved,  this  approach  should  be  considered,  be¬ 
cause  of  its  several  interesting  features.  These  include:  the  convenient 
launching  of  particular  modes  or  polarization,  operation  unlimited  by 
temperature  per  se  (e.  g.  ,  no  curie  point  limit),  noncontact,  and  flexible 
geometrical  constraints.  Thompson  has  shown  that  the  method  is  well 
suited  for  certain  measurements  of  ultrasonic  phase  and  group  velocity, 
and  possibly  for  detecting  defects  in  large  steel  pipes.  These  applications 
represent  both  laboratory  and,  ultimately,  rugged  field  tests.  If  such 
field  applications  prove  practical,  then  the  extension  of  this  approach  to 
measurements  in  ducts  or  other  situations  may  be  a  worthwhile  pursuit. 

PRECISION,  BIAS,  STABILITY  AND  RESPONSE  TIME 


The  precision  of  the  electronic  phase  measurement,  if  done  in  a 
digital  manner,  is  +  1  count.  This  may  be  +  0.  1  out  of  a  range  of  Hr  180  . 
Even  finer  resolution  is  obtainable  on  commercial  phase  meters,  e.  g.  , 

+  0.  01  .  However,  other  factors  will  surely  limit  precision  so  that  elec¬ 
tronic  efforts  to  measure  phase  to  better  than  ~0.  1%  are  considered  un¬ 
warranted,  for  the  anticipated  applications. 

The  bias  in  the  optimum  design  will  probably  be  dominated  by  path 
or  propagation  differences  in  the  upstream  versus  downstream  directions. 

To  the  extent  that  the  paths  are  noncongruent,  and  to  the  extent  that  the 
propagation  "constant"  is  not  actually  constant,  even  over  the  brief  inter¬ 
rogation  interval  of  ~1  ms,  there  will  not  be  a  complete  cancellation  of 
the  effects  of  c.  Thus,  the  v  result  will  be  biased  by  the  uncanceled 
influence  of  c. 

Depending  on  the  actual  flow  and  temperature  variations  encoun¬ 
tered  in  the  PWT  facilities,  appropriate  selection  of  the  time  constant 
may  average  out  much  of  the  short-term  c  bias.  As  mentioned  elsewhere, 
the  time  constant  may  also  determine,  to  some  extent,  the  influence  of 
modulation  noise  (phase-jitter)  that  the  measuring  system  does  not  suppress. 
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Since  the  stated  response  time  in  this  program  is  given  as  <  1  second, 
for  Requirement  2,  and  this  is  essentially  the  same  as  required  in  the 
Eustis  fuel  flow  program  (0.  5  second)  it  is  expected  that  the  same  com¬ 
putational  times  will  suffice  in  the  present  gas  flow  metering  project. 

INSTALLATION,  OPERATION  AND  MAINTAINABILITY 

The  optimized  electronics  proposed  for  the  PWT  facilities  tests 
would  be  maintained  the  same  as  other  electronic  test  equipment  of 
similar  complexity. 

The  electrostatic  transducers  proposed  for  the  PWT  facilities 
tests  have  survived  repeated  starts  and  runs  (in  some  cases  over  ~100 
runs)  in  the  wind  tunnel  duct  used  in  our  bench  tests.  Other  electro¬ 
static  transducers  have  passed  similar  tests  (Schmidt  and  Tilmann,  1969). 
These  type  transducers  can  be  repaired  in  some  cases.  However,  the 
simplest  procedure  will  probably  be  their  direct  replacement,  should 
failure  occur.  The  two  main  failure  modes  appear  to  be:  (a)  tearing 
due  to  shear  stresses  near  the  wall;  (b)  electrical  breakdown  of  the 
dielectric. 

As  mentioned  elsewhere,  in  the  transducer  section,  piezoelectric 
materials  have  been  used  over  wide  pressure  ranges,  and  over  some 
temperature  ranges,  with  negligible  maintenance. 

Electromagnetic  transducers  would  probably  exhibit  mechanical 
problems  similar  to  the  electrostatics,  and  would  be  sensitive  to  coil 
burnout,  but  relatively  insensitive  to  dielectric  breakdown. 

Magnetostrictive  transducers  appear  to  be  the  most  rugged,  and 
would  therefore  be  the  easiest  to  maintain,  in  the  sense  that  they  would 
generally  require  the  least  maintenance.  Such  transducers,  for  example, 
have  been  used  for  thousands  of  hours  in  ultrasonic  thermometers  for 
nuclear  reactors. 

CALIBRATION 


Calibration,  historically,  has  often  been  a  difficult  area  to  re¬ 
solve,  for  ultrasonic  flowmeters.  In  some  cases  this  is  due  to  the 
differences  between  ultrasonic  readings  and  reference  readings  being 
attributable  to  errors  in  the  reference. 

For  the  anticipated  PWT  facilities  tests,  we  expect  to  compare 
the  ultrasonic  readings  with  venturi  or  other  standard  measures  of  flow 
velocity  and  mass  flow  rate.  Such  reference  measures  are  understood 
to  be  available  at  AEDC  for  use  in  the  PWT  facilities  tests. 
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DATA  ACQUISITION  AND  RECORDING 


At  present,  the  5  MHz  liquid  fuel  "steady  state"  flow  velocimeter 
system  we  are  developing  for  the  Eustis  Directorate  provides  a  digital 
readout  of  mass  flow  rate  M.  (This  system  contains  several  improve¬ 
ments  over  that  described  in  the  appended  1972  publication.  )  By  rela¬ 
tively  simple  changes  in  or  additions  to  the  basic  program,  the  calculator 
chip  in  this  system  can  yield  digital  computation  results  for  other  param¬ 
eters  such  as  c,  p  v,  v/c,  etc.  Totalizing,  M  =  2M(t)  At,  can  also  be 
accomplished,  which  is  convenient  for  weigh  tank  calibrations  and  other 
purposes.  Conversion  to  analog  readout  can,  of  course,  be  accomplished 
with  D/A  modules.  The  optimized  low  frequency  system,  designed  to 
operate  at,  say,  a  50  kHz  carrier  frequency,  could  be  expected  to  provide 
similar  outputs,  digital  and/or  analog.  The  digital  response  time  of  the 
5  MHz  system  presently  is  ~  1  second.  This  could  be  shortened  to  ~0.  1 
second  if  required. 

In  another  design  developed  in  our  laboratory,  unsteady  flow  has 
been  measured  with  a  dynamic  flowmeter  having  a  response  time  of  ~10 
ms.  One  form  of  this  system  is  described  in  the  1973  appended  publica¬ 
tion.  It  is  seen  that  a  continuous  analog  signal  is  available  for  recording. 
For  example,  in  oscillatory  flow  tests  on  liquid  nitrogen,  conducted 
February  1974,  response  out  to  50  Hz  was  demonstrated  by  recording 
output  voltages  using  a  light  beam  oscillograph,  and  also  using  an  oscil¬ 
loscope  and  camera.  Similarly,  a  chart  recorder  and  oscilloscope  and 
camera  were  used  to  record  unsteady  or  transient  ultrasonic  flow  mea¬ 
surements,  wherein  a  column  of  water,  approximately  75  mm  diameter  x 
nearly  2  m  long,  was  allowed  to  drop  under  the  acceleration  of  gravity,  for 
intervals  up  to  ~0.  5  seconds. 

It  is  to  be  understood  that  gas  flowmeters  which  may  be  developed 
for  the  PWT  facilities  or  other  AEDC  facilities,  could  be  built  so  their 
output  would  be  compatible  with  the  particular  recording  equipment  avail¬ 
able  at  the  test  site  of  interest. 
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EXPERIMENTAL  RESULTS 

Limited  bench  tests  were  carried  out  to  obtain  experimental 
information  on  the  statistics  of  signals  transmitted  across  a  25  cm 
(10  inch)  path  in  turbulent  gases.  The  precision  of  ultrasonic  flow 
velocity  measurements  depends  on  maximizing  signal-to-noise>  and 
processing  out  the  phase  and  amplitude  jitter  as  well  as  the  spectral 
broadening  of  the  transmitted  signal.  The  effectiveness  of  coherent 
detection  in  improving  the  precision  of  such  measurements  was  dem¬ 
onstrated.  Finally,  electrostatic  transducers,  particularly  suitable 
to  the  present  work,  were  tested.  The  transducers  were  mounted 
nonintrusively  with  minor  disturbance  to  the  wind  tunnel  wall.  The 
resulting  transducers  have  excellent  isolation  from  the  wall,  can  be 
driven  over  a  broad  frequency  range,  and  can  match  power  into  the 
gas  more  efficiently  than  other  high  frequency  transducers.  Thus, 
key  ideas  for  an  ultrasonic  flowmeter  have  been  reduced  to  practice 
and  serve  as  a  basis  for  extrapolating  to  higher  Mach  number 
conditions. 

A  simple  wind  tunnel  (Fig.  22)  was  constructed  to  create  a 
path  in  a  turbulent  air  stream  25.  4  cm  long,  flowing  at  speeds  be¬ 
tween  Ms  0.  1  and  0.  2.  A  standard  15  HP  blower  was  used  to  take 
air  from  outside  the  building,  drive  it  through  7  feet  of  10  inch  square 
cross  section  wind  tunnel  duct,  and  exit  the  air  back  into  the  atmo¬ 
sphere.  The  test  section  was  mounted  at  the  end  of  a  3  foot  length 
of  1  0  inch  cross  section  duct  work  attached  to  the  blower  as  shown 
in  Fig.  22.  The  dynamic  pressure  and  air  temperature  were  deter¬ 
mined  6  inches  downstream  from  the  test  section. 

Since  the  test  section  is  only  four  duct  diameters  downstream 
from  the  entrance,  the  boundary  layer  is  only  about  1  inch  thick  and 
the  freestream  turbulence  is  determined  by  blower  turbulence  at  the 
duct  entrance.  It  was  felt  that  any  attempt  to  control  or  measure 
freestream  turbulence  was  beyond  the  scope  of  the  limited  tests  en¬ 
visioned.  However,  some  speculation  on  the  state  of  the  flow  at  the 
duct  inlet  might  be  appropriate.  The  blower  is  designed  to  avoid 
large  scale  turbulence  which  causes  equipment  vibration.  However, 
a  considerable  volume  of  turbulent  gas  is  generated  by  the  boundary 
layer  on  the  blower  housing.  The  outer  rim  of  the  blower  blade  is 
traveling  roughly  at  M  0.  25  relative  to  the  gas.  On  the  8  ft  of  blower 
housing  between  cutoff  and  the  end  of  the  housing  a  turbulent  boundary 
layer  of  about  3  inches  is  developed.  Typically,  the  velocity  profile 
of  such  a  blower  peaks  at  the  upper  edge  of  this  boundary  layer  and 
drops  to  30%  of  its  peak  value  at  the  blower  wall  near  the  cutoff.  The 
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Fig.  22,  Wind  tunnel  for  bench  tests 
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sharp  constriction,  at  point  A  (Fig,  22)  displaces  this  maximum  towards 
the  center  of  the  duct  and  reduces  the  variation  in  the  mean  flow  across 
the  duct  entrance.  The  discontinuity  at  flange  A  also  causes  separation 
of  the  flow  and  subsequent  thickening  of  the  turbulent  boundary  layer 
downstream.  Thus,  the  flow  in  the  duct  freestream  is  highly  turbulent. 
A  large  cross  flow  component  is  present  due  to  the  right  angle  turn  at 
the  blower  intake.  In  addition,  there  is  a  wake-like  periodicity  (period 
'VT  sec)  in  the  freestream  which  is  readily  observed  in  the  fluctuations 
of  the  ultrasonic  transmissions  and  the  dynamic  head  readings.  These 
nearly  periodic  variations  of  flow  amount  to  about  5%  of  the  mean  flow.* 


The  qualitative  flow  description  can  be  translated  into  variations 
of  amplitude  and  phase  of  the  ultrasonic  transmissions  across  the  duct. 

u 


A  5%  fluctuation 


(5)‘ 


in  a  mean  flow  U  of  Mach  number  —  = 

c 


0.  1  is 


equivalent  to  a  sound  speed  variation  of  0.  5%.  That  is,  when  the  ve¬ 
locity  fluctuation  u  is  in  the  same  direction  as  the  propagation  direction, 
the  sound  speed  is  increased  from  c  to  c  +  u.  Similarly,  an  opposing 
velocity  fluctuation  decreases  the  sound  speed  to  c  -  u.  Here  c  is  the 
sound  speed  in  a  still  gas.  Assuming  the  largest  vortex  to  be  ~10  cm, 
the  magnitude  of  the  phase  fluctuations  caused  by  these  fluctuating  flows 
can  be  estimated.  The  maximum  possible  phase  shift  fluctuation  can  be 
obtained  by  multiplying  the  fluctuation  in  sound  speed  times  the  number 
wavelengths  in  the  path  times  360  .  The  results  are  given  in  Table  5 
for  room  temperature  air  flow  in  the  present  wind  tunnel. 


Table  5.  Propagation  parameters  for  25.  4  cm  path  in 
wind  tunnel. 


Frequency 

kHz 

Number  of  Wavelengths 
in  25.  4  cm 

A  6 

max 

Wavelength 

cm 

150 

114 

205 

0.  233 

100 

75.  8 

13.6 

0.  335 

50 

37.  9 

68.  2 

0.  670 

10 

7.  58 

13.  6 

8.  35 

An  additional  (~0.  1%)  index  of  refraction  fluctuation  results  from 
the  convection  of  air  from  the  warm  boundary  layers  into  the  main  stream 


*It  should  be  noted  here  that  no  periodicity  is  observed  in  the  noise  which 
is  dominated  by  the  shear  in  the  wall  boundary  layer. 
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by  secondary  flows.  In  fact,  the  largest  phase  shift  fluctuation  observed 
at  50  kHz  was  30  .  As  can  be  seen  by  the  oscillogram  in  Fig.  25a,  the 
average  fluctuation  was  less  than  10  at  50  kHz.  The  measured  phase 
fluctuations  can  be  used  to  characterize  the  size  and  intensity  of  the  free- 
stream  vorticity  in  some  detail  (Schmidt  and  Tilmann,  1969,  1972;  G.  Witten, 
1971;  Dinkelacker  and  Stiewitt,  1971).  However,  it  is  sufficient  here  to  note 
that  the  amplitude  and  phase  jitter  of  sound  transmitted  through  the  present 
wind  tunnel  are  largely  due  to  the  freestream  flow  fluctuations  discussed 
above. 


The  acoustic  noise  is  determined  by  the  shear  on  the  wind  tunnel 
wall.  The  effects  of  the  wall  shear  should  extend  into  the  flow  a  distance 
roughly  equal  to  the  equivalent  flat  plate  boundary  layer  thickness  (for  a 
completely  turbulent  freestream) 

6  =  0.37,  (v)1/5- 

Taking  f  to  be  the  distance  from  the  duct  entrance  to  the  test  section, 

U  the  freestream  flow  velocity,  and  v  the  kinematic  viscosity  for  300°K 
air,  the  boundary  layer  thickness  is  about  5  cm  (2  inch).  Thus,  the 
boundary  layers  extend  over  roughly  40%  the  total  duct  cross  section. 

The  slope  of  the  noise  versus  frequency  allows  the  displacement  thick¬ 
ness  to  be  estimated.  As  described  below,  this  was  carried  out  and  gave 
the  expected  value  within  ~2%.  This  shows  that  the  boundary  layer  does 
indeed  dominate  the  flow  noise. 

Two  test  sections,  with  10  inch  square  cross  sections  to  match 
the  duct,  were  made.  The  test  section  was  joined  to  the  duct  as  indicated 
in  Fig.  22  by  sheet  metal  flanges.  Grooves  ~l/4"  wide  and  1/8"  deep 
were  formed  on  the  inside  of  the  duct  by  the  radius  of  curvature  of  joint 
between  the  flanges  and  the  duct  side  walls.  This  disturbance  in  the  duct 
wall  is  much  smaller  than  the  turbulent  boundary  layer.  In  addition,  no 
difference  inflow  noise  level  was  observed  for  transducers  located  differ¬ 
ent  distances  from  this  groove. 

Electrostatic  transducers  were  "flush"  mounted  on  the  duct  walls. 
A  construction  was  employed  which  emphasized  the  adaptability  of  the 
electrostatic  transducers  for  model  instrumentation.  The  transducers 
are  completely  flexible  and  can  be  mounted  as  easily  on  curved  surfaces 
as  on  flat  surfaces.  The  transducers  match  acoustically  into  the  gas 
better  than  into  the  solid  walls.  Therefore,  no  acoustic  feedthrough  in 
the  walls  is  observed.  The  overall  height  of  the  transducers  above  the 
wind  tunnel  walls  is  0.  01 6"  for  the  one  of  the  design  employed  and  0.  004" 
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(0.  1  mm)  for  the  other.  The  major  perturbation  is  the  coaxial  line 
required  to  carry  the  signal.  Even  this  could  have  been  done  with 
strip  line  construction  without  broaching  the  wind  tunnel  walls. 

The  detailed  construction  of  the  transducers  is  given  in  Figs. 

23a  and  b.  A  large  diameter  transducer  was  built  with  paper  spacers 
to  allow  operation  at  low  pressures.  The  displacement  required  to 
produce  a  given  pressure  p^  in  the  gas  is  related  to  the  frequency 
according  to  . 

i  =  -a- 

pcu 

where  p  and  c  are  density  and  sound  speed  of  the  gas  and  u  is  the 
angular  frequency  (2  it  f).  Thus,  a  low  frequency  transducer  needs  a 
large  amplitude  swing  to  create  a  given  pressure.  To  maintain  a  high 
output  p"  at  low  frequencies  requires  sufficient  spacing  between  the 
membranes  to  accommodate  $  .  The  paper  washers  provide  the  spac¬ 
ing  to  maintain  output  down  to  9  kHz. 

The  transducer  consists  of  two  aluminized  mylar  layers.  The 
aluminum  coating  on  the  middle  layer  is  the  high  voltage  electrode.  A 
dc  voltage  (200V)  on  the  high  voltage  electrode  would  cause  the  entire 
assembly  to  draw  together.  Both  the  aluminum  film  and  the  outer  mylar 
film  and  the  wind  tunnel' wall  are  grounded.  Superposition  of  an  alter¬ 
nating  field  (50V)  on  the  membrane  causes  the  two  mylar  films  to  vibrate 
and  emit  an  ultrasonic  wave.  The  connections  to  the  aluminum  films  are 
accomplished  by  silver  epoxy  cement  and  the  paper  spacers  are  bonded 
to  the  middle  mylar  film.'  The  aluminum  is  removed  from  the  outer  1/4" 
of  the  middle  mylar  to  prevent  short  circuiting.  The  outer  rim  of  both 
mylars  are  bonded  with  a  nonconducting  cement  to  the  wind  tunnel  walls. 
The  active  area  of  the  larger  transducer  is  roughly  equal  to  the  area  of 
the  aluminum  film  (31.  6  cm^)  on  the  middle  mylar.  Referring^to  Table  5 
and  Fig.  23c,  the  ratio  of  transducer  diameter  to  wavelength  varies 

between  1  and  30  for  the  larger  transducer. 

The  small  transducer  used  in  the  present  work  does  not  have  any 
spacers.  This  transducer  operates  on  the  air  trapped  in  surface  rough¬ 
ness.  The  roughness  of  the  wall  and  the  middle  mylar  aluminum  coat  is 
exaggerated  to  illustrate  this  point  in  Fig.  23b.  The  output  of  the  smaller 
transducer  falls  off  at  about  40  kHz.  The  active  area  of  this  transducer 
is  2.  54  cm  (1  inch)  giving  a(^-  )  of  ~0.  3-10.  As  discussed  below,  the 
small  j  transducer  was  used  to  provide  acoustic  rays  over  a  large  an¬ 
gular  spread. 
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LEGEND: 

[  1  ALUMINUM  COATING 


MYLAR  DIELECTRIC 
PAPER  WASHERS  * 


SILVER  EPOXY 


Fig.  23a.  Low  cutoff  frequency  large  diameter  electrostatic 
transducer. 


COAXIAL  CONNECTOR 


Fig.  23b.  High  cutoff  frequency  small  diameter  electrostatic 
transducer. 
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Transducers  were  mounted  on  opposing  walls  as  shown  in  the 
wind  tunnel  diagram,  Fig.  22.  Two  sets  of  large  diameter  (Fig.  23a) 
transducers  were  mounted  on  the  vertical  side  walls.  One  pair  was 
mounted  on  walls  directly  opposing  each  other* **  (x  =  0).  Thus,  when 
the  flow  was  turned  on  the  signal  dropped.  A  second  pair  of  trans¬ 
ducers  was  mounted  with  an  offset  of  x  =  1  inch.  The  beam  drift  in 
the  present  wind  tunnel  flow  is  roughly  one  inch  (see  beam  drift  above). 
These  transducers  provide  isignal-to-noise  information  appropriate  to 
a  beam  drift  flow  velocity  measurement  (page  10). 

A  second  set  of  small  (Fig.  23b)  transducers  was  mounted  on  the 
horizontal  walls  with  a  separation  x  =  3".  The  ^  j  of  these  transducers 
was  made  sufficiently  small  to  assure  an  upstream  ray.  The  signal-to- 
noise  information  obtained  with  these  transducers  is  representative  of 
transmission  (upstream-downstream  phase  shift)  measurement  of  flow 
velocity  (page  il  6). 

The  coherent  detection  system  used  for  the  tests  is  shown  in 
Fig.  24.  The  main  burden  of  the  bench  test  electronics  is  carried  by  the 
PAR*'1'  phase-lock  amplifier  Model  HR-8.  All  of  the  components  within 
the  dotted  lines  of  Fig.  24  are  contained  in  the  PAR  HR- 8:  A  cw  sine 
wave  oscillator  is  amplified  and  used  to  drive  an  electrostatic  transducer. 
The  signal  was  picked  up  by  a  similar  receive  transducer,  amplified  and 
fed  back  into  the  PAR  HR-8  signal  channel.  The  first  stage  of  the  HR- 8 
is  a  tuned  amplifier  gang  tuned  with  the  cw  sine  wave  generator  with  a  Q 
of  25.  The  output  of  this  stage  was  detected  and  recorded  to  give  signal 
amplitude.  The  phase-lock  output  was  also  recorded.  This  output  al¬ 
lowed  the  product  of  amplitude  and  angle  to  be  recorded  at  various  noise 
bandwidths  between  1  kHz  and  0.  1  Hz  with  and  without  flow. 

•  The  signal -to -noise  was  determined  by  measuring  the  amplitude 
and  phase-lock  outputs  with  and  without  a  signal  on  the  drive  transducer. 
The  signal  amplitude  to  noise  amplitude  showed  a  100:1  improvement 
when  the  frequency  was  raised  from  10  kHz  to  100  kHz,  as  expected. 

The  tradeoff  between  bandwidth  and  signal-to-nois  e  was  dramatically 
demonstrated  using  the  phase-lock  output. 

The  frequency  and  flow  velocity  scaling  of  the  acoustic  flow  noise 
and  turbulent  scattering  is  well-known  from  recent  work  on  sonar  and 
atmospheric  propagation  studies.  These  scaling  laws  were  used  to  pre¬ 
dict  system  performance  across  a  25  cm  duct  at  elevated  Mach  numbers. 
The  detailed  results  are  presented  below. 


*x  is  defined  in  Fig.  22  as  the  streamwise  distance  between  the  centers 
of  the  transducers. 

**PAR  -  Princeton  Applied  Research,  Princeton,  New  Jersey. 


69 


Fig.  24.  Coherent  detection  system  for  bench  tests 
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The  measured  transmission  showed  maxima  and  minima  every 
200  Hz  as  the  frequency  was  varied.  This  is  due  to  multiple  reflections 
within  the  duct.  There  ape  only  certain  modes  of  sound  propagation 
which  can  satisfy  the  boundary  conditions.  These  modes  can  be  thought 
of  as  systems  of  plane  waves  propagating  at  small  angles  with  respect 
to  the  duct  walls  (Reschevkin,  1963).  The  dominant  mode  is  plane  wave 
propagation  normal  to  the  duct  walls.  Other  modes  reflect  back  and 
forth  at  angles.  The  observed  maxima  are  very  closely  spaced  at  high 
frequencies,  150  kHz,  because  there  are  many  low  angle  modes  to  re¬ 
verberate.  The  modes  begin  to  disappear  at  about  250  kHz  because  of 
absorption  and  overlap. 

There  are  two  techniques  for  avoiding  the  amplitude  buildup  due 
to  the  reverberation  process  described  above.  First,  the  system  can 
be  pulsed  and  reverberations  gated  out.  Second,  the  operating  frequency 
may  be  chosen  sufficiently  high  to  absorb  most  of  the  energy  (under  flow 
conditions)  in  one  round  trip. 

The  only  effect  the  reverberations  had  during  the  present  tests 
was  to  increase  the  phase  jitter.  The  flow  variations  disturbed  the 
resonance  of  the  mode,  causing  large  phase  shifts  (Mungar,  1969; 

Eversman,  1972). 

The  signals  transmitted  across  the  gas  path  for  the  small  trans¬ 
ducer  at  50  kHz  are  shown  in  Figs.  25  and  26.  The  signal  and  noise  under 
flow  conditions  are  compared  with  the  no  flow  signal  in  Fig.  25  for  the 
case  where  the  transmitter  is  located  upstream.  The  signal  transmitted 
under  no  flow  is  smaller  than  the  signal  with  the  flow  on  because  the  small 
transducers  are  offset  by  7.  62  cm.  The  small  transducers  are  ~1.  90X  in 
radius  at  50  kHz.  Referring  to  Fig.  14,  the  angle  at  which  transducer 
output  is  reduced  by  one-half  is  0.  132  radians.  The  angle  between  the 
centers  of  the  two  transducers  is  0.  253  radians.  Thus,  under  no  flow  one 
would  expect  a  signal  level  reduced  by  about  0.  23  from  the  maximum.  Under 
flow  conditions  an  upstream  ray  offset  angle  equal  to  the  Mach  number  (0.  12 
radians)  reaches  the  receive  transducer.  The  signal  at  0.  12  radians  off 
beam  center  should  be  about  0.  580  of  maximum.  Thus,  the  ratio  of  the  flow 
to  no  flow  signals"'  is 


The  ratio  observed  from  Fig.  25  is 

R  =  3.  66  . 

This  indicates  virtually  no  attenuation. 


*The  square  is  due  to  angle  of  incidence. 
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No  Flow  -  Receive  Signal  -  500  (x  V  Sensitivity  Flow  Signal  -  1  mV  Sensitivity 


Noise  Under  Flow  -  50  p.V  Sensitivity  Noise  Under  Flow  -  1  mV  Sensitivity 

Fig.  25.  Transmitted  signals  across  a  25.4  cm  duct  with  M  =  0.  15.  2.  5  cm  transducers,  receiver 

center  7.  5  cm  downstream  from  transmitter  center. 
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Noise  Under  Flow  -  100  (jl  V  Sensitivity  Noise  Under  Flow  -  200  p  V  Sensitivity 

Fig.  26.  Transmitted  signals  across  a  25.4  cm  duct  with  M  =  0.  15.  2.  5  cm  transducers,  receiver 

center  7.  5  cm  upstream  from  transmitter  center. 
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Similar  signals  are  shown  in  Fig.  26  for  the  upstream  direction. 
The  signal  under  no  flow  is  larger  in  this  case  than  the  signal  under 
flow.  The  ray  which  reaches  the  receiver  under  flow  conditions  is  off¬ 
set  by  the  0.  25  3  radians  due  to  transducer  location  plus  0.  15  radians 
due  to  Mach  number  from  the  beam  center.  In  fact,  the  ratio  of  flow 
to  no  flow  signal  from  Fig.  26  is  0.  400. 

The  attenuation  due  to  thermal  absorption  and  scattering  over 
the  25  cm  pathlength  used  is  negligible.  The  attenuation  is  entirely 
due  to  beam  spread.  These  losses  can  be  reduced  by  optimum  deploy¬ 
ment  of  transducers.  For  instance,  an  array  of  downstream  receivers* 
can  be  used  to  give  temperature.  A  receiver  directly  across  from  the 
transmitter  can  be  used  to  receive  the  upstream  ray  containing  flow 
information.  The  "temperature"  receivers  would  be  arranged  so  that 
one  of  them  would  be  close  to  the  beam  drift  position  at  each  Mach 
number.  Thus,  the  attenuation  for  these  receive  positions  would  be 
less  than  observed  in  the  bench  tests.  The  receiver  directly  across 
from  the  transmitter  would  see  a  ray  between  0.  2  and  0.  7  radians  off 
the  transmitter  beam  axis  for  Mach  numbers  between  M  =  0.  2  and 
M  =  0.  8.  Thus,  the  signal  attenuation  for  any  of  these  transducers 
would  be  similar  to  that  obtained  in  the  bench  tests  reported  above. 

The  signal-to-noise  versus  frequency  was  also  measured.  The 
noise  versus  frequency  is  shown  in  Fig.  27  for  the  one  inch  transducers. 
The  noise  falls  off  as  the  frequency  squared  over  the  entire  frequency 
range  10  kHz-150  kHz.  On  the  other  hand,  the  signal  for  the  transducers 
falls  off  much  more  rapidly  as  the  beam  spread  becomes  too  small  to 
provide  a  ray  at  the  necessary  angle.  Measurements  with  the  larger 
transducers  (described  above,  page  68)  offset  by  the  beam  drift  distance, 
showed  a  signal-to-noise  which  increased  as  f^  over  the  entire  range. 

A  signal-to-noise  ratio  of  40  for  the  downstream  transducer 
and  8  for  the  upstream  signal  at  50  kHz  under  the  M  =  0.  15  flow  was 
observed  in  the  present  bench  tests.  Proper  location  of  the  transducers 
should  allow  signals  equal  to  the  downstream"'  transmission  to  be  ob¬ 
tained  at  least  up  to  Mach  0.  5.  Since  the  noise  increases  at  the  Mach 
number  to  the  sixth  power,  a  signal-to-noise  ratio  of  3  x  10”^  would  be 
obtained  at  M  =  0.  5  in  a  25  cm  duct  at  5  0  kHz  and  a  bandwidth  of  2  kHz 
(Q  =  25).  A  bandwidth  of  ~1  Hz  can  be  obtained  in  the  phase-locked 
mode,  giving  a  signal-to-noise  of  600  to  1  at  M  =  0.  5  in  a  25  cm 
diameter  duct. 


^Different  transducers  would  provide  the  signal  for  different  Mach 
numbers. 
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Fig.  27.  Measured  signal -to-noise  versus  frequency. 
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The  signal  levels  used  in  the  present  work  were  10  atmo¬ 
spheres.  This  was  done  to  allow  a  wide  band  of  frequencies  to  be  used 
in  the  present  tests.  There  is  a  tradeoff  between  gain  and  bandwidth. 

The  signal  levels  can  be  raised  two  orders  of  magnitude  (Kuhl,  1954) 
by  redesigning  the  transmitter.  A  signal-to-noise  level  of  200  to  1 
could  be  obtained  this  way. 

Transmissions  were  obtained  over  the  frequency  region  between 
10  and  500  kHz.  Some  of  these  results  are  summarized  in  Appendix  A. 

At  ~50  kHz,  the  observed  signal-to-noise  ratio  and  phase  stability  were 
optimum,  and  the  electrostatic  transducer  used,  worked  well  at  this 
frequency. 

Extrapolation  of  the  experimental  results  obtained  in  the  present 
program  at  higher  Mach  numbers  indicates  that  a  signal-to-noise  ratio 
of  600  to  1  can  be  obtained  over  a  25  cm  path  in  M  =  0.  5  flow  at  one  atm. 

A  series  of  one  inch  diameter  receive  transducers  operated  at  50  kHz 
would  be  suitable  for  the  measurements.  A  bandwidth  of  1  Hz  in  the 
phase-locked  mode  would  give  about  a  600:1  to  1  signal-to-noise  ratio.* 
This  is  based  on  obtaining  a  10"^  atm  output  from  properly  designed 
electrostatic  transducers.  The  signal-to-noise  ratio  would  be  decreased 
to  35:1  at  M  ~  0.  9.  Thus,  there  is  sufficient  signal  amplitude  to  make  a 
flow  measurement  in  a  ten  inch  duct  essentially  over  the  entire  flow  range 
of  interest  as  defined  by  Requirement  2,  page  12. 

For  comparison  purposes,  and  to  illustrate  the  use  of  a  broadband 
piezoelectric  transducer,  we  include  Fig.  28.  This  shows  transmission 
across  a  ~60  cm  air  path  at  elevated  temperature  and  pressure. 


*The  limitations  of  equipment  available  during  the  present  bench  tests, 
prevented  our  taking  advantage  of  a  bandwidth  as  narrow  as  1  Hz. 
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Drive  waveform 


Received  signal  on 
delayed  sweep 


Fig.  28.  Oscillogram  showing  transmission  of  0.  5  MHz  rf  burst 
across  ~60  cm  static  (no  flow)  air  path,  at  ~40  atm  (585  psig)  and 
slightly  elevated  temperature  (50°C).  Measurement  utilized  a 
pair  of  broadband  piezoelectric  transducers.  Transit  time  =  1.  87 
ms.  Sweep,  5ps/div. 
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OPTIMUM  DESIGN  FOR  FUTURE  TESTS 
IN  PWT/F ACILITIES  AT  AEDC 

The  flow  measuring  system  to  be  used  in  the  PWT  facilities 
must  provide  an  overall  accuracy  in  the  range  0.  5%  to  1%.  It  must 
be  capable  of  operation  in  the  presence  of  very  high  ambient  noise 
levels,  Mach  numbers  between  0.  1  and  0.  5,  and,  of  course,  highly 
turbulent  flow  conditions,  in  which  gradients  are  expected  to  occur. 

Before  describing  the  actual  measurement  system,  the  effects 
of  the  above  environmental  factors  will  be  briefly  described. 

The  acoustic  noise  types  to  be  encountered  in  the  measure¬ 
ment  are  of  two  distinct  types.  These  are  (1)  additive  noise,  and 
(2)  modulation  noise.  The  additive  noise  is  that  which  appears  at 
a  receiving  transducer  due  to  acoustic  radiation  from  within  the 
turbulent  gas,  and  due  to  frictional  effects  of  the  gas  flowing  by  the 
transducer  location.  This  noise  is  present  whether  or  not  an  acoustic 
signal  is  being  transmitted  across  the  duct.  As  shown  earlier  in  this 
report,  the  noise  power  per  hertz  decreases  with  frequency  at  high 
frequencies  (above  ~20  kHz)  and  is  relatively  independent  of  frequency 
at  lower  frequencies. 

The  second  noise  type,  which  we  have  named  modulation  noise, 
results  in  the  spreading  of  the  spectrum  of  an  otherwise  narrow  band 
transmitted  acoustic  wave,  and  is  due  to  the  varying  properties  (tem¬ 
poral  and  spatial)  of  the  turbulent  medium  being  probed.  Because  of 
these  perturbations,  multipath  effects  will  give  rise  to  rapid  fading 
(amplitude  modulation)  of  the  received  signal.  In  addition,  local  varia¬ 
tions  in  sound  speed  and  flow  velocity  will  result  in  phase  modulation  of 
the  received  wave.  These  modulation  effects  become  more  severe  as 
the  frequency  of  the  transmitted  wave  is  increased  because  more  acoustic 
wavelengths  must  be  traversed  by  the  transmitted  wave  across  the  duct, 
which  results  in  increased  phase  shifts  as  well  as  increased  multipath 
interference.  If  we  characterize  the  medium  as  being  permeated  with 
moving  scattering  centers  (the  inhomogeneities),  it  is  easy  to  visualize 
how  scattering  by  these  centers  can  result  in  interference  effects  due 
to  small-angle  forward  scattering.  Such  interference  will  increase  with 
frequency. 

In  addition  to  the  above,  acoustic  reverberation  of  the  transmitted 
wave  within  the  duct  can  occur.  Such  reverberation  would  not  only  inher¬ 
ently  reduce  measurement  accuracy,  but  would  also  greatly  enhance  the 
modulation  noise  effects  discussed  above. 
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From  the  above  discussion,  it  is  clear  that,  in  order  to 
overcome  additive  noise  effects,  we  should  attempt  to  operate  at 
high  frequency,  utilizing  very  narrow  bandwidth.  Coherent  opera¬ 
tion  would  be  ideal.  However,  such  high  frequency  operation  would 
result  in  unacceptable  phase  and  amplitude  modulation  of  the  trans¬ 
mitted  wave.  Depending  on  the  pathlength,  use  of  high  frequencies 
may  also  result  in  excessive  attenuation.  Thus,  there  would  appear 
to  be  impasse.  However,  the  flowmetering  situation  can  be  greatly 
improved  if  the  following  features  are  incorporated  in  the  measuring 
system: 

(1)  Utilize  a  high  frequency  carrier  wave  to  provide  operation 
in  the  spectral  region  of  low  additive  noise. 

(2)  Modulate  the  above  carrier  in  such  a  way  as  to  reduce  or 
eliminate  the  effects  due  to  modulation  noise  and  rever¬ 
berations. 

(3)  Simultaneously  transmit  both  upstream  and  downstream  to 
eliminate  first  order  errors  due  to  variations  in  sound  speed 
from  the  measurement.  (It  may  be  necessary  to  use  beam 
steering  to  efficiently  accomplish  this  over  a  common  path.  ) 

The  measurement  technique  described  below  meets  the  three 
requirements. 

A  coherent  carrier  wave  in  the  range  50  kHz  to  100  kHz  is 
amplitude  modulated  by  a  rectangular  pulse  whose  time  duration  is 
approximately  half  the  transit  time  across  the  duct  (on  the  order  of 
half  a  millisecond  for  a  one  foot  duct).  This  waveform  is  applied 
simultaneously  to  both  the  upstream  and  the  downstream  transducers. 

At  an  appropriate  time  (~1  millisecond  after  transmission  begins),  the 
two  transducers  are  connected  to  separate  receivers.  The  receivers 
are  connected  for  ~1  millisecond.  This  sequence  is  repeated  after  a 
time  sufficiently  long  (~10  milliseconds)  to  insure  that  energy  due  to 
reverberations  in  the  duct  will  not  enter  the  receivers.  The  modulation 
waveform  (assuming  a  one  foot  duct  and  sound  speed  c  -  10^  ft/ sec)  is 
shown  in  Fig.  29,  along  with  the  timing  of  the  receiver  gates  and  the 
arrivals  of  the  upstream  and  downstream  waveforms  at  the  appropriate 
receivers.  (A  block  diagram  is  shown  in  Fig.  30.  ) 

-4 

In  the  timing  scheme  of  Fig.  29  a  modulation  period  of  5  x  10 
sec  and  utilization  of  a  50  kHz  carrier  provide  25  complete  cycles  (not 
shown  in  the  figure)  in  each  transmission.  The  total  transmitted  wave¬ 
form  can  be  mathematically  shown  to  be  the  Fourier  series  of  the 
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Fig.  29.  Timing  diagram  showing 

(a)  Simultaneous  upstream/ downstream  transmission 

(b)  Simultaneous  enabling  of  upstream/ downstream  receivers 

(c)  Downstream  received  signal.  Delay  = 

(d)  Upstream  received  signal.  Delay  = 
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modulation  is  just  the  reciprocal  of  the  repetition  rate  (100  Hz  in  the 
example).  This  frequency  can  be  recovered  in  the  receivers  by  means 
of  amplitude  demodulation  and  narrow  band  filtering.  The  times  Tj 
and  T ^  (see  Fig.  30)  can  be  independently  determined  by  means  of 
separately  measuring  the  phase  shift  of  each  of  the  100  Hz  fundamentals 
in  the  upstream  and  downstream  receivers.  It  can  be  shown  that  the 
flow  velocity,  averaged  over  the  path,  is  given  by 


where  l  =  axial  distance  between  transducers.  Thus,  accurate  mea¬ 
surements  of  both  T2  and  Tj  insure  an  accurate  determination  of  v. 

This  measurement  system  is  considered  to  be  the  optimum 
design  for  meeting  those  contractually-specified  parts  of  Requirement 
2,  that  can  be  simulated  in  the  PWT  facilities.  It  is  a  low-frequency 
(e.  g.  ,  50  kHz)  version  of  a  5  MHz  flow  velocimeter  that  has  been  built 
for  the  Eustis  Directorate  for  use  with  aviation  fuels.  Based  on  pre¬ 
liminary  tests  with  the  5  MHz  system,  and  in  view  of  the  analysis  and 
test  results  in  the  present  phase  of  this  program,  it  appears  that  the 
optimum  design  will  provide  a  practical  solution  for  the  PWT  facilities' 
and  similar  applications. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

In  view  of  Requirement  2  in  the  Scope  of  Work,  and  in  parti¬ 
cular,  with  regard  to  tests  anticipated  in  the  PWT  facilities  at  AEDC, 
it  is  concluded  that  the  transmission  approach  presently  offers  the 
most  promise*  with  regard  to  accuracy  and  response  time,  for  flows 
up  to  at  least  Mach  0.  5,  and  possibly  as  high  as  ~0.  9.  Theory  and 
experiments  indicate  that,  for  a  path  of  ~25  cm,  a  suitable  carrier 
frequency  is  50  kHz.  Based  on  the  present  requirements,  and  taking 
note  of  recent  developments  using  coherent  detection  of  modulated  cw 
transmission  across  liquids  such  as  aviation  fuel,  water  and  liquid 
nitrogen,  it  appears  reasonable  to  apply  a  low-frequency  (e.  g. ,  50 
kHz)  version  of  that  type  of  flow  measuring  system  to  the  gas  flow¬ 
metering  situation.  Electrostatic  transducers  offer  advantages  of 
impedance  matching  and  small  perturbation  of  the  duct,  and  so  may 
be  preferred  over  the  more  usual  piezoelectric  types,  or  even  over 
the  more  rugged  magnetostrictive  types  which  staffer  from  lower 
sensitivity. 

It  is  recommended  that  the  "optimized  design"  outlined  in  this 
report,  be  fabricated  and  tested  in  the  PWT  facilities  at  AEDC.  The 
object  of  such  tests  would  be  to  determine,  under  flow  conditions  of 
interest,  the  actual  limits  of  signal/noise  ratio,  and,  correspondingly, 
the  limits  on  accuracy  and  response  time.  It  is  expected  that  accuracy 
of  ~0.  5%,  and  response  times  of  <  1  second,  can  be  achieved. 


*That  is  to  say,  the  ultrasonic  transmission  approach,  preferably  executed 
using  simultaneous  launching  of  upstream  and  downstream  waves  over  a 
common  path,  is  to  be  preferred  over  the  other  ultrasonic  approaches  con¬ 
sidered,  namely,  reflection/ doppler  scatter,  beam  drift,  correlation,  noise 
and  vortex  shedding. 
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APPENDIX  A 

ACOUSTIC  TRANSMISSION  ACROSS  25  CM  DUCT  AT  MACH  0.  15 

The  propagation  of  sound  waves  through  a  Mach  0.  15  flow 
was  investigated  over  a  wide  range  of  frequencies  and  with  several 
transducer  configurations  as  indicated  in  the  text.  The  parameters 
which  led  to  the  choice  of  50  kHz  as  our  operating  frequency  can  be 
understood  more  clearly  by  studying  signal  and  noise  parameters 
over  a  wide  range  of  conditions.  The  four  major  test  configurations 
are  displayed  in  Fig.  Al.  The  direction  of  flow,  the  boundary  layer 
thickness,  and  transmission  direction  are  indicated  in  the  figure. 

d  ^ 

Also  shown  is  the  near  field  limit  —  for  frequencies  at  which  data  was 

taken.  Roughly,  the  sound  beam  propagates  with  no  spreading  out  to 
the  near  field  limit.  Beyond  the  near  field  limit  spreading  occurs 
and  the  intensity  is  a  function  of  the  angle  of  propagation.  The  near 
field  is  shown  by  parallel  lines  in  Ifig.  Al  to  indicate  the  parallel 
phase  fronts.  The  far  field  pattern  is  indicated  by  an  intensity  versus 
angle  pattern.  The  intensity  is  proportional  to  the  radius  vector  drawn 
from  the  center  of  the  transducer  to  the  far  field  pattern. 

The  first  set  of  photographs  (A2-A5)  shows  a  small  transducer 
transmitting  to  a  downstream  transducer  (upper  left  in  Fig.  Al).  The 
flow  offsets  the  beam  in  such  a  way  as  to  increase  the  received  signal 
as  compared  to  the  no  flow  signal.  The  second  set  of  photographs 
(A6-A9)  shows  a  small  transducer  transmitting  to  an  upstream  trans¬ 
ducer  (upper  right  in  Fig.  Al).  The  flow  moves  the  beam  in  such  a 
.way  as  to  decrease  the  received  signal  as  compared  to  the  no  flow  sig¬ 
nal.  The  no  flow  signal  with  flow  and  flow  noise  is  shown  for  each 
frequency.  The  flow  noise  was  obtained  by  turning  off  the  transmit 
signal  and  recording  the  receive  output.  The  receiver  noise  with  no 
flow  was  orders  of  magnitude  below  the  flow  noise  displayed  in  the 
photographs. 

The  signal -to- noise  for  the  upstream  transmitter  series  is 
less  than  1  at  10  kHz,  increasing  to  about  50  at  40  kHz  and  reaching 
~150  at  80  kHz.  However,  the  phase  coherence  of  the  signal  is  optimum 
at  ~50  kHz.  The  jitter  in-phase  in  Fig.  A4  (40  kHz)  is  about  45°  while 
in  Fig.  A5  (80  kHz)  the  jitter  increases  to  180  .  The  phase  jitter  in¬ 
crease  in  the  wavelength  becomes  smaller  relative  to  either  the  large 
scale  vorticity  or  the  receiver  transducer  diameter. 

The  signal-to-noise  ratio  and  phase  jitter  show  a  similar  pattern 
for  the  downstream  transmitter  series.  However,  the  signal-to-noise 
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ratio  is  4  at  40  kHz  and  5  at  80  kHz.  This  is  due  to  the  fact  that  the 
flow  pushes  the  beam  further  away  from  the  transducer.  The  signal- 
to-noise  ratio  does  not  improve  as  quickly  as  one  would  expect  with 
frequency  because  at  80  kHz  the  near  field  extends  about  75%  of  the 
distance  across  the  tube.  Thus,  the  rays  sufficiently  angled  upstream 
to  reach  the  receiver  have  a  very  low  intensity  at  80  kHz  compared  to 
50  kHz. 


The  final  set  of  photographs  shows  signals  obtained  with  large 
transducers  offset  by  the  distance  the  beam  is  expected  to  drift  due  to 
flow,  2.  54  cm.  The  geometry  for  transmission  with  (lower  left)  and 
against  (lower  right)  the  flow  is  shown  in  Fig.  Al.  Note  the  signal -to- 
noise  is  3  at  10  kHz  for  the  large  transducers  since  we  are  operating 
very  close  to  the  beam  axis.  In  addition,  the  variation  of  signal  ampli¬ 
tude  with  flow  direction  is  very  small  at  low  frequencies  since  beam 
spread  makes  up  for  beam  drift.  However,  at  higher  frequencies  the 
beam  is  narrow  and  the  downstream  signal  is  larger  than  the  upstream 
signal. 


Over  a  hundred  photographs  of  acoustic  transmissions  were 
taken,  in  addition  to  recorder  traces  of  frequency  response  and  signal - 
to-noise  ratio.  Since  the  absorption  over  large  acoustic  paths  (>  10  cm) 
drastically  reduces  the  signal  above  500  kHz,  we  generally  concentrated 
on  frequencies  below  500  kHz.  In  addition,  for  turbulent  flow  large 
pathlengths  include  large  vortices.  -This  fact,  together  with  the  large 
number  of  wavelengths  in  the  path,  lead  one  to  expect  phase  jitter.  In 
practice,  at  100  kHz  across  the  25.  4  cm  duct,  180°  of  phase  jitter  is 
clearly  observable.  However,  the  signal-to-noise  ratio  is  relatively 

high.  At  lower  frequencies,  ~50  kHz,  the  phase  jitter  is  much  smaller 
than  the  100  kHz  result,  but  the  signal-to-noise  ratio  is  poorer.  Since 
phase  jitter  (or  modulation  noise)  is  undesirable  for  the  proposed  flow¬ 
meter  application,  the  50  kHz  operating  frequency  is  preferred  to.  In 
smaller  ducts  a  higher  frequency  should  be  used  to  maximize  signal- 
to-noise  ratio,  consistent  with  not  exceeding  allowable  phase  jitter. 

Based  on  the  observed  propagation  over  the  wide  range  of  fre¬ 
quencies  obtained  in  the  present  tests,  together  with  our  previous  work 
at  megahertz  frequencies,  we  were  able  to  select  the  optimized  measure¬ 
ment  parameters  for  the  25  cm  duct  flowmeter  problem.  Furthermore, 
we  are  able  to  recommend  the  system  described  on  pages  70-74.  This 
system  is  expected  to  achieve ~  1%  accuracy  and  ~  1  second  response 
time,  for  the  PWT  or  equivalent  facilities. 
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1 00  p  V/cm 
No  Flow  Signal 


100  p  V/  cm 
Flow  Signal 


1 00  p  V/cm 
Flow  Noise 


Fig.  A.  2.  2.  54  cm  transducers  -  transmitter  upstream  by  7.  62  cm 

from  receiver  (10  kHz). 
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100  p  V  /  cm 
No  Flow  Signal 


100  p  V  /  cm 
Flow  Signal 


Fig.  A.  3.  2.  54  cm  transducers  -  transmitter  upstream  by  7.  62  cm 

from  receiver  (20  kHz). 
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250  n  V/cm 
No  Flow  Signal 


500  p  V/cm 
Flow  Signal 


5  0  p  V  /  cm 
Flow  Noise 


Fig.  A.  4.  2.  54  cm  transducers  -  transmitter  upstream  by  7.  62  cm 

from  receiver  (40  kHz). 
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No  Flow  Signal 


250  jji  V /cm 
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1 0  |x  V  /  cm 
Flow  Noise 


Fig.  A.  5.  2.  54  cm  transducers  -  transmitter  upstream  by  7.  62  cm 

from  receiver  (80  kHz). 
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25  |jl  V/ cm 
No  Flow  Signal 


250  (jl  V/cm 
Flow  Signal 


Fig.  A.  6.  2.  54  cm  transducers  -  transmitter  downstream  by  7.  62  cm 

from  receiver  (10  kHz). 
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250  p  V  /cm 
No  Flow  Signal 


100  p  V/cm 
Flow  Signal 


Fig.  A.  7.  2.  54  cm  transducers  -  transmitter  downstream  by  7.  62  cm 

from  receiver  (20  kHz). 
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500  jiV/cm 
No  Flow  Signal 


100  jiV/cm 
Flow  Signal 


100  jiV/cm 
Flow  Noise 


Fig.  A.  8.  2.  54  cm  transducers  -  transmitter  downstream  by  7.  62  cm 

from  receiver  (40  kHz). 
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500  |X  V/cm 
No  Flow  Signal 


20  |x  V/cm 
Flow  Noise 


Fig.  A.  9.  2.  54  cm  transducers  -  transmitter  downstream  by  7.  62  cm 

from  receiver  (80  kHz). 
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5  mV/ cm 
No  Flow  Signal 


5  mV/ cm 
Flow  Signal 


5  mV/ cm 
Flow  Noise 


Fig.  A.  10.  6.  35  cm  transducers  -  transmitter  upstream  by  2.  54  cm 

from  receiver  (10  kHz). 
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5  mV/ cm 
No  Flow  Signal 


5  mV/ cm 
Flow  Signal 


S2  mV /  cm 

Flow  Noise 

Fig.  A.  11.  6.  35  cm  transducers  -  transmitter  upstream  by  2.  54  cm 

from  receiver  (20  kHz). 
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5  mV  /  cm 
No  Flow  Signal 


5  mV/ cm 
Flow  Signal 


1  mV/ cm 
Flow  Noise 


Fig.  A.  12.  6.  25  cm  transducers  -  transmitter  upstream  by  2.  54  cm 

from  receiver  (50  kHz). 
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1  mV  /  cm 
No  Flow  Signal 


50  p.  V/ cm 
Flow  Noise 


Fig.  A.  13.  6.  25  cm  transducers  -  transmitter  upstream  by  2.  54  cm 

from  receiver  (82.  2  kHz). 
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2  mV/ cm 
No  Flow  Signal 


2  mV /cm 
Flow  Signal 


50  (x  V/ cm 
Flow  Noise 


Fig.  A.  14.  6.  25  cm  transducers  -  transmitter  upstream  by  2.  54  cm 

from  receiver  (99.  4  kHz). 
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5  mV/ cm 
No  Flow  Signal 


5  mV/ cm 
Flow  Signal 


5  mV/ cm 
Flow  Noise 


Fig.  A.  15.  6.  25  cm  transducers  -  transmitter  downstream  by  2.  54  cm 

from  receiver  (10  kHz). 
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1 0  mV/ cm 
No  Flow  Signal 


1 0  mV/ cm 
Flow  Signal 


Fig.  A.  1  6.  6.  25  cm  transducers  -  transmitter  downstream  by  2.  54  cm 

from  receiver  (20  kHz). 
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2  mV/ cm 
No  Flow  Signal 


2  mV/ cm 
Flow  Signal 


200  jjl  V/ cm 
Flow  Noise 


Fig.  A.  17.  6.  25  cm  transducers  -  transmitter  downstream  by  2.  54  cm 

from  receiver  (50  kHz). 
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5  mV/cm 
No  Flow  Signal 


2  mV/cm 
Flow  Signal 


50  [iV/cm 
Flow  Noise 


I 
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Fig.  A.  18.  6.  25  cm  transducers  -  transmitter  downstream  by  2.  54  cm 

from  receiver  (82.  1  kHz). 
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2  mV/ cm 
No  Flow  Signal 


Fig.  A.  19.  6.  25  cm  transducers  -  transmitter  downstream  by  2.  54  cm 

from  receiver  (99.  4  kHz). 
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ULTRASONIC  MASS  FLOWMETER 

L.  C.  Lynnworth  and  N.  E.  Pederaen 
Panametrics,  Inc.,  Waltham,  Massachusetts  02154 

ABSTRACT.  Development  began  on  a  new  type  of  ultrasonic  mass  flowmeter  for  fuel  flow  in 
gas  turbine  engines.  The  fuel  flowmeter  consists  of  a  flow  velocimeter,  densitometer,  time  inter- 
valometer,  and  a  metering  section  containing  nonintrusive  transducers.  The  velocimeter  measures 
the  phase  difference  between  two  coherent  individually  coded  5  MHz  waves  continuously  transmitted 
obliquely  upstream  and  downstream  across  a  common  path,  yielding  a  phase  difference  proportional 
to  v/c^,  where  v  =  area-averaged  flow  velocity  and  c  =  sound  speed.  The  densitometer  consists  of 
a  stepped-diameter  probe;  the  difference  in  echo  amplitudes  at  the  wetted  end  and  at  the  dry  step  is 
a  function  of  the  acoustic  impedance  p  c  of  the  fuel  where  p  is  the  fuel  density.  The  time  interval- 
ometer  measures  the  time  T  between  two  successive  echoes,  which  time  is  proportional  to  1/c.  The 
mass  flow  rate  M  is  M  =  (K')(v/c^)  (pc)  (1/T)  where  Kii  a  constant.  The  complete  system  was 
tested  on  various  liquids  at  rates  up  to  ~2300  kg/hr.  It  was  operated  during  and  recalibrated  after 
104-hr  exposure  to  a  contaminated  fluid  flowing  at  ~860  kg/hr.  Response  time  was  20  ms.  The 
flowmeter  can  operate  in  laminar,  transitional  and  turbulent  flow,  and  uses  a  novel  method  of  weight¬ 
ing  the  profile.  It  is  concluded  that  refinement  of  the  present  approach  should  be  able  to  attain  M 
accuracy  of  1%. 

Introduction  no  rotating  parts.  Accordingly,  it  is  not  subject  to 

clogging. 

This  paper  is  mainly  concerned  with  the  ultra¬ 
sonic  measurement  of  the  mass  flow  rate  of  aviation  In  the  following,  we  describe  a  novel  ultrasonic 

fuels,  presently  for  diagnostic  purposes,  and  ulti-  approach  to  this  mass  flowmetry  problem,  with  con- 

mately,  for  control  purposes  too.  Those  who  are  sideration  of  the  flow  profile,  and  we  summarize 

familiar  with  biomedical  applications  only  may  be  tests,  resuits  and  conclusions  based  on  this  first 

surprised  to  hear  that  over  100  ultrasonic  flowmeters  year's  work.  A  more  complete  account  is  available.3 

have  been  used  in  industrial  applications.  Further, 

industrial  test,  measurement  and  control  applications  Ultrasonic  Approach 

of  ultrasonics  cover  over  a  dozen  different  param¬ 
eters,  with  standard  instruments  available  from  many  In  analyzing  ultrasonic  flow  velocimeter  ap- 

sources,  some  of  which  are  listed  in  Table  I.1*2  proaches,  one  finds  that  a  measurement  of  upstream 

minus  downstream  transit  time  (or  phase)  yields  a 
From  the  viewpoint  of  those  responsible  for  term  proportional  to  v/c^,  where  v  *  flow  velocity 

diagnosing  and  controlling  the  performance  of  air-  and  c  =  sound  speed,  whereas  a  measurement  of  the 

craft  engines,  it  is  important  to  know  accurately  the  difference  between  two  frequencies  which  are  recip- 

value  of  the  fuel  mass  flow  rate  M.  The  heating  rocals  of  these  times  responds  to  v,  independent  of  c. 

values  of  typical  fuels  such  as  JF-4,  JP-S,  Diesel'  Doppler  shifts ,  however,  are  proportional  to  v/c. 

fuel  and  Avgas  grades  from  80  to  145  octane  are  with-  These  outputs  may  be  combined  with  Quid  density  (p), 

in  about  1%  of  19000  Btu/lb.  Thus,  the  maximum  characteristic  impedance  (pc)  and  sound  speed  mea- 

power  that  can  be  extracted  from  these  fuels  is  pro-  surements  to  yield  M,  the  mass  flow  rate,  as  shown 

portional  to  M.  in  Table  II. 

In  engine  diagnostic  studies,  one  is  interested  A  two-frequency  Doppler  approach  was  selected 

in  measuring  engine  performance  from  start-up,  at  at  the  beginning  of  this  program,  because  at  that  time 

relatively  low  flow  rates  of  about  20  kg/hr,  up  to  it  appeared  to  offer,  at  least  theoretically,  the  best 

maximum  running  conditions  where  M  =  1000  kg/hr  in  way  of  measuring  average  flow  velocity  with  minimum 
some  engines.  Thus,  a  range  of  ~50:1  is  of  interest.  disturbance  of  simple  pipelines.  Calculations  of  the 

Accuracy  of  better  than  1%  is  desirable  at  high  flow  power  scattered  at  20  MHz  from  fuel  contaminated 

and  ~  5%  at  the  low  end.  For  control,  response  times  according  to  MIL-E-5007C  showed  that  scattering 
of  ~  10-30  ms  are  desired.  Continuing  to  define  the  should  be  adequate  even  in  fuels  containing  ~1%  of 

problem,  the  fuel  composition  may  be  an  unknown  the  specified  contaminant.  However,  tests  on  uncon- 

mixture  of  two  or  more  fuels.  The  fuel  temperature  taminated  liquid  did  not  provide  enough  scatter  to  be 
may  range  between  -55°C  to  +70°C.  The  fuel  may  be  detected  with  our  equipment.  Therefore  we  aban- 
contaminated  by  iron  oxide,  sand  and  dirt  particles  doned  the  Doppler  approach  and  adopted  a  transmission 

having  sizes  up  to  420  pm,  and  by  corrosive  agents.  approach.  Interrogation  was  accomplished  by  trans- 

The  engine  environment  includes  noise  and  vibration.  mitting  a  pair  of  5  MHz  continuous  plane  waves  which 

The  flowmeter  measuring  cell,  nonintrusive  trans-  ideally  remain  undistorted  in  the  flowing  liquid.  These 

ducers  and  electronics  should  be  small,  lightweight  waves  simultaneously  traversed  the  same  path,  and 

and  rugged.  sampled  the  entire  cross  section  of  flowing  liquid 

such  that  v/c^  was  determined.  To  separate  up- 
The  main  problem  with  previous  (turbine)  flow-  stream  from  downstream  waves,  each  wave  was 
meters  has  been  clogging  within  ten  hours  under  con-  phaBe-coded  with  a  pseudo- random-noise  code.  We 
taminated  fuel  flow.  The  present  ultrasonic  flowmeter  think  this  is  the  first  such  separation  in  an  ultrasonic 
employs  recessed,  nonintrusive  transducers,  and  has  flowmeter. 

Paper  L6  IEEE  Ultrasonics  Symposium  Boston  October  4-7,  1972 
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Table  1.  Examples  of  Industrial  process  control  and  quality  control  parameters 
and  vendors  of  standard  ultrasonic  equipment  for  their  measurement. 

Parameter  Vendor 

Flow . Alma  Instruments,  Badger,  Controlotronics,  NUSonics,  Saratoga  Systems,  Scarpa 

Labs,  Tokyo  Precision  Instruments,  Westinghouse. 

Temperature . J.  Bell/Univ.  of  Aston  fwire),  Hewlett-Packard  (quartz),  Honeywell  (fluidic), 

Panametrics  (wire). 

Pressure . Columbia,  Endevco,  Kristal,  PCB  Piezotronics,  Wallace  i<  Tiernan. 

Viscosity . Bendix,  Unipan  Scientific  Instruments. 

Level . Bindicator,  Controlotronics,  Delavan,  Endress  fc  Hauser,  National  Sonics,  Raytheon. 

Density . .  .  ITT-Barton,  National  Sonics. 

Thickness . Automation  Industries,  Branson,  Krautkramer,  Magnaflux,  Nortec,  Panametrics. 

Modulus  .  . . Magnaflux,  H.  M.  Morgan,  Panametrics 

Flaws . See  "Thickness";  and  Holosonics,  Shurtronics,  Sonic  Instruments,  Tehran, 

Ultrasonoscope. 

Fluids  Composition.  .  .  NUSonics,  Tracor. 

Acoustic  Emission  .  .  .  Dunegan,  Endevco,  Nortec,  Panametrics,  Trodyne. 

Leaks . Dawe,  Deleon  (Hewlett-Packard). 

Bubbles . Vero  Precision  Engineering  Ltd. 


Table  II.  Comparison  of  output  products  of  six 

combinations  of  continuous  wave  flow  ve- 
locimeters  and  ultrasonic  densitometers. 


~~~^-~^Dens  itomete  r 

Reflection 

Resonant 

Flow 

Coefficient 

Frequency 

Velocimeter 

■7m 

Doppler(v/c) 

M 

M/c 

Phase  Difference(v/c2) 

M/c 

Frequency  Difference(v) 

Me 

M 

To  determine  fuel  density,  we  utilized  the  re¬ 
flection  of  normally  incident  longitudinal  wave  pulses. 
At  normal  incidence  the  reflection  coefficient  for  these 
plane  waves  depends  only  on  the  characteristic  acous¬ 
tic  impedance  of  the  adjacent  media  at  the  boundary. 
Since  the  fuel's  characteristic  impedance  equals  pc, 
if  the  probe's  Impedance  Zj  is  known,  measurements 
of  the  reflection  coefficient  at  the  boundary,  together 
with  a  measurement  of  sound  speed  c  (or  time  inter¬ 
val  T,  where  T  is  proportional  to  1/c)  enable  one  to 
determine  the  fuel  density  p.  A  separate  transmis¬ 
sion  path  normal  to  the  fuel  flow  direction  was  uti¬ 
lized  to  measure  a  time  interval  T,  which  is  inverse¬ 
ly  proportional  to  c,  and  independent  of  v.  Combining 
these  measurements,  we  have  M  =  (K* )  (v/c2)  (pc) 
(1/T)  where  K'  depends  on  the  area  of  the  cell's  flow 
channel,  and  on  the  units  of  measure. 

Reynolds  Number  and  Flow  Profile 

If  one  interrogates  the  fluid  with  a  wave  that 
averages  over  the  diameter,  but  not  over  the  whole 
area,  one  measures  vj  =  va/K,  where  K  is  less  than 
unity  and  dtpends  on  the  Reynolds  number  Re.  The 
diameter-averaged  flow  velocity  v^  is  larger  than  the 
area-averaged  value  va  by  25%  for  laminar  flow 
(Re  <  2000).  K“  0.  881  +  0.  011  logjQ  Re  for  turbulent 
flow  for  Re  >  104.  Fractional  differences  (v^-Vj)/ 
v j  =  1  •  K  are  listed  in  Table  III  as  functions  of  Re. 

In  many  industrial  situations,  where  flow  is 
turbulent,  temperature  and  composition  of  the  process 
fluid  are  sufficiently  constant  so  that  Reynolds  num¬ 
ber  can  be  estimated  from  va.  Since  K  is  relatively 


Table  III.  Potential  errors  due  to  flow  profile.* 


Re 

<2  x  103 

104 

105 

106 

107 

K 

.750 

.925 

.936 

.947 

.958 

1-K 

25% 

7.5% 

6.  4% 

5.3% 

4.2% 

insensitive  to  Re,  the  estimated  Re  is  sufficient  to 
determine  K  to  a  small  fraction  of  1%,  which  means 
va  can  be  determined  to  that  accuracy  too.  In  our 
case,  however,  the  combination  of  unknown  fuel  mix¬ 
tures  and  wide  temperature  range  prevent  one  from 
determining  K  to  less  than  1%  from  measurements  of 
v^j  alone.  Can  one  measure  va  directly?  Yes,  pro¬ 
vided  one  properly  weights  the  entire  flowing  cross 
section.  To  a  good  approximation,  it  appears  that 
one  way  this  condition  is  met  is  by  using  a  square 
flow  channel,  interrogated  obliquely  in  the  near  field 
of  a  rectangularly-collimated  transducer  (Fig.  la,b). 

Flow  Velocimeter  -  CW  Phase  Measurement 

Similar  to  previous  ultrasonic  transmission 
flowmeters,  the  system  to  be  described  measures  the 
phase  difference  between  waves  transmitted  upstream 
and  downstream.  The  two  cw  waves  traverse  virtually 
the  same  path  at  virtually  the  same  time.  The  two 
waves  have  the  same  carrier  frequency,  yet  can  be 
separated,  or  identified,  by  the  manner  in  which  they 
are  each  modulated  and  then  demodulated. 

In  Fig.  2,  a  5  MHz  crystal  controlled  oscillator 
signal  is  divided  into  two  isolated  channels.  The 
waves  in  the  upper  and  lower  channels  are  phase 
coded  by  means  of  a  dual  coherent  pseudo-random- 
noise  code  generator,  producing  a  sequence  of  l's 
and  0's  such  that  the  resulting  code  displays  all  the 
characteristics  of  perfect  randomness.  The  length  of 
the  code,  as  well  as  the  clock  rate  are  arbitrarily 
variable.  Each  "1"  produces  a  +90°  phase  shift  of 
the  carrier,  while  each  "0"  produces  a  -90°  phase 
shift.  The  codes  in  the  upper  .and  lower  channels  are 
different.  These  pseudo- random  sequences  are 
orthogonal.  The  two  5  MHz  waves  are  independently 
phase  code  modulated  (Cj  and  Cj)  by  the  DBM's  at 
the  left.  These  are  then  fed,  via  circulators,  to  the 
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SQUARE 

FLOW  CHANNEL 


Fig.  1  (a)  Ultrasonic  mass  flowmeter  measuring  cell;  (b)  Transducer  locations. 


two  transducers.  The  received  waves  are  fed  (via  the 
circulators)  to  the  top  and  bottom  DBM's  on  the  right. 
The  demodulating  codes,  C2<j  and  Cld  are  also  fed  to 
these  DBM's.  The  outputs  of  the  DBM's  are  two  cw 
signals,  as  well  as  wide  band  noise.  The  "noise" 
results  from  noncorrelation  of  unwanted  signals  from 
either  the  wrong  channel  or  from  stray  reflections  and 
reverberations  within  the  pipe.  Only  that  wave  which 
is  due  to  a  reverberation-free  straight-through- 
transmission  results  in  a  cw  wave  at  the  demodulator 
DBM's,  We  next  go  through  a  pair  of  tuned  5  MHz 
amplifiers  which  strip  off  most  of  the  nncorrelated 
information,  and  finally  feed  into  a  pair  of  phase- 
locked  loops. 

Instead  of  phase-locked  loops  we  used  a  stan¬ 
dard  phase  meter.  Another  simplification  involved 
using  just  one  code  for  both  directions. 


Tests  on  Fuel  and  Contaminated  Fluid 

The  main  purpose  of  these  tests,  conducted  at 
Avco-Lycoming  between  June  and  August  1972,  was 
to  determine  the  accuracy,  reproducibility  and  lin¬ 
earity  of  the  ultrasonic  mass  flowmeter  system,  and 
to  demonstrate  operation  on  uncontaminated  fuel  such 
as  JP-4  and  JP-5  as  well  as  on  a  fluid  contaminated 
in  accordance  with  MIL-E-5007C.  Response  time  was 
not  determined  in  these  tests,  but  was  calculated  from 
the  system  bandwidth  to  be  20  ms. 

Test  data  for  four  successive  runs  in  JP-5  are 
plotted  in  Fig.  3.  These  data  exhibit  linearity  over 
the  range  indicated.  Significantly,  for  3  of  4  runs, 
the  data  are  within^  4%  of  one  straight  line.  Disre¬ 
garding  the  anomalous  fourth  run,  these  data  are  quite 
encouraging,  in  that  they  demonstrate  the  basic 
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linearity  and  reproducibility  of  the  M  system.  In  a 
later  teat  on  JP-5,  where  hi  ranged  to  ~2300  kg/hr, 
(beyond  the  range  of  the  present  multiplier),  we  found 
'  that  the  measured  phase  angle  could  be  approximated 
as  a  piecewise  linear  function.  Departures  from  lin¬ 
earity  are  apparently  due  mainly  to  eddies  in  the 
measuring  zone.  It  is  expected  that  these  can  be 
minimized  using  a  window,  wedge,  screen  or  other 
means.  Nonideal  inlet  conditions  may  also  contribute 
to  the  observed  nonlinearity.  Another  deficiency  was 
that  undesirable  circuit  unbalances  were  produced  by 
fluctuations  in  fuel  temperature  of  a  few  degrees,  or 
by  changes  in  the  acoustic  radiation  impedance  due  to 
changes  in  impedance  of  the  fuel  or  transducer  back¬ 
ing.  To  overcome  this  problem,  one  should  choose  a 
transducer  with  minimal  variation  in  properties  due 
to  temperature,  and  should  isolate  the  transducer 
from  the  fuel,  using,  for  example,  an  attenuating 
acoustic  window  as  a  pad. 

Following  JP-4  and  JP-S  tests,  the  cell 
(Fig.  la)  was  operated  with  a  contaminated  fluid. 
Exposure  totaled  104  hrs  at  ~860  kg/hr,  plus  16  hrs 
of  filtered  contaminant  at  this  rate,  plus  soaks  at 
no-flow  totaling  428  hrs.  Signal  levels  degraded 
somewhat  during  these  testa,  but  not  due  to  any  fun¬ 
damental  deficiency.  Endurance  appears  to  be  far  in 
excess  of  the  two  2.  5-hour  teats  required  for  instru¬ 
ments  in  MIL-F-3615.  The  complete  system  was 
then  retested  on  JP-4  at  25  C  and  33  C,  and  again 
did  not  exhibit  deficiencies  attributable  to  the  con¬ 
taminant.  Due  to  accuracy  limitations  of  the  phase 
.meter  available  for  these  tests,  performance  of  the 
flowmeter  was  correspondingly  limited.  However, 
phase  meters  are  now  available  which  should  reduce 
this  source  of  errors  to  less  than  0.5%  of  range,  and 
possibly  to  0.25%.  Assuming  comparable  accuracy 
for  the  densitometer,  a  system  accuracy  of  1%  or 
better  for  iJt  appears  attainable. 

Conclusions 

e  A  pseudo- random- noise  phase  coding  technique  has 
been  demonstrated  which  separates  two  continuous 
ultrasonic  waves  of  the  same  frequency  transmitted 
upstream  and  downstream  over  a  common  path. 

This  technique  provided  a  flow  velocimeter 


response  time  of  20  ms,  a  factor  of  25  faster  than 
the  contractual  objective  of  0. 5  s,  and  a  range  up  to 
~  2300  kg/hr. 

e  The  ultrasonic  mass  flowmeter  measuring  cell  has 
been  demonstrated  to  operate  without  clogging,  both 
during  and  after  exposure  to  a  fluid  contaminated 
essentially  in  conformance  to  M1L-E-5007C.  This 
exposure  included  over  100  hr  at  ~860  kg/hr  and 
over  400  hr  soaking  at  zero  flow. 

e  Immunity  to  uncertainties  in  Reynolds  number  or 
flow  profiles  appears  achievable  using  the  square 
flow  channel  wherein  100%  of  the  fuel  cross  section 
is  interrogated  by  a  rectangular  plane  wave  beam  of 
ultrasound.  This  is  the  basis  for  accurate,  linear 
response  to  flow  velocity,  independent  of  fuel  com¬ 
position,  temperature,  density  or  viscosity.  In  the 
present  measuring  cell,  the  observed  departures 
from  linearity  appear  to  be  due  to  nonideal  inlet 
conditions,  probably  correctable  with  a  flow 
fltraightener,  and  to  eddies  at  the  measuring  region, 
probably  correctable  with  flush-mounted  acoustic 
windows,  wedges,  screens  or  other  means. 

e  Based  on  the  encouraging  first  results  obtained  with 
coded  single  fixed  frequency  waves,  it  appears  that 
the  use  of  two  variable  frequencies  should  work  too, 
to  measure  v  independent  of  c.  This  is  the  basis 
for  an  improved  ultrasonic  mass  flowmeter,  where 
both  v  and  p  are  each  measured  independent  of  c. 
The  product  of  their  outputs  is  thus  proportional  to 
M.  Accuracies  of  1%  of  range  appear  achievable. 
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APPENDIX  C 


ULTRASONICS  SYMPOSIUM 


NONINTRUSIVE  DYNAMIC  FLOWMETER 

N.  E.  Pedersen  and  L.  C.  Lynnworth 
Panametrics,  Inc.  ,  Waltham.  Massachusetts  02154 

ABSTRACT.  An  ultrasonic  dynamic  flowmeter  which  combines  nonintrusiveness. 
fast  response,  high  accuracy  and  high  resolution,  was  designed  for  use  with  cryogenic 
liquids  and  water.  The  flow  meter  was  designed  to  measure  to  1%  accuracy,  the  dynamic 
(1  to  100  Hz)  as  well  as  the  steady  flow  velocity,  averaged  over  the  pipe  area.  The  non- 
intrusive  requirement  is  satisfied  by  using  the  vertically  polarized  component  of  obliquely 
incident  shear  waves,  generated  by  transducers  mounted  outside  the  conduit.  Compared 
to  obliquely  incident  longitudinal  waves,  the  choice  of  obliquely  incident  shear  waves  typ¬ 
ically  doubles  the  solid/liquid  energy  transmission  coefficient,  and  doubles  the  refracted 
angle  in  the  liquid.  Calculations  showed  that  one  path  was  sufficient  to  obtain  the  ca - 

averaged  measurement  of  dynamic  flow  velocity  to  1%  accuracy.  The  electror.c  measuring 
system  employs  two  different  frequencies,  each  near  1  MHz.  These  frequenc5  -  arc  trans¬ 
mitted  continuously,  one  upstream  and  the  other  downstream,  over  a  common  pa'.h.  The 
frequencies  arc  servoed  such  that  the  phase  in  each  direction  is  maintained  cont-ant.  By 
means  of  subtracting  the  two  voltages  corresponding  to  the  two  frequency  deviations,  an 
output  signal  is  produced  which  Is  proportional  only  to  fluid  velocity  and  is  independent  of 
variations  in  sound  speed.  ' 


lnt  roduction 

According  to  references  1  and  2,  many 
liquid-propellant  rocket  vehicles  have  experienced 
longitudinal  vibration  because  of  an  instability  aris¬ 
ing  from  interaction  of  the  vehicle  structure  with 
the  propulsion  system.  The  vibrations,  nicknamed 
"pogo"  after  the  jumping  stiek,  have  occurred  prin¬ 
cipally  in  the  first  longitudinal  structural  mode  dur¬ 
ing  operation  of  the  first  liquid-propellant  stage  of  a 
launch  vehicle.  Figure  1  Illustrates  a  generalized 
representation  of  the  effects  of  a  typical  instability. 
The  vibration  begins  spontaneously,  intensifies,  and 
then  dies  away  -  typically  in  a  period  of  10  to  40 
seconds.  The  frequency  of  vibration  tracks  that  of 
the  first  structural  mode.  The  frequency  increases 
as  propellant  decreases.  Less  often,  pogo  vibra¬ 
tions  have  occurred  in  higher  modes  of  longitudinal 
vibration,  also,  multiple  periods  of  instability, 
each  involving  a  different  mode  of  vibration,  have 
occurred  during  operation  of  a  single  vehicle  stage. 
Vibrations  have  occurred  in  the  frequency  range  of 
5  to  60  Hz,  with  vibration  amplitudes  (zero  to  peak) 
as  high  as  17  g's  at  the  input  to  payloads  and  34  g' s 
at  an  engine. 


The  brief  account  given  below  presents  the  ultrasonic 
and  electronic  main  design  considerations.  Details 
are  available  in  another  report,  NASA  CR-112313.3 

Principal  Design  Requirements  (LOX  Region  1) 

1.  Nonintrusive  -  no  measurable  perturbation  of 
flow  profile  by  transducers. 

2.  Frequency  response:  1  to  100  Hz. 

3.  Accuracy:  errors  <  1%  for  area-average. 

4.  Resolution:  1%  of  oscillatory  flow  component,  or 

3  mm/ a,  whichever  is  greater 
(steady  flow  velocity  <  5.  5  m/s). 

5.  Application:  LOX,  ~93K,  J.  76  MN/m2  (~17  atm), 

pipe  ID  =  30  cm. 

Transducer  Mode  and  Orientation 

In  this  section  we  consider  the  selection  of 
the  transducer  vibration  mode  (longitudinal  vs  shear) 
and  transducer  orientation. 
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Fig.  1.  Typical  occurrence  of  pogo  vibration. 


Refraction  angles  are  readily  computed  from 
Snell's  Law,  provided  sound  speeds  are  known 
(Fig.  2).  Historically,  the  problem  of  calculating 
sound  energy  transmission  coefficients  for  solid/ 
liquid  interfaces  and  oblique  incidence  may  be  traced 
to  Knott.4  F rom  the  work  of  Firestone  and  Ling,5 
Mayer^  and  Garber,^  it  is  seen  that  obliquely  incident 
shear  waves  often,  but  not  always,  have  an  advantage 
regarding  efficiency  of  energy  transfer.  In  the  cases 
cited,  shear  offers  the  transmission  coefficient  ad¬ 
vantage  at  oblique  incidence  when  the  solid  has  a 
higher  characteristic  impedance  than  the  liquid. 

Using  the  Fortran  IV  G  "ERRWAB"  program 
listing  provided  to  us  as  a  courtesy  by  Garber  and 
Union  Carbide  Corporation,  Nuclear  Division,  Oak 
Ridge,  we  computed  the  energy  ratios  for  LOX/ 
aluminum  and  LOX/steel  interfaces. 


1973  Ultrasonics  Symposium  Proceedings, 
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DENSITY,  g/cm3 


Fig.  2.  Velocity  of  sound  plotted  against  density 
for  saturated  liquid  and  compressed  fluid  oxygen, 
courtesy  G.  C.  Straty  and  B.  A.  Younglove,  NBS. 
J.  Chem.  Thermodynamics  5^  305  (1973). 

I 

For  the  present  design  purposes  it  is  suffi¬ 
cient  to  plot  transmission  coefficients  for  an  L  wave 
incident  from  the  liquid  onto  the  solid  (Fig.  3).  One 
need  only  compare  these  energy  curves  to  see  that 
an  incident  shear  mode  would  be  coupled  twice  as 
efficiently  as  longitudinal.  Further,  comparing  the 
LOX /pipe  transmitted  energy  ratios  for  shear  inci¬ 
dent  at  60  ,  the  choice  of  aluminum  ovsr  stainless 
steel  yields  an  advantage  of  4.  3  dB  per  side.  This 
advantage  in  signal,  together  with  an  observed 
"acoustic  short  circuit"  noise  advantage  of  over  20 
dB,  strongly  favor  the  choice  of  an  aluminum  pipe 
over  steel.  (See  also:  LCL  1971  Patent  #3,575,050.  ) 

In  the  present  design,  it  turned  out  to  be  con¬ 
venient  to  separate  the  transmitter  and  receiver 
transducers  on  each  side  of  the  common  path,  as 
shown  in  Fig.  4.  For  oblique  incidence,  we  again 
can  make  a  comparison  of  shear  versus  longitudinal. 
Both  modes  share  the  advantage  of  avoiding  the  driv¬ 
ing  signal  (except  via  stray  capacitive  pick-up)  on 
the  receiving  transducer.  However,  if  the  shear 
wave  elements  are  oriented  at  t45  and  -45  to  the 
plane  of  Incidence,  such  that  neighboring  elements 
arc  vibrating  orthogonally  to  one  another,  then  theo¬ 
retically  there  is  no  acoustic  coupling  between  them. 


Fig,  3.  LOX/ aluminum  energy  coefficients. 


Fig.  4.  Pair  of  orthogonally  oriented  Bhear  wave 
transducers,  on  each  aide  of  fluid,  may  be  coupled 
to  fluid  under  test  acoustically  but  separated  elec¬ 
trically  from  each  other  by  the  sandwiched  insulator. 

Dynamic  Flow  Profile  Analysis 

The  purpose  of  this  analysis  was  to  deter¬ 
mine  the  value  of  the  correction  factor  K  required  to 
convert  an  ultrasonic  measurement  averaged  over  a 
single  path  whose  projection  coincides  with  a  diam¬ 
eter,  to  the  desired  area-averaged  flow  velocity. 

Thus.  K  vdia  e  varea.  The  analysis  was  conducted 
by  Professor  Forbes  T.  Brown  of  Lehigh  University. 
While  the  "static  K"  is  readily  found  in  terms  of 
Reynold's  number  to  be  about  0.  95  at  Re  =  10  ,  com¬ 
putation  of  the  dynamic  K  is  mors  complicated  and 
depends  on  frequency.®  For  the  Region  1  LOX  condi¬ 
tions,  and  subject  to  assumptions  stated  in  Ref.  3,  it 
was  found  that  K  >  0.  99,  1-100  He,  and  K  >  0.  995 
for  10-100  He.  Thus,  the  dynamic  profile  is  signif¬ 
icantly  flatter  than  the  steady-state  profile.  More¬ 
over  the  area-averaged  dynamic  flow  velocity  can 
be  taken  as  0.  995  vdia  with.errors  not  exceeding  0.5%. 
over  the  dynamic  frequency  range  of  interest,  1-100 Hz. 

Flow  Cell  Design 

Based  on  electronic,  ultrasonic  and  profile 
considerations,  the  cell  design  was  configured  as  in 
Fig.  5.  Measurements  are  made  simultaneously 
over  a  common  path.  Four  nonintrusive  shear  wave 
transducers  are  used,  oriented  _t  45  from  the  plane 
of  incidence,  with  60°  angles  of  incidence  obtained 
by  the  wedges. 


Electronic  System  Design 


To  obtain  a  voltage  output  which  is  indepen¬ 
dent  of  sound  speed  but  is  proportional  to  fluid  ve¬ 
locity,  one  can  set  the  downstream  and  upstream 
phases  equal  to  different  integer  constants: 


1  1 

=  <2N-»f- 

(1) 

sin  9  c  +  v  cos  8 

2-f  D 

=  (2M-l)f-. 

1  1 

(2) 

sin  6  c  -  v  cos  8 
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Solving  for  f  and  f^> 


fl  =  Tl^  (Z  N_1  >  (c  +  v  coa  6  ) 


f2  =  (2  M-l)  (c-  vcosfl). 


The  frequencies  fj  and  f2  are  each  composed  of  a 
fixed  and  variable  component: 


£10+F1:£2 


£20  +  F2 


where  fjg  and  f2g  are  preset  and  fixed,  and  Fj  and 
F2  are  voltage  controlled  variables.  This  set  of 
conditions  is  obtainable  in  a  number  of  ways,  e.  g.  , 
through  the  use  of  a  pair  of  Model  31  00A  Monsanto 
Digital  Frequency  Synthesizers.  By  means  of  the 
circuit  of  Fig.  6,  the  frequency  F,  is  controlled,  via 
a  first  feedback  loop,  to  maintain  the  condition  given 
by  Eq.  (1),  and  the  frequency  F2  is  controlled,  via  a 
second  feedback  loop,  to  maintain  the  condition  givsn 
by  Eq.  (2).  3  We  have,  then: 

F1  *  ^4^  (2  N-l)  (c  +  v  cos  B  )  -  f,0  (6) 

f2  =  *^0  (2  M-l)  (c  -  v  cosfl)  -  f2|)  (?) 

I  MHz  THICK  _ _ 


The  frequencies  F.  and  are  linear!/ 
proportional  to  thsir  respective  control  voltages  V 
and  V2.  respectively  Thus  F^  c  and  F^  = 

V_  where  K  and  are  front  panel  adjustments  on 
the  Monsanto  3I00A  synthesizers.  If  we  adjust  the  re¬ 
lationship  between  K  and  such  that  K  =  K  (2M-  1)/ 
(2N-1),  we  can  rewrite  (6)  and  (7)  in  terms  of  V  and 
V  ,  respectively: 

2  f 

V1  =  K1Td(2N'11  (ctvco»«)--^  (8) 


V  =  — 
2  K. 


(2  N- 1)  (c  -  v  cos  9  ) 


f20  /  2  N-  l\ 
"Kj  \  2M-  1 }' 


<V1-V2>  = 


Finally,  subtracting  (9)  from  (8),  we  obtain 
2  N-l  sinfl  cosfl  1  f,  (Z  N- 1\  ,  1 


v41[fio-(r^i)  £2o]<«>> 


Transducers  mounted  on  pipe  section. 


With  regard  to  the  term  on  the  RHS  of  (10) 
we  note  that  f  and  f  are  constants  to  within  10"° 
in  the  Monsanto  units;  M  and  N  are  invariants  as  isKj. 
Therefore,  this  term  is  a  constant  which  can  be  sub¬ 
tracted  from  under  zero  flow  conditions  in  the 

case  of  measuring  the  steady  state  flow  component,  or 
neglected  entirely  if  only  a  dynamic  flow  measurement 
is  to  be  made.  Variations  in  sound  speed  can  also  be 
measured  if  desired  by  adding  Eqs.  (8)  and  (9): 

The  reason  for  providing  a  significant  (~109J 
separation  between  the  frequencies  of  the  upstream 
and  downstream  waves  is  that  such  a  separation  eli¬ 
minates  the  need  for  a  phase  code  modulatlon/demo- 
dulation  scheme  as  is  summarized  in  ref.  9.  More 
conventional  techniques  suffice  to  isolate  the  received 
and  transmitted  waves  at  the  transducer  locations. 


Although  the  contract  acknowledged  below 
called  for  a  design  only,  not  tssts,  it  was  agreed  to 
conduct  a  limited  experimental  program.  The  most 
significant  experimental  result  was  obtained  using  a 
four  transducer  arrangement  connected  electrically 
as  in  Fig.  6,  and  mounted  in  a  cell  similar  to  that 
used  in  an  earlier  program.^  This  test  demonstrated 
stable,  independent  and  simultaneous  frequency 
servoing  of  both  the  upstTeam  and  downstream  waves. 
Also,  upon  repetitive  intentional  interruption  of 
either  received  signal,  the  system  immediately 
locked  onto  the  same  total  phase  when  the  signal 
was  reconnected.  By  filtering,  one  can  separately 
read  out  both  the  dynamic  (1-100  Hz)  and  dc  flow 
components.  These  experimental  results  prove  that 
a  servoed  frequency  system  can  in  fact  be  imple¬ 
mented.  The  number  and  type  of  constraints  imposed 
upon  the  phases  can  be  varied  somewhat  depending  on 
the  details  of  the  application.  For  example,  in  an¬ 
other  program,10  constraints  are  chosen  to  provide 
an  output  voltage  proportional  to  v/c,  the  Mach  num¬ 
ber.  This  number,  when  multiplied  by  the  fluid's 
impedance  p  c,  yields  a  final  output  proportional  to 
mass  flow  rate.  In  general,  one  can  arrange  for  the 
flow  "velocimeter"  to  yield  v,  v/c,  or  v/c^. 
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LEGEND:  PA  =  Power  Amp’lf'er;  B  -  Bu'fer  Amp'.ifer;  DC  =  DC  Amplifier;  IF  =  IF  Amplifier. 
Fig.  6.  Block  diagram  of  dynamic  flowmeter  electronic  deaign. 


We  are  now  building  a  complete  ayatem  for 
initial  dynamic  flow  tcating  on  liquid  nitrogen  at 
NBS-Bouider  (pipe  diameter  ~10  cm),  with  aubse- 
quent  testing  at  NASA-Huntsvillc  (pipe  diameter 
~30  cm).  Future  applicationa  may  Include  ground 
teata  for  the  apace  shuttle  Orbiter,  where  dynamic 
flow  velocity  needs  to  be  meaaured  in  liquid  oxygen 
and  hydrogen  at  various  locations  (NASI -1 1 756). 
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ISA-authorizcd  reprint  from:  posium,  conducted  1971  in  Pittsburgh,  PA.  Sponsored 

by  American  Society  of  Mechanical  Engineers,  Amcri- 
Ftow  -  Its  Measurement  and  Control  in  Science  uid  can  Institute  of  Physics,  National  Bureau  of  Standards, 
Industry,  Volume  1.  and  ISA. 

NONINTRUSIVE  ULTRASONIC  MEASUREMENT 
OF  FLOW  VELOCITY  AND  MASS  FLOW  RATE 

L.C.  LYNNWORTH,  N.E.  PEDERSEN,  E.P.  PAPADAKIS*  AND  J.H.  BRADSHAW** 

Panametrics 
Waltham,  Massachusetts 


Four  nonintrusive  methods  are  analyzed  for  introducing  ultrasonic  waves  into  fluids  within  a  pipe,  without 
mechanically  perturbing  the  inside  of  the  pipe.  In  one  method,  obliquely  incident  shear  waves  are  mode 
converted  to  longitudinal  waves  in  the  fluid.  In  a  second  method,  wavefronts  are  diffracted  into  fluids  at  an 
angle.  A  third  method  involves  beaming  plane  or  focused  waves  across  two  separated  paths,  and  measuring  the 
correlation  time  between  fluctuations  in  propagation  over  the  two  paths.  A  fourth  method  would  use  a  fixed 
transducer  array  to  control  the  sound  distribution  in  the  fluid.  Other  nonintrusive  methods  such  as  inserts, 
recesses  or  grooves  require  that  tha  pipe  be  penetrated  or  its  inside  perturbed.  Techniques  are  described  which 
weight  flow  velocity  proportional  to  cross-sectional  area.  Theoretically,  this  yields  acoustic  flow  velocity 
readings  which  essentially  do  not  require  profile  corrections  for  laminar,  transitional  or  turbulent  flaw.  One 
application  of  nonintrusive  ultrasonic  measurement  of  flow  velocity,  together  with  nonintrusive  ultrasonic 
density  compensation,  is  the  measurement  of  mass  flow  rate  M  of  liquid  fuels  in  a  helicopter  engine.  Here,  we 
demonstrated  a  transmission  technique  which  employed*  phase-coded  continuous  waves  and  a  redaction 
coefficient  densitometer  The  coded  waves  were  separated  and  their  phase  difference  was  measured  after 
upstream  and  downstream  oblique  transmission  across  square-channeled  flow,  over  a  common  path  between  a 
pair  of  recessed,  rectangularly-collimated  transducers. 


INTRODUCTION 

This  paper  is  concerned  with  the  analysis  and 
design  of  flowmeters  employing  nonintrusive 
ultrasonic  transducers  for  the  measurement  of 
fluid  flow,  sometimes  without  any  penetration 
or  perturbation  of  the  container  of  the  fluid. 
We  shall  limit  the  scope  to  analytic  and 
experimental  work  done  in  our  laboratory  these 
past  five  years,  and  shall  cover  only  the  main 
points. 

The  user  of  a  flowmeter  usually  seeks  a  precise, 
and  sometimes  accurate,  measurement  averaged 
over  the  full  cross-sectional  area  of  the  fluid, 
and  averaged  over  a  tune  interval  which  may 
range  from  10  ms  in  cases  of  aircraft  engine 
control1  or  Pogo  (flow  oscillation)  diag¬ 
nostics,3  to  24  hr  in  the  case  of  pollution 
monitoring  of  stack  exhausts.3  There  are  also 
medical  applications  where  flow  velocity  profile 
is  to  be  determined.4  and  wind  tunnel  applica¬ 
tions  where  profiles  of  both  flow  velocity  and 
density  are  to  be  measured,  to  determine  mass 
flow  rates.5  Accuracy  requirements  are  usually 
in  the  0.1  to  10  percent  range.  Mach  numbers 

'Present  address,  ford  Motor  Co.,  Detroit, 
Michigan. 

"Present  address;  Gillette  Co.,  Boston,  Massa¬ 
chusetts. 


range  from  less  than  I0T3  to  greater  than  10; 
Reynolds  number  Re  ranges  from  less  than  100 
to  greater  than  107.  Diameters  range  from 
~l  mm  to  tens  of  meters.  Fluids  are  not  always 
single-phase  or  homogeneous,  and  flow  is  not 
always  unidirectional  or  steady. 

Ultrasonics  may  offer  the  preferred  approach  to 
problems  containing  much  of  this  complexity. 
One  of  the  first  details  to  be  considered  deals 
with  data  averaging,  both  area-  and  time-wise. 
Details  of  nonintrusiveness,  precision,  accuracy, 
reliability,  response  time,  equipment  versality, 
complexity,  maintainability,  cost  and  other 
factors  are  also  important. 

Average  Flow  Velocity 

The  ultrasonic  measurement  of  flow  velocity 
averaged  across  a  single  path  which  is  a  small 
fraction  of  the  circular  pipe  area  usually 
exceeds  the  true  area-averaged  flow  velocity  v 
by  less  than  8  percent  for  Re  >  I04 ,  but  by  25 
percent  for  laminar  flow.1’*’7  We  have  found 
that  one  way  to  properly  average  the  flow 
velocity  profile  is  to  measure  transmissions 
across  the  entire  fluid  cross  section  in  a 
rectangular  or  square  flow  channel,  using  rec¬ 
tangularly-  or  square-collimated  plane  waves.1 
Another  way  is  to  scatter  a  uniform  plane  wave 
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off  scatterers  which  are  uniformly  distributed,  OBLIQUELY  INCIDENT 

and  measure  the  doppler  shift.1  These  methods  SHEAR  WAVES 

avoid  assumptions  about  the  flow  profile.  In 

some  cases  of  unsteady  flow,  the  major  part  obliquely  incident  vertically  polarized  ultra- 
of  the  flow, has  a  profile  so  nearly  flat  that  sonjc  shear  waves  can  be  mode-converted  to 
measurements  over^a  single  path  should  provide  longitudinal  waves,  and  then  transmitted 
I  percent  accuracy.  diagonally  across  the  conduit  such  that  there  is 

To  achieve  the  highest  accuracy  in  a  given  time  a  substantial  component  of  the  longitudinal 
interval,  continuous  waves  potentially  are  best,  velocity  in  a  direction  parallel  to  the  axis  of 
Various  electronic  methods  such  as  phase-code  flow  of  the  fluid.9  When  this  is  achieved  with 
modulation,  dual  frequencies  and  fm  have  now  transducers  mounted  external  to  the  conduit 
been  demonstrated,  to  separate  upstream  from  there  is  no  perturbation  of  the  flow  of  the  fluid 
downstream  waves  transmitted  simultaneously  stream,  other  than  possibly  thermal 
over  an  essentially  common  path  between  a  perturbation, 
given  pair  of  transducers.1  •*  >® 

In  any  particular  application,  the  optimum  demonstrated  this  noninvasive  technique  in 
electronic  approach  generally  cannot  be  water  contained  in  a  Cu  pipe  ~  1,3  cm  diam- 
selected  independent  of  the  details  of  the  eler-  and  in  a  steel  pipe  2S  cm  diameter,  in  the 
transducer  mounting  material  and  geometry  -  ,ate  sixties.  Present  applications  involve  cryo- 
and  vice  versa.  The  selection  also  must  rccog-  genic  liquids,  water  and  liquid  sodium, 
nize  that  different  electronic  flow  velocimeter 

approaches  respond  primarily  to  v,  v/c  or  v/cl .  R  ,  compares  paths  for  longitudinal  and 
where  c  is  the  sound  speed  in  the  flurd.  shear  waves  ind(Jem  „  4S  deg  jn  waler  Note 

We  next  present  four  nonintrusive  transducer  that  the  probes  are  flush-mounted  inserts.  If  the 
designs  which,  to  our  knowledge,  originated  in  fluid  is  water,  and  the  insert  probe  (wedge)  is 
our  laboratory.  made  of  acrylic,  (a)  Ll  =  0.44D,  and  (b) 


Fig.  1.  Comparison  of  obliquely  incident  (45  degrees)  longitudinal  and  shear  waves. 
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Ls  =  5.67D,  at  room  temperature.  In  (c)  the 
focused  incident  wave  permits  the  axial- 
projccted  path  to  remain  constant  by  means  of 
temperature-adaptive  ray  selection.  In  (d)  we 
show  a  hybrid  case.  For  some  combinations  of 
liquids  and  pipe  wall  materials,  the  45  deg 
obliquely-incident  longitudinal  wave  in  the  still 
liquid  is  refracted  and  mode-converted  to  a 
shear  wave  in  the  pipe  wall.  This  shear  wave 
reconverts  to  a  single  longitudinal  wave  in  the 
flowing  liquid.  When  the  sound  speeds  are  equal 
in  the  liquids  in  the  transducer  cavity  and  in  the 
pipe,  L=  D. 

It  can  be  seen  that  the  shear  waves  are 
mode-converted  by  refraction  at  the  interface 
into  longitudinal  waves,  which  are  now  trans¬ 
mitted  at  a  much  higher  angle,  typically  twice 
that  achievable  with  an  incident  longitudinal 
wave,  through  the  fluid  material.  Since  (he 
ultrasonic  waves  being  transmitted  through  the 
fluid  have  a  significant  component  of  velocity 
in  a  direction  parallel  to  the  axis  of  flow,  the 
resolution  of  the  measurement  is  increased  by  a 
factor  of  two  or  more.  Further,  obliquely 
incident  shear,  compared  to  obliquely  incident 
longitudinal,  may  be  coupled  more  efficiently 
into  the  fluid.10  In  the  pipe  itself,  the  slower, 
zigzagging  shear  wave  is  less  troublesome  than  a 
zigzagging  longitudinal  and  its  mode-converting 
reverberations.  At  high  temperature,  the  avail¬ 
able  transducer  material  (e.g.,  LiNbOj)  may  be 
more  efficient  in  shear  than  in  longitudinal 
mode.  Another  advantage  of  using  shear  waves 
obliquely  incident  at  a  solid/liquid  interface  is 
that  upstream  and  downstream  waves  can  be 
generated  and  received  by  pairs  of  transducers 
polarized  orthogonally  at  ±45  degrees  to  the 
plane  of  incidence.  This  facilitates  separation  of 
the  upstream  and  downstream  waves  even  if 
they  are  transmitted  continuously  over  a 
common  path,  at  the  same  frequency. 

DIFFRACTION  GRATING 

A  second  nonintrusive  method  of  introducing 
ultrasound  into  the  pipe  (duct,  channel)  and  to 
detect  it  is  by  means  of  ultrasonic  diffraction 
gratings.13  This  method,  used  in  dispersive 
delay  lines,  is  advantageous  from  the  point  of 
view  of  sensitivity  in  one  measurement  scheme 
because  the  flow  distorts  and  translates  the 
Huygens  wavelets  radiated  into  the  fluid  in  such 
a  way  as  to  steer  the  diffracted  beams.  Consider 


the  grating  in  Fig.  2.  The  circular  Huygens 
wavelets  in  the  no-flow  condition  are  warped 
by  laminar  flow  or  translated  (with  distortion) 
by  turbulent  flow.  The  +1  diffracted  order, 
radiating  downstream  as  drawn,  has  its  Poynting 
vector  P  steered  more  downstream,  making 
its  reception  at  the  opposite  side  of  the  channel 
occur  later  than  it  would  have  under  no-flow 
conditions.  This  result  is  opposite  to  the 
ordinary  result  in  which  the  flow  velocity 
causes  the  wave  to  arrive  earlier  in  time.  When 
angle  <t>  is  large  (approaching  90  degrees), 
the  steering  effect  can  be  large. 

Reception  is  accomplished  by  a  second  diffrac¬ 
tion  grating,  longer  in  extent  and  of  the  same 
spacing.  The  diffraction  grating  flowmeter 
could  be  operated  in  a  pulsed  rf  mode  with 
time-of- flight  as  the  measured  parameter. 
Another  mode  of  operation  is  fm.  The  beam 
can  be  steered  from  the  transmitter  grating  to  a 


HuyE'"*  for  no  flow  The  Poynting  vector  P 

ta  generated  at  angle  v  defusing  ike  path  ef  energy  flow 


Waveleta  for  laminar  flow 


Wavelata  (or  turbulent  now 

1'ig.  2.  Analysis  of  wavelets  launched  by  diffraction 
grating. 
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particular  location  on  the  receiver  grating  (or  to 
a  very  small  receiver)  by  changing  the  input 
frequency  to  compensate  for  the  shifting  of  the 
Huygens  wavelets  by  the  flow.  Readout  would 
be  the  frequency  generating  the  maximum 
received  signal. 

For  thin-walled  pipe,  tubing,  and  ducts,  the 
grating  may  be  on  the  outside  of  the  wall.  Here, 
“thin-walled”  is  less  than  a  wavelength.  For 
thick-walled  vessels,  the  grating  should  be  on 
the  inside  (this  may  not  always  be  “nonin- 
trusive”).  Screw  threads  are  adequate  gratings 
for  transmission  and  reception  in  circular  pips. 
While  the  grating  diffraction  method  of  intro¬ 
ducing  a  wave  into  a  pipe  (duct,  channel)  can 
be  utilized  for  beam  steering  in  the  particular 
flowmeter  scheme  described  above,  it  can  also 
be  used  to  introduce  ultrasonic  waves  non- 
intrusively  for  other  purposes.  An  endfire  array 
(0  =  90  degrees)  can  be  used  for  ordinary 
time-of-flight  measurements  or  for  doppler 
measurements.  A  principal  limitation  on  this 
method  is  its  extreme  sensitivity  to  the  bound¬ 
ary  layer  right  next  to  the  wall. 

NORMALLY  INCIDENT 
LONGITUDINAL  WAVES 

This  nonintrusive  method  involves  beaming 
plane  or  focused  waves  across  two  axially- 
separated  diameters.16  Flow  transit  lime 
between  these  two  diameters  is  taken  as  the 
correlation  time  between  fluctuations  in  propa¬ 
gation  (velocity,  scattering  and/or  attenuation) 
over  said  diameters.1 7 

The  method  is  illustrated  schematically  in 
Fig.  3.  Provided  the  fluid  wets  or  contacts  the 
pipe  with  sufficient  intimacy  so  ultrasound  can 
be  coupled  across  the  interface,  the  acoustic 
part  is  relatively  simple.  One  chooses  a  fre¬ 
quency  high  enough  to  yield  directivity,  and  to 
attenuate  the  potentially  troublesome  acoustic 
short  circuit  around  the  pipe  perimeter,  yet  not 
so  high  that  attenuation  in  the  fluid  itself  is 
excessive. 

To  reduce  memory  requirements  in  the  correla¬ 
tion  detector,  one  can  periodically  introduce  a 
thermal  or  other  disturbance.  Naturally,  within 
the  spirit  of  the  present  “nonintrusive” 
approach,  we  suggest  this  perturbation  be  intro¬ 
duced  noninvasively.  For  example,  a  high- 


Fig.  3.  Correlation  lime  between  longitudinal 
waves  provides  a  measure  of  flow  velocity. 


intensity  focused  pulse  introduced  upstream 
will  cause  a  slight  temperature  rise,  and  possibly 
cavitation  near  the  axis. 

As  with  other  tagging  techniques  or  correlation 
flowmeters  reported  in  this  Symposium,  care 
must  be  exercised  with  regard  to  weighting  and 
tag  smearing. 

FIXED  TRANSDUCER  ARRAY 

In  this  method  of  nonintrusive  flow  measure- 
ment,.sound  would  be  introduced  into  a  section 
of  pipe  and  detected  by  an  array  of  Fixed 
transducers.1 8  The  acoustic  wavelength  and 
pipe  geometry  determine  the  resultant  mode. 
This  mode  ideally  would  weight  each  area 
element  of  the  fluid’s  cross  section  equally.  The 
array  requires  that  a  fixed  acoustic  wavelength, 
X,  be  maintained  within  the  pipe.  The  average 
velocity  of  the  flow  within  the  pipe  is  deter¬ 
mined  by  measuring  the  frequency  shift  f  ±  F 
required  to  maintain  X  (and  applying  profile 
corrections).  Since  the  shifted  frequency 
f  =  (c  ±  v)/X  due  to  flow  velocity  can  be 
measured  in  either  direction,  the  difference 
between  the  two  can  be  extracted  directly  and 
is  equal  to  Af  =  (c  +  v)/X  -  (c  —  v)/X  =  2  v/X. 
The  array  is  composed  of  a  row  of  transducers 
spaced  and  polarized  such  that  the  response  due 
to  sound  “short-circuiting”  through  the  walls  of 
the  pipe  will  be  negligible  around  the  operating 
frequency.  The  last  transducer  on  each  end  is 
used  for  the  transmitter  of  the  up-  and  down¬ 
stream  sound  artd  a  sing-around  oscillator  is 
provided  for  each  direction.  The  two  fre¬ 
quencies  are  mixed  to  obtain  the  difference,  Af. 
Figure  4  illustrates  this  system  schematically. 

There  are  two  conditions  which  should  be 
considered  when  using  this  technique  of  flow 
measurement.  First,  the  effective  sound 
velocity  in  the  fluid  (cf)  must  be  different  from 
that  in  the  pipe  (Cp).  That  is,  if  2n  is  the 


121 


AEDC-TR-74-77 


NONINTRUSIVE  ULTRASONIC  MEASUREMENT 


Fig.  4.  Transducer  array  and  sing-around 
flowmeter. 

number  of  receive  transducers  and  cp/cf  =  n'/n 
then  n'  should  not  be  harmonically  related  to  n 
(i.e.,  ri  ^  mn)  within  the  range  cf  =  c  ±vmax. 
Secondly,  the  flow  profile  should  be  reasonably 
defined,  consistent  with  accuracy  requirements. 
(The  fluid  itself  need  not  be  homogeneous.) 
The  type  of  profile  correction  is  somewhat 
flexible,  as  the  choice  of  acoustic  mode 
(depends  on  pipe  diameter  to  wavelength)  or 
other  considerations  may  make  this  correction 
unnecessary,  or  less  critical  than  in  other 
sing-around  flowmeters  which  dp  not  sample 
the  entire  fluid  cross  section.  The  authors  have 
not  yet  demonstrated  this  method. 

RECESSED  TRANSDUCERS 

Recessed,  nonintrusive  transducers  have  been 
used  by  several  earlier  investigators.  In  our 
work  on  the  helicopter  engine  fuel  mass  flow¬ 
meter,  we  used  this  approach,  but  not  exactly 
in  its  usual  form.  Rather,  we  recessed  the 
transducers  behind  a  rectangular  “sound  port”, 
to  achieve  a  uniform  weighting  of  the  square- 
channeled  flow  intersected  by  the 
rectangularly-collimated  sound  beam.1 

The  hybrid  designs  of  Fig.  5  retain  the  advan¬ 
tages  of  a  transducer  mounted  outside  the 
conduit,  yet  retain  a  beam  direction  which,  like 
in  the  case  of  the  usual  recessed  transducer,  can 
be  substantially  or  completely  independent  of 
sound  speed  in  the  fluid. 

MISCELLANEOUS  INTERROGATION 
METHODS 


b) 


Fig.  5.  Use  of  vee  grooves  to  transmit  over  oblique 
paths. 

fluid  flow  through  permeable  reflectors 
(Fig.  6c),n  or  between  heat  exchanger  tubes 
offering  limited  access  (Fig.  6d),  are  examples 
of  special  applications  where  ultrasonics  offers 
a  unique  approach.  Figure  7  illustrates  a  more 
general  potential  application,  where  ultrasound 
could  measure  v  plus  one  or  more  other  process 
variables.1 3 
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Fig.  6a.  Approximate  measure  of  subsonic  gas  flow 
velocity  obtained  at  combustor  exit  by  ratio  of  probe 
offset  to  gas  path  length. 
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Ultrasonic  measurement  of  Mach  number 
(Fig.  6a),1 1  of  the  drawing  speed  of  wire, 
particularly  magnetostrictive  wire  (Fig.  6b),  of 


Fig.  6b.  Measurement  of  material  velocity  using 
guided  exiensional  waves.  Waves  are  introduced  and 
detected  magnetostrictively,  without  contact  between 
coil  and  material. 
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Fig.  6c.  Reflectors  for  bidirectional  interrogation 
of  flow  region  using  remote  transducers. 


Fig.  6d.  Correlation  lime  between  transmissions 
beamed  over  two  parallel  paths  a  distance  L  apart, 
provides  a  measure  of  flow.  Mulli-bcvclcd  probes  arc 
shown  located  in  triangular  cusps.  Analogously,  a 
single  wire  or  rod  may  contain  temperature-sensing 
zones  spaced  L  apart.  One  would  measure  correlation 
lime  between  1|  and  Ij. 

MASS  FLOWMETER  FOR 
AVIATION  FUELS 

From  the  viewpoint  of  those  responsible  for 
diagnosing  and  controlling  the  performance  of 
aircraft  engines,  it  is  important  to  know 
accurately  the  value  of  the  fuel  mass  flow  rate 
M.1  In  1972,  we  demonstrated  a  fuel  flowmeter 
consisting  of  a  flow  velocimeter,  densitometer, 
time  intervalometer,  and  a  metering  section 


containing  nonintrusive  transducers.  The 
velocimeter  measured  the  phase  difference 
between  two  coherent  individually  coded  5 
MHz  waves  continuously  transmitted  obliquely 
upstream  and  downstream  across  a  common 
path,  yielding  a  phase  difference  proportional 
to  v/cJ.  The  densitometer  consisted  of  a 
stepped-diampter  probe;  the  difference  in  echo 
amplitudes  at  the  wetted  end  and  at  the  dry 
step  is  a  function  of  the  acoustic  impedance  pc 
of  the  fuel  where  p  is  the  fuel  density.  The  time 
intervalometer  measured  the  time  r  between 
two  successive  echoes,  which  time  is 
proportional  to  1/c.  The  mass  flow  rate  is 
M  =(K)(v/c2)  (pc)(l/r)  where  K  is  a  constant. 
The  complete  system  was  tested  on  various 
liquids  at  rates  up  to  ~2300  kg/hr.  It  was 
operated  during  and  recalibrated  after  104-hr 
exposure  to  a  contaminated  fluid  flowing  at 
~860  kg/hr.  Response  time  was  20  ms.  The 
flowmeter  can  operate  in  laminar,  transitional 
and  turbulent  flow.  It  was  concluded  that 
refinement  of  this  approach  should  be  able  to 
attain  M  accuracy  of  1  percent,  despite  a  wide 
range  of  Reynolds  numbers  and  other 
variables.1  Work  directed  to  attaining  this 
performance  is  now  in  progress. 

CONCLUSIONS 

Nonintrusive  flow  measurements  are  attractive 
because  of  their  potential  application  to  a 
variety  of  flowmetry  problems  where  it  is 
inconvenient  or  undesirable  to  penetrate  the 
container,  or  to  locate  a  sensor  in  the  flow 
itself.  A  variety  of  ultrasonic  transducer  designs 
are  now  available  which  do  not  obslruct  the 
flow.  In  some  cases  the  transducers  may  be 
mounted  outside  the  conduit,  so  that  they  do 
not  disturb  the  flow  profile  in  any  significant 
way. 

Depending  on  the  details  of  the  application,  a 
transmission  or  a  scattering  approach  may  be 
preferred.  It  may  even  be  desirable  to  incor¬ 
porate  both  approaches  in  one  modular  system, 
to  handle  two-phase  applications  where  scatter¬ 
ing  is  not  always  sufficient  to  rely  on  a  doppler 
system,  yet  may  sometimes  be  so  large  as  to 
interfere  with  a  transmission  measurement. 
Other  factors  such  as  area-averaging,  range/ 
resolution  selection  by  frequency  shift  or  multi¬ 
bounce  transmission,  accuracy  and  response 
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Control  Signals  to  Pumps,  Molars,  Mistrs, 
Sproys,  Hooters,  Tensioners,  Etc. 


Hg.  7.  Schematic  of  ultrasonic  instrumentation  system  for  fluid  process  measurement  and  control. 


time  requirements  enter  into  the  optimization 
of  available  flowmeter  approaches  for  specific 
applications. 

Analysis  of  specific  applications  such  as 
measuring  the  steady  flow  of  aviation  fuels,  and 
the  oscillatory  flow  of  cryogenic  liquids,  to 
absolute  accuracies  of  I  percent  or  better  with 
a  response  time  of  about  10  ms,  showed  that,  at 
least  for  these  cases,  the  optimum  approach 


should  include  the  use  of  cw  or  rf  bursts,  not 
pulses,  transmitted  upstream  and  downstream 
over  a  common  path. 

At  this  occasion  of  the  First  Symposium  on 
Flow,  it  is  instructive  to  look  back  ten  years  to 
a  review  which  compared  ultrasonic  and  other 
flowmeters  with  regard  to  maximum  flow  rate, 
accuracy,  media,  temperature,  maximum 
pressure  and  viscosity.1 4  The  past  decade  has 
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witnessed  experiments  on  a  number  of  ultra¬ 
sonic  flowmeters,  including  successful  applica¬ 
tions  in  gases,  liquids  and  airborne  solids, 
primarily  for  industrial  process  control. 
Looking  ahead  to  the  seventies,  we  may 
anticipate  ultrasonics’  expanding  role  in  human 
and  ecological  terms,  based  on  measuring  flow 
not  only  in  a  variety  of  industrial  applications, 
but  also  in  fluids  in  the  human  body  and  in  the 
ecosphere. 
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LIST  OF  SYMBOLS 

a  absorption  coefficient  for  dry  air;  angle  between  direction  of 

propagation  and  normal  to  phase  front 

a  classical  sound  absorption  coefficient  due  to  thermal  conduc¬ 

tivity  and  viscosity 

a  angle  of  wave  propagation  relative  to  normal  to  transducer  face 

y  specific  heat  ratio 

6*  displacement  thickness  of  turbulent  boundary  layer 

rj  viscosity 

0  angle 

X  wavelength 

|x  magnetic  permeability 

|x  index  of  refraction  rms  fluctuation 

(  particle  displacement  in  sound  wave 

3 

p  density  gm/cm 

T  correlation  delay 

0/  time  correlation  function 

xy 

w  angular  frequency  (2  it  f)  of  sound  wave 

convection  frequency  characteristic  of  boundary  layer 
a  electrical  conductivity 

a^  correlation  length  of  burbulence 

B  velocity  correlation  coefficient 

B,j,  thermal  fluctuation  correlation  coefficient 

c  sound  speed 
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LIST  OF  SYMBOLS  (cont'd) 

2 

area  density^  gm/cm  ;  transducer  diameter 
duct  diameter 
frequency  of  sound  wave 

frequency  difference  between  two  sound  waves 

doppler  frequency 

wave  number  of  sound  wave  - 

\ 

pathlength,  length  of  shock  formation 
sound  path 

loss  due  to  acoustic  impedance  mismatch 
loss  due  to  absorption 
Mach  number 

pressure  fluctuation  in  sound  wave 

maximum  pressure  acoustic  wave  which  can  be  propagated 
a  distance  l  without  shock  formation 

Prandtl  number 

time 

Temperature 

mean  square  temperature  fluctuation  of  turbulence 

maximum  flow  velocity 

flow  velocity  (particle  velocity) 

mean  square  velocity  fluctuation  of  turbulence 

streamwise  distance  and  beam  drift  distance 

normal  distance  from  transducer  face 
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